Program: gdcm
Module: $RCSfile: gdcmUtil.cxx,v $
Language: C++
- Date: $Date: 2005/01/25 15:35:15 $
- Version: $Revision: 1.123 $
+ Date: $Date: 2005/04/19 09:58:19 $
+ Version: $Revision: 1.148 $
Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
l'Image). All rights reserved. See Doc/License.txt or
#include <stdio.h> //only included in implementation file
#if defined(_MSC_VER) || defined(__MINGW32__)
- #include <winsock.h> // for gethostname & gethostbyname
+ #include <winsock.h> // for gethostname and gethostbyname
#undef GetCurrentTime
#else
#ifndef __BORLANDC__
// For GetCurrentThreadID()
#ifdef __linux__
-#include <sys/types.h>
-#include <linux/unistd.h>
+ #include <sys/types.h>
+ #include <linux/unistd.h>
#endif
#ifdef __sun
-#include <thread.h>
+ #include <thread.h>
#endif
namespace gdcm
{
+//-------------------------------------------------------------------------
+const std::string Util::GDCM_UID = "1.2.826.0.1.3680043.2.1143";
+std::string Util::RootUID = GDCM_UID;
+const uint16_t Util::FMIV = 0x0001;
+uint8_t *Util::FileMetaInformationVersion = (uint8_t *)&FMIV;
+
+//-------------------------------------------------------------------------
+// Public
/**
* \brief Provide a better 'c++' approach for sprintf
* For example c code is:
- * sprintf(trash, "%04x|%04x", group , elem);
+ * char result[200]; // hope 200 is enough
+ * sprintf(result, "%04x|%04x", group , elem);
*
* c++ code is
* std::ostringstream buf;
* buf.str();
*
* gdcm style code is
- * Format("%04x|%04x", group , elem);
+ * string result;
+ * result = gdcm::Util::Format("%04x|%04x", group , elem);
*/
-
std::string Util::Format(const char *format, ...)
{
char buffer[2048];
/**
* \brief Because not available in C++ (?)
+ * @param str string to check
+ * @param tokens std::vector to receive the tokenized substrings
+ * @param delimiters string containing the character delimitors
+
*/
void Util::Tokenize (const std::string &str,
std::vector<std::string> &tokens,
- const std::string& delimiters)
+ const std::string &delimiters)
{
std::string::size_type lastPos = str.find_first_not_of(delimiters,0);
std::string::size_type pos = str.find_first_of (delimiters,lastPos);
/**
* \brief Because not available in C++ (?)
* Counts the number of occurences of a substring within a string
+ * @param str string to check
+ * @param subStr substring to count
*/
int Util::CountSubstring (const std::string &str,
const std::string &subStr)
{
- int count = 0; // counts how many times it appears
- std::string::size_type x = 0; // The index position in the string
+ int count = 0; // counts how many times it appears
+ std::string::size_type x = 0; // The index position in the string
do
{
- x = str.find(subStr,x); // Find the substring
- if (x != std::string::npos) // If present
+ x = str.find(subStr,x); // Find the substring
+ if (x != std::string::npos) // If present
{
- count++; // increase the count
- x += subStr.length(); // Skip this word
+ count++; // increase the count
+ x += subStr.length(); // Skip this word
}
}
- while (x != std::string::npos); // Carry on until not present
+ while (x != std::string::npos);// Carry on until not present
return count;
}
return r;
}
+unsigned int Util::GetCurrentThreadID()
+{
+// FIXME the implementation is far from complete
+#if defined(_MSC_VER) || defined(__BORLANDC__) || defined(__MINGW32__)
+ return (unsigned int)GetCurrentThreadId();
+#endif
+#ifdef __linux__
+ return 0;
+ // Doesn't work on fedora, but is in the man page...
+ //return (unsigned int)gettid();
+#endif
+#ifdef __sun
+ return (unsigned int)thr_self();
+#else
+ //default implementation
+ return 0;
+#endif
+}
+
+unsigned int Util::GetCurrentProcessID()
+{
+#if defined(_MSC_VER) || defined(__BORLANDC__) || defined(__MINGW32__)
+ // NOTE: There is also a _getpid()...
+ return (unsigned int)GetCurrentProcessId();
+#else
+ // get process identification, POSIX
+ return (unsigned int)getpid();
+#endif
+}
+
+/**
+ * \brief tells us if the processor we are working with is BigEndian or not
+ */
+bool Util::IsCurrentProcessorBigEndian()
+{
+#if defined(GDCM_WORDS_BIGENDIAN)
+ return true;
+#else
+ return false;
+#endif
+}
+
/**
* \brief Create a /DICOM/ string:
* It should a of even length (no odd length ever)
return s1_even == s2_even;
}
-/**
- * \brief tells us if the processor we are working with is BigEndian or not
- */
-bool Util::IsCurrentProcessorBigEndian()
-{
-#ifdef GDCM_WORDS_BIGENDIAN
- return true;
-#else
- return false;
-#endif
-}
-
-
-
#ifdef _WIN32
-typedef BOOL(WINAPI * pSnmpExtensionInit) (
- IN DWORD dwTimeZeroReference,
- OUT HANDLE * hPollForTrapEvent,
- OUT AsnObjectIdentifier * supportedView);
-
-typedef BOOL(WINAPI * pSnmpExtensionTrap) (
- OUT AsnObjectIdentifier * enterprise,
- OUT AsnInteger * genericTrap,
- OUT AsnInteger * specificTrap,
- OUT AsnTimeticks * timeStamp,
- OUT RFC1157VarBindList * variableBindings);
-
-typedef BOOL(WINAPI * pSnmpExtensionQuery) (
- IN BYTE requestType,
- IN OUT RFC1157VarBindList * variableBindings,
- OUT AsnInteger * errorStatus,
- OUT AsnInteger * errorIndex);
-
-typedef BOOL(WINAPI * pSnmpExtensionInitEx) (
- OUT AsnObjectIdentifier * supportedView);
+ typedef BOOL(WINAPI * pSnmpExtensionInit) (
+ IN DWORD dwTimeZeroReference,
+ OUT HANDLE * hPollForTrapEvent,
+ OUT AsnObjectIdentifier * supportedView);
+
+ typedef BOOL(WINAPI * pSnmpExtensionTrap) (
+ OUT AsnObjectIdentifier * enterprise,
+ OUT AsnInteger * genericTrap,
+ OUT AsnInteger * specificTrap,
+ OUT AsnTimeticks * timeStamp,
+ OUT RFC1157VarBindList * variableBindings);
+
+ typedef BOOL(WINAPI * pSnmpExtensionQuery) (
+ IN BYTE requestType,
+ IN OUT RFC1157VarBindList * variableBindings,
+ OUT AsnInteger * errorStatus,
+ OUT AsnInteger * errorIndex);
+
+ typedef BOOL(WINAPI * pSnmpExtensionInitEx) (
+ OUT AsnObjectIdentifier * supportedView);
#endif //_WIN32
-
+/// \brief gets current M.A.C adress (for internal use only)
int GetMacAddrSys ( unsigned char *addr )
{
#ifdef _WIN32
WSADATA WinsockData;
- if (WSAStartup(MAKEWORD(2, 0), &WinsockData) != 0)
+ if ( (WSAStartup(MAKEWORD(2, 0), &WinsockData)) != 0)
{
std::cerr << "This program requires Winsock 2.x!" << std::endl;
return -1;
HANDLE PollForTrapEvent;
AsnObjectIdentifier SupportedView;
- UINT OID_ifEntryType[] = { 1, 3, 6, 1, 2, 1, 2, 2, 1, 3 };
- UINT OID_ifEntryNum[] = { 1, 3, 6, 1, 2, 1, 2, 1 };
+ UINT OID_ifEntryType[] = { 1, 3, 6, 1, 2, 1, 2, 2, 1, 3 };
+ UINT OID_ifEntryNum[] = { 1, 3, 6, 1, 2, 1, 2, 1 };
UINT OID_ipMACEntAddr[] = { 1, 3, 6, 1, 2, 1, 2, 2, 1, 6 };
AsnObjectIdentifier MIB_ifMACEntAddr = {
sizeof(OID_ipMACEntAddr) / sizeof(UINT), OID_ipMACEntAddr };
// implementation for POSIX system
-#ifdef __sun
+#if defined(CMAKE_HAVE_NET_IF_ARP_H) && defined(__sun)
//The POSIX version is broken anyway on Solaris, plus would require full
//root power
struct arpreq parpreq;
// We should investiage the use of SIZEOF_ADDR_IFREQ
//
#ifdef HAVE_SA_LEN
-#ifndef max
-#define max(a,b) ((a) > (b) ? (a) : (b))
-#endif
-#define ifreq_size(i) max(sizeof(struct ifreq),\
- sizeof((i).ifr_name)+(i).ifr_addr.sa_len)
+ #ifndef max
+ #define max(a,b) ((a) > (b) ? (a) : (b))
+ #endif
+ #define ifreq_size(i) max(sizeof(struct ifreq),\
+ sizeof((i).ifr_name)+(i).ifr_addr.sa_len)
#else
-#define ifreq_size(i) sizeof(struct ifreq)
+ #define ifreq_size(i) sizeof(struct ifreq)
#endif // HAVE_SA_LEN
if( (sd = socket(AF_INET, SOCK_DGRAM, IPPROTO_IP)) < 0 )
#endif //__sun
}
+/**
+ * \brief Mini function to return the last digit from a number express in base 256
+ * pre condition data contain an array of 6 unsigned char
+ * post condition carry contain the last digit
+ */
+inline int getlastdigit(unsigned char *data)
+{
+ int extended, carry = 0;
+ for(int i=0;i<6;i++)
+ {
+ extended = (carry << 8) + data[i];
+ data[i] = extended / 10;
+ carry = extended % 10;
+ }
+ return carry;
+}
+
/**
* \brief Encode the mac address on a fixed lenght string of 15 characters.
* we save space this way.
// http://groups-beta.google.com/group/comp.unix.solaris/msg/ad36929d783d63be
// http://bdn.borland.com/article/0,1410,26040,00.html
unsigned char addr[6];
- uint64_t n = 0;
-
+
int stat = GetMacAddrSys(addr);
if (stat == 0)
{
- // Horner evaluation
- for(int i=0; i<6; i++)
+ // We need to convert a 6 digit number from base 256 to base 10, using integer
+ // would requires a 48bits one. To avoid this we have to reimplement the div + modulo
+ // with string only
+ bool zero = false;
+ int res;
+ std::string sres;
+ while(!zero)
{
- n *= 256;
- n += addr[i];
+ res = getlastdigit(addr);
+ sres.insert(sres.begin(), '0' + res);
+ zero = (addr[0] == 0) && (addr[1] == 0) && (addr[2] == 0)
+ && (addr[3] == 0) && (addr[4] == 0) && (addr[5] == 0);
}
- // we fit on 15 bytes maximum < 256^6.
-#if defined(_MSC_VER) || defined(__BORLANDC__)
- return Format("%I64u", n);
-#else
- return Format("%llu", n);
-#endif
+ return sres;
}
else
{
- gdcmVerboseMacro("Problem in finding the MAC Address");
+ gdcmWarningMacro("Problem in finding the MAC Address");
return "";
}
}
-/**
- * \brief Return the IP adress of the machine writting the DICOM image
- */
-std::string Util::GetIPAddress()
-{
- // This is a rip from
- // http://www.codeguru.com/Cpp/I-N/internet/network/article.php/c3445/
-#ifndef HOST_NAME_MAX
- // SUSv2 guarantees that `Host names are limited to 255 bytes'.
- // POSIX 1003.1-2001 guarantees that `Host names (not including the
- // terminating NUL) are limited to HOST_NAME_MAX bytes'.
-# define HOST_NAME_MAX 255
- // In this case we should maybe check the string was not truncated.
- // But I don't known how to check that...
-#if defined(_MSC_VER) || defined(__BORLANDC__) || defined(__MINGW32__)
- // with WinSock DLL we need to initialize the WinSock before using gethostname
- WORD wVersionRequested = MAKEWORD(1,0);
- WSADATA WSAData;
- int err = WSAStartup(wVersionRequested,&WSAData);
- if (err != 0)
- {
- // Tell the user that we could not find a usable
- // WinSock DLL.
- WSACleanup();
- return "127.0.0.1";
- }
-#endif
-
-#endif //HOST_NAME_MAX
-
- std::string str;
- char szHostName[HOST_NAME_MAX+1];
- int r = gethostname(szHostName, HOST_NAME_MAX);
-
- if( r == 0 )
- {
- // Get host adresses
- struct hostent *pHost = gethostbyname(szHostName);
-
- for( int i = 0; pHost!= NULL && pHost->h_addr_list[i]!= NULL; i++ )
- {
- for( int j = 0; j<pHost->h_length; j++ )
- {
- if( j > 0 ) str += ".";
-
- str += Util::Format("%u",
- (unsigned int)((unsigned char*)pHost->h_addr_list[i])[j]);
- }
- // str now contains one local IP address
-
-#if defined(_MSC_VER) || defined(__BORLANDC__) || defined(__MINGW32__)
- WSACleanup();
-#endif
-
- }
- }
- // If an error occur r == -1
- // Most of the time it will return 127.0.0.1...
- return str;
-}
-
/**
* \brief Creates a new UID. As stipulate in the DICOM ref
* each time a DICOM image is create it should have
*/
std::string Util::CreateUniqueUID(const std::string &root)
{
- std::string prefix = root;
+ std::string prefix;
std::string append;
if( root.empty() )
{
- // No root was specified use "GDCM" then
- // echo "gdcm" | od -b
- // 0000000 147 144 143 155 012
- prefix = "147.144.143.155"; // special easter egg
+ // gdcm UID prefix, as supplied by http://www.medicalconnections.co.uk
+ prefix = RootUID;
+ }
+ else
+ {
+ prefix = root;
}
- // else
+
// A root was specified use it to forge our new UID:
append += ".";
append += Util::GetMACAddress();
return prefix + append;
}
-unsigned int Util::GetCurrentThreadID()
+void Util::SetRootUID(const std::string &root)
{
-// FIXME the implementation is far from complete
-#if defined(_MSC_VER) || defined(__BORLANDC__) || defined(__MINGW32__)
- return (unsigned int)GetCurrentThreadId();
-#endif
-#ifdef __linux__
- return 0;
- // Doesn't work on fedora, but is in the man page...
- //return (unsigned int)gettid();
-#endif
-#ifdef __sun
- return (unsigned int)thr_self();
-#else
- //default implementation
- return 0;
-#endif
-}
-
-unsigned int Util::GetCurrentProcessID()
-{
-#if defined(_MSC_VER) || defined(__BORLANDC__) || defined(__MINGW32__)
- // NOTE: There is also a _getpid()...
- return (unsigned int)GetCurrentProcessId();
-#else
- // get process identification, POSIX
- return (unsigned int)getpid();
-#endif
-
+ if( root.empty() )
+ RootUID = GDCM_UID;
+ else
+ RootUID = root;
}
-/**
- * \brief
- * @param os ostream to write to
- * @param val val
- */
-template <class T>
-std::ostream &binary_write(std::ostream &os, const T &val)
+const std::string &Util::GetRootUID()
{
- return os.write(reinterpret_cast<const char*>(&val), sizeof val);
+ return RootUID;
}
+//-------------------------------------------------------------------------
/**
* \brief binary_write binary_write
* @param os ostream to write to
*/
std::ostream &binary_write(std::ostream &os, const uint16_t &val)
{
-#ifdef GDCM_WORDS_BIGENDIAN
+#if defined(GDCM_WORDS_BIGENDIAN) || defined(GDCM_FORCE_BIGENDIAN_EMULATION)
uint16_t swap;
- swap = ((( val << 8 ) & 0x0ff00 ) | (( val >> 8 ) & 0x00ff ) );
+ //swap = ((( val << 8 ) & 0xff00 ) | (( val >> 8 ) & 0x00ff ) );
+ //save CPU time
+ swap = ( val << 8 | val >> 8 );
+
return os.write(reinterpret_cast<const char*>(&swap), 2);
#else
return os.write(reinterpret_cast<const char*>(&val), 2);
*/
std::ostream &binary_write(std::ostream &os, const uint32_t &val)
{
-#ifdef GDCM_WORDS_BIGENDIAN
+#if defined(GDCM_WORDS_BIGENDIAN) || defined(GDCM_FORCE_BIGENDIAN_EMULATION)
uint32_t swap;
- swap = ( ((val<<24) & 0xff000000) | ((val<<8) & 0x00ff0000) |
- ((val>>8) & 0x0000ff00) | ((val>>24) & 0x000000ff) );
+// swap = ( ((val<<24) & 0xff000000) | ((val<<8) & 0x00ff0000) |
+// ((val>>8) & 0x0000ff00) | ((val>>24) & 0x000000ff) );
+// save CPU time
+ swap = ( (val<<24) | ((val<<8) & 0x00ff0000) |
+ ((val>>8) & 0x0000ff00) | (val>>24) );
return os.write(reinterpret_cast<const char*>(&swap), 4);
#else
return os.write(reinterpret_cast<const char*>(&val), 4);
return os.write(val.c_str(), val.size());
}
+/**
+ * \brief binary_write binary_write
+ * @param os ostream to write to
+ * @param val value
+ * @param len length of the 'value' to be written
+ */
+std::ostream &binary_write(std::ostream &os, const uint8_t *val, size_t len)
+{
+ // We are writting sizeof(char) thus no need to swap bytes
+ return os.write(reinterpret_cast<const char*>(val), len);
+}
+
+/**
+ * \brief binary_write binary_write
+ * @param os ostream to write to
+ * @param val val
+ * @param len length of the 'value' to be written
+ */
+std::ostream &binary_write(std::ostream &os, const uint16_t *val, size_t len)
+{
+// This is tricky since we are writting two bytes buffer.
+// Be carefull with little endian vs big endian.
+// Also this other trick is to allocate a small (efficient) buffer that store
+// intermediate result before writting it.
+#if defined(GDCM_WORDS_BIGENDIAN) || defined(GDCM_FORCE_BIGENDIAN_EMULATION)
+ const int BUFFER_SIZE = 4096;
+ static char buffer[BUFFER_SIZE];
+ uint16_t *binArea16 = (uint16_t*)val; //for the const
+
+ // how many BUFFER_SIZE long pieces in binArea ?
+ int nbPieces = len/BUFFER_SIZE; //(16 bits = 2 Bytes)
+ int remainingSize = len%BUFFER_SIZE;
+
+ for (int j=0;j<nbPieces;j++)
+ {
+ uint16_t *pbuffer = (uint16_t*)buffer; //reinitialize pbuffer
+ for (int i = 0; i < BUFFER_SIZE/2; i++)
+ {
+ *pbuffer = *binArea16 >> 8 | *binArea16 << 8;
+ pbuffer++;
+ binArea16++;
+ }
+ os.write ( buffer, BUFFER_SIZE );
+ }
+ if ( remainingSize > 0)
+ {
+ uint16_t *pbuffer = (uint16_t*)buffer; //reinitialize pbuffer
+ for (int i = 0; i < remainingSize/2; i++)
+ {
+ *pbuffer = *binArea16 >> 8 | *binArea16 << 8;
+ pbuffer++;
+ binArea16++;
+ }
+ os.write ( buffer, remainingSize );
+ }
+ return os;
+#else
+ return os.write(reinterpret_cast<const char*>(val), len);
+#endif
+}
+
+//-------------------------------------------------------------------------
+// Protected
+
+//-------------------------------------------------------------------------
+// Private
+/**
+ * \brief Return the IP adress of the machine writting the DICOM image
+ */
+std::string Util::GetIPAddress()
+{
+ // This is a rip from
+ // http://www.codeguru.com/Cpp/I-N/internet/network/article.php/c3445/
+#ifndef HOST_NAME_MAX
+ // SUSv2 guarantees that `Host names are limited to 255 bytes'.
+ // POSIX 1003.1-2001 guarantees that `Host names (not including the
+ // terminating NUL) are limited to HOST_NAME_MAX bytes'.
+#define HOST_NAME_MAX 255
+ // In this case we should maybe check the string was not truncated.
+ // But I don't known how to check that...
+#if defined(_MSC_VER) || defined(__BORLANDC__) || defined(__MINGW32__)
+ // with WinSock DLL we need to initialize the WinSock before using gethostname
+ WORD wVersionRequested = MAKEWORD(1,0);
+ WSADATA WSAData;
+ int err = WSAStartup(wVersionRequested,&WSAData);
+ if (err != 0)
+ {
+ // Tell the user that we could not find a usable
+ // WinSock DLL.
+ WSACleanup();
+ return "127.0.0.1";
+ }
+#endif
+
+#endif //HOST_NAME_MAX
+
+ std::string str;
+ char szHostName[HOST_NAME_MAX+1];
+ int r = gethostname(szHostName, HOST_NAME_MAX);
+
+ if( r == 0 )
+ {
+ // Get host adresses
+ struct hostent *pHost = gethostbyname(szHostName);
+
+ for( int i = 0; pHost!= NULL && pHost->h_addr_list[i]!= NULL; i++ )
+ {
+ for( int j = 0; j<pHost->h_length; j++ )
+ {
+ if( j > 0 ) str += ".";
+
+ str += Util::Format("%u",
+ (unsigned int)((unsigned char*)pHost->h_addr_list[i])[j]);
+ }
+ // str now contains one local IP address
+
+#if defined(_MSC_VER) || defined(__BORLANDC__) || defined(__MINGW32__)
+ WSACleanup();
+#endif
+
+ }
+ }
+ // If an error occur r == -1
+ // Most of the time it will return 127.0.0.1...
+ return str;
+}
+
+//-------------------------------------------------------------------------
} // end namespace gdcm