「mrcImage(API)」の版間の差分
提供: Eospedia
行2: | 行2: | ||
==定数== | ==定数== | ||
− | === | + | ===マスク処理=== |
typedef unsigned int mrcStatusMask; | typedef unsigned int mrcStatusMask; | ||
typedef unsigned int mrcStatusType; | typedef unsigned int mrcStatusType; | ||
行22: | 行22: | ||
#define mrcStatusMask_ReferMemory ((mrcStatusMask)0x00020000L) | #define mrcStatusMask_ReferMemory ((mrcStatusMask)0x00020000L) | ||
#define mrcStatusMask_ReferFile ((mrcStatusMask)0x00010000L) | #define mrcStatusMask_ReferFile ((mrcStatusMask)0x00010000L) | ||
− | + | ||
+ | ===mrcファイルの状態=== | ||
#define MRC_Status_OK ((mrcStatusType)0x00000000L) | #define MRC_Status_OK ((mrcStatusType)0x00000000L) | ||
行41: | 行42: | ||
extern char* mrcImageGeneralTypeList[]; | extern char* mrcImageGeneralTypeList[]; | ||
− | === | + | ===ヘッダ情報=== |
====Header Size==== | ====Header Size==== | ||
#define MRC_HEADER (1024) | #define MRC_HEADER (1024) | ||
行47: | 行48: | ||
====Image Data Type==== | ====Image Data Type==== | ||
− | + | =====for compatibility===== | |
− | typedef mrcImageParaTypeInteger mrcImageMode; | + | typedef mrcImageParaTypeInteger mrcImageMode; |
− | #define MRC_char_image ((mrcImageMode)0) | + | #define MRC_char_image ((mrcImageMode)0) |
#define MRC_short_image ((mrcImageMode)1) | #define MRC_short_image ((mrcImageMode)1) | ||
#define MRC_float_image ((mrcImageMode)2) | #define MRC_float_image ((mrcImageMode)2) | ||
行62: | 行63: | ||
#define MRC_complex_long_fft ((mrcImageMode)1103) | #define MRC_complex_long_fft ((mrcImageMode)1103) | ||
#define MRC_complex_double_fft ((mrcImageMode)1104) | #define MRC_complex_double_fft ((mrcImageMode)1104) | ||
− | #define MRC_NotSuppotedMode ((mrcImageMode)0xffffffffL) | + | #define MRC_NotSuppotedMode ((mrcImageMode)0xffffffffL) |
− | + | ||
− | #define mrcCharImage ((mrcImageMode)0) | + | #define mrcCharImage ((mrcImageMode)0) |
#define mrcShortImage ((mrcImageMode)1) | #define mrcShortImage ((mrcImageMode)1) | ||
#define mrcFloatImage ((mrcImageMode)2) | #define mrcFloatImage ((mrcImageMode)2) | ||
行70: | 行71: | ||
#define mrcComplexFloatFT ((mrcImageMode)4) | #define mrcComplexFloatFT ((mrcImageMode)4) | ||
#define mrcComplexFloatOneDimFT ((mrcImageMode)5) | #define mrcComplexFloatOneDimFT ((mrcImageMode)5) | ||
− | + | ||
+ | =====Change of IMOD===== | ||
#define mrcUShortImage ((mrcImageMode)6) | #define mrcUShortImage ((mrcImageMode)6) | ||
#define mrcCharImageRGB ((mrcImageMode)16) | #define mrcCharImageRGB ((mrcImageMode)16) | ||
− | + | ||
+ | =====Change for IMOD: (6,7) -> (26,27)===== | ||
#define mrcFloatRadonImage ((mrcImageMode)26) | #define mrcFloatRadonImage ((mrcImageMode)26) | ||
#define mrcFloatRadonFT ((mrcImageMode)27) | #define mrcFloatRadonFT ((mrcImageMode)27) | ||
行84: | 行87: | ||
#define mrcComplexLongFFT ((mrcImageMode)1103) | #define mrcComplexLongFFT ((mrcImageMode)1103) | ||
#define mrcComplexDoubleFFT ((mrcImageMode)1104) | #define mrcComplexDoubleFFT ((mrcImageMode)1104) | ||
− | #define mrcNotSuppotedMode ((mrcImageMode)0xffffffffL) | + | #define mrcNotSuppotedMode ((mrcImageMode)0xffffffffL) |
− | + | ||
#define MRC_MAX_LABEL_N (10) | #define MRC_MAX_LABEL_N (10) | ||
#define MRC_MAX_LABEL_LEN (80) | #define MRC_MAX_LABEL_LEN (80) | ||
行91: | 行94: | ||
#define MRC_MAX_EXTRA_CCP4 (15) | #define MRC_MAX_EXTRA_CCP4 (15) | ||
− | === | + | ===トレイラ情報=== |
#define mrcImageTailerMode2DProjection 0 | #define mrcImageTailerMode2DProjection 0 | ||
行97: | 行100: | ||
typedef enum mrcImageInformationMode { | typedef enum mrcImageInformationMode { | ||
meanOfAll = 0, | meanOfAll = 0, | ||
− | + | meanOfSphereSurface = 1, | |
− | + | meanOfCylinderSurface = 2, | |
− | + | meanOfEdge = 3, | |
− | + | meanOfCentre = 4, | |
− | + | meanOfSparse = 5, | |
meanOf2DEdgeY = 6, | meanOf2DEdgeY = 6, | ||
RMSofAllPixels = 7 | RMSofAllPixels = 7 | ||
行115: | 行118: | ||
mrcPixelPhase | mrcPixelPhase | ||
} mrcPixelDataType; | } mrcPixelDataType; | ||
− | + | ||
typedef enum mrcPixelDataHowToGet { | typedef enum mrcPixelDataHowToGet { | ||
mrcPixelHowNearest = 0, | mrcPixelHowNearest = 0, | ||
行124: | 行127: | ||
==構造体== | ==構造体== | ||
− | === | + | ===マスク処理=== |
typedef union MRC_Status_t { | typedef union MRC_Status_t { | ||
mrcStatusType all; | mrcStatusType all; | ||
行149: | 行152: | ||
typedef unsigned char mrcImageParaTypeCharacter; | typedef unsigned char mrcImageParaTypeCharacter; | ||
+ | 座標情報(int) | ||
typedef struct mrcImageParaTypeIntegerCoord { | typedef struct mrcImageParaTypeIntegerCoord { | ||
mrcImageParaTypeInteger x; | mrcImageParaTypeInteger x; | ||
行155: | 行159: | ||
} mrcImageParaTypeIntegerCoord; | } mrcImageParaTypeIntegerCoord; | ||
+ | 座標情報(実数) | ||
typedef struct mrcImageParaTypeRealCoord { | typedef struct mrcImageParaTypeRealCoord { | ||
mrcImageParaTypeReal x; | mrcImageParaTypeReal x; | ||
行160: | 行165: | ||
mrcImageParaTypeReal z; | mrcImageParaTypeReal z; | ||
} mrcImageParaTypeRealCoord; | } mrcImageParaTypeRealCoord; | ||
− | + | ||
− | === | + | ===ヘッダ情報=== |
+ | ====デフォルト==== | ||
typedef struct _mrcImageHeader { | typedef struct _mrcImageHeader { | ||
mrcImageParaTypeIntegerCoord N; /* Data Numbers */ | mrcImageParaTypeIntegerCoord N; /* Data Numbers */ | ||
行180: | 行186: | ||
mrcImageParaTypeInteger NSYMBT; /* Number of bytes */ | mrcImageParaTypeInteger NSYMBT; /* Number of bytes */ | ||
mrcImageParaTypeReal EXTRA[MRC_MAX_EXTRA]; /* Extra, user defined storage space */ | mrcImageParaTypeReal EXTRA[MRC_MAX_EXTRA]; /* Extra, user defined storage space */ | ||
− | + | mrcImageParaTypeReal OriginX; | |
− | + | mrcImageParaTypeReal OriginY; | |
− | + | mrcImageParaTypeInteger LabelN; | |
− | + | char Label[MRC_MAX_LABEL_N][MRC_MAX_LABEL_LEN]; | |
} _mrcImageHeader; | } _mrcImageHeader; | ||
− | + | ====CCP4 2000==== | |
typedef struct _mrcImageHeaderCCP4 { | typedef struct _mrcImageHeaderCCP4 { | ||
mrcImageParaTypeIntegerCoord N; /* Data Numbers */ | mrcImageParaTypeIntegerCoord N; /* Data Numbers */ | ||
行208: | 行214: | ||
mrcImageParaTypeReal SKWTRN[3]; /* Skew translation */ | mrcImageParaTypeReal SKWTRN[3]; /* Skew translation */ | ||
mrcImageParaTypeReal EXTRA[MRC_MAX_EXTRA_CCP4]; /* Extra, user defined storage space */ | mrcImageParaTypeReal EXTRA[MRC_MAX_EXTRA_CCP4]; /* Extra, user defined storage space */ | ||
− | + | mrcImageParaTypeCharacter MAP[4]; | |
− | + | mrcImageParaTypeInteger MARCHST; | |
− | + | mrcImageParaTypeReal ARMS; | |
− | + | mrcImageParaTypeInteger LabelN; | |
− | + | char Label[MRC_MAX_LABEL_N][MRC_MAX_LABEL_LEN]; | |
} _mrcImageHeaderCCP4; | } _mrcImageHeaderCCP4; | ||
− | + | ====IMOD==== | |
typedef struct _mrcImageHeaderIMOD { | typedef struct _mrcImageHeaderIMOD { | ||
mrcImageParaTypeIntegerCoord N; /* Data Numbers */ | mrcImageParaTypeIntegerCoord N; /* Data Numbers */ | ||
行250: | 行256: | ||
mrcImageParaTypeReal yorg; /* Extra, user defined storage space */ | mrcImageParaTypeReal yorg; /* Extra, user defined storage space */ | ||
mrcImageParaTypeReal zorg; /* Extra, user defined storage space */ | mrcImageParaTypeReal zorg; /* Extra, user defined storage space */ | ||
− | + | mrcImageParaTypeCharacter MAP[4]; | |
− | + | mrcImageParaTypeInteger MARCHST; | |
− | + | mrcImageParaTypeReal ARMS; | |
− | + | mrcImageParaTypeInteger LabelN; | |
− | + | char Label[MRC_MAX_LABEL_N][MRC_MAX_LABEL_LEN]; | |
} _mrcImageHeaderIMOD; | } _mrcImageHeaderIMOD; | ||
− | + | ||
− | + | ====共用体==== | |
typedef union mrcImageHeader { | typedef union mrcImageHeader { | ||
unsigned char All[MRC_HEADER]; | unsigned char All[MRC_HEADER]; | ||
行267: | 行273: | ||
} mrcImageHeader; | } mrcImageHeader; | ||
− | + | ====FEI==== | |
#define FEI_EXTENDED_HEADER_EACH (128) | #define FEI_EXTENDED_HEADER_EACH (128) | ||
行299: | 行305: | ||
} FEIextendedHeader; | } FEIextendedHeader; | ||
− | === | + | ===トレイラ情報=== |
typedef struct _mrcImageTailer { | typedef struct _mrcImageTailer { | ||
char Code[4]; | char Code[4]; | ||
行309: | 行315: | ||
mrcImageParaTypeInteger Mirror; | mrcImageParaTypeInteger Mirror; | ||
} _mrcImageTailer; | } _mrcImageTailer; | ||
− | + | ||
− | + | ||
typedef union mrcImageTailer { | typedef union mrcImageTailer { | ||
unsigned char All[MRC_TAILER]; | unsigned char All[MRC_TAILER]; | ||
行317: | 行322: | ||
_mrcImageTailer Cont; | _mrcImageTailer Cont; | ||
} mrcImageTailer; | } mrcImageTailer; | ||
− | + | ||
− | + | ||
typedef char mrcImageSymmetryOperator; | typedef char mrcImageSymmetryOperator; | ||
− | ===mrcImage | + | ===[[mrcImage]]の構造=== |
typedef struct mrcImage { | typedef struct mrcImage { | ||
− | + | mrcImageHeader Header; | |
− | + | size_t BytePerImage; | |
− | + | mrcImageParaTypeInteger BytePerBand; | |
− | + | mrcImageParaTypeInteger BandPerPixel; | |
− | + | size_t PixelPerImage; | |
− | + | mrcStatusType status; | |
− | + | char* Image; | |
− | + | unsigned char* ImageCharImage; | |
− | + | unsigned short* ImageShortImage; | |
− | + | unsigned long* ImageLongImage; | |
− | + | float* ImageFloatImage; | |
− | + | double* ImageDoubleImage; | |
/* Tailer */ | /* Tailer */ | ||
行356: | 行360: | ||
mrcImageSymmetryOperator* SymmetryOperator; | mrcImageSymmetryOperator* SymmetryOperator; | ||
} mrcImage; | } mrcImage; | ||
− | + | ||
+ | ===[[mrcImage]]の情報=== | ||
typedef struct mrcImageInformation { | typedef struct mrcImageInformation { | ||
− | + | double mean; | |
− | + | double rootMeanSquare; | |
− | + | double meanOfSphereSurface; | |
− | + | double sdOfSphereSurface; | |
− | + | double seOfSphereSurface; | |
− | + | double meanOfCylinderSurface; | |
− | + | double sdOfCylinderSurface; | |
− | + | double seOfCylinderSurface; | |
− | + | double meanOfEdge; | |
− | + | double sdOfEdge; | |
− | + | double seOfEdge; | |
− | + | double max; | |
mrcImageParaTypeRealCoord maxCoord; | mrcImageParaTypeRealCoord maxCoord; | ||
− | + | double min; | |
mrcImageParaTypeRealCoord minCoord; | mrcImageParaTypeRealCoord minCoord; | ||
− | + | double sd; | |
− | + | double se; | |
double rmax; | double rmax; | ||
long flagrmax; | long flagrmax; | ||
行397: | 行402: | ||
} mrcImageInformation; | } mrcImageInformation; | ||
+ | ===メンバー呼び出しの省略表記を定義=== | ||
#define HeaderN Header.Cont.N | #define HeaderN Header.Cont.N | ||
#define HeaderMode Header.Cont.Mode | #define HeaderMode Header.Cont.Mode | ||
行607: | 行613: | ||
mrcPixelDataType mode); | mrcPixelDataType mode); | ||
− | ==== | + | ====インライン型関数の画素データの設定==== |
extern double mrcPixelDataSet(mrcImage* mrc, | extern double mrcPixelDataSet(mrcImage* mrc, | ||
mrcImageParaTypeReal x /* Pixel Unit */, | mrcImageParaTypeReal x /* Pixel Unit */, | ||
行681: | 行687: | ||
extern int mrcImageByteSwap(mrcImage* mrc, long mode); | extern int mrcImageByteSwap(mrcImage* mrc, long mode); | ||
− | === | + | ===投影=== |
extern void lmrcImageXSection(mrcImage* proj, mrcImage* img, double Y); | extern void lmrcImageXSection(mrcImage* proj, mrcImage* img, double Y); | ||
extern void lmrcImageYSection(mrcImage* proj, mrcImage* img, double X); | extern void lmrcImageYSection(mrcImage* proj, mrcImage* img, double X); | ||
行700: | 行706: | ||
extern void lmrcImageInformation(mrcImageInformation* info, mrcImage* img); | extern void lmrcImageInformation(mrcImageInformation* info, mrcImage* img); | ||
− | ===mrcImageCopy=== | + | ===データコピー(mrcImageCopy)=== |
extern void lmrcImageCopy(mrcImage* dst, mrcImage* src, mrcImageParaTypeRealCoord to); | extern void lmrcImageCopy(mrcImage* dst, mrcImage* src, mrcImageParaTypeRealCoord to); | ||
− | === | + | ===画像モードのチェック=== |
extern int IsFloatImage(mrcImage* mrc, char* message, mrcImageParaTypeInteger mode); | extern int IsFloatImage(mrcImage* mrc, char* message, mrcImageParaTypeInteger mode); | ||
extern int IsFloatFT(mrcImage* mrc, char* message, mrcImageParaTypeInteger mode); | extern int IsFloatFT(mrcImage* mrc, char* message, mrcImageParaTypeInteger mode); | ||
行709: | 行715: | ||
extern int IsFT(mrcImage* mrc, char* message, mrcImageParaTypeInteger mode); | extern int IsFT(mrcImage* mrc, char* message, mrcImageParaTypeInteger mode); | ||
extern int IsFFT(mrcImage* mrc, char* message, mrcImageParaTypeInteger mode); | extern int IsFFT(mrcImage* mrc, char* message, mrcImageParaTypeInteger mode); | ||
+ | |||
+ | ===画像情報の設定=== | ||
extern void mrcImageInfoSet(mrcImage* mrc, FILE* fptIn, FILE* fptOut, long mode); | extern void mrcImageInfoSet(mrcImage* mrc, FILE* fptIn, FILE* fptOut, long mode); | ||
===mrcImageToIntImage=== | ===mrcImageToIntImage=== | ||
extern void lmrcImageToIntImage(mrcImage* dst, mrcImage* src, mrcImageParaTypeInteger mode); | extern void lmrcImageToIntImage(mrcImage* dst, mrcImage* src, mrcImageParaTypeInteger mode); | ||
− | + | ||
+ | ===FFT関連=== | ||
extern void lfft2d(mrcImage* fft, mrcImage* img); | extern void lfft2d(mrcImage* fft, mrcImage* img); | ||
2014年5月13日 (火) 06:11時点における版
DataManip/mrcImageは、mrcImage型の画像を取り扱うための低水準APIです。
目次
- 1 定数
- 2 構造体
- 3 API
- 3.1 初期化
- 3.2 ファイルの読み込み
- 3.3 ファイルの書き出し
- 3.4 画素のデータの取得
- 3.5 画素のデータの設定
- 3.6 画像処理一般のライブラリ
- 3.7 エンディアン関連
- 3.8 投影
- 3.9 データコピー(mrcImageCopy)
- 3.10 画像モードのチェック
- 3.11 画像情報の設定
- 3.12 mrcImageToIntImage
- 3.13 FFT関連
- 3.14 mrcError
- 3.15 lmrcFFTInfo
- 3.16 mrcImageCheckSameSize
- 3.17 mrcImageCheckFFT
- 3.18 lmrcImageSmoothing
- 3.19 lmrcImageCVE
- 3.20 lmrcImageTfunction
- 3.21 lmrcImageHighlighting
- 3.22 lmrcImageLowPassFilter
- 3.23 lmrcImageHighPassFilter
- 3.24 lmrcImageBandPassFilter
- 3.25 lmrcImageShrink
- 3.26 lmrcImageShift
- 3.27 mrcImagePrint
- 3.28 lmrcImageMultiCTFCompensation
- 3.29 lmrcImageCTFObservation
- 3.30 lmrcImageCTFCompensation
- 3.31 lmrcImageCTFSN
- 3.32 lmrcImageFourierPowerSpectrum
- 3.33 mrcImageSectionGet
- 3.34 mrcImageSectionSet
定数
マスク処理
typedef unsigned int mrcStatusMask; typedef unsigned int mrcStatusType;
#define mrcStatusMask_ExitFailure ((mrcStatusMask)0x80000000L) #define mrcStatusMask_Warning ((mrcStatusMask)0x40000000L) #define mrcStatusMask_Information ((mrcStatusMask)0x20000000L) #define mrcStatusMask_Message ((mrcStatusMask)0x10000000L) #define mrcStatusMask_1B ((mrcStatusMask)0x08000000L) #define mrcStatusMask_1A ((mrcStatusMask)0x04000000L) #define mrcStatusMask_19 ((mrcStatusMask)0x02000000L) #define mrcStatusMask_18 ((mrcStatusMask)0x01000000L) #define mrcStatusMask_17 ((mrcStatusMask)0x00800000L) #define mrcStatusMask_16 ((mrcStatusMask)0x00400000L) #define mrcStatusMask_15 ((mrcStatusMask)0x00200000L) #define mrcStatusMask_14 ((mrcStatusMask)0x00100000L) #define mrcStatusMask_ReferMRCInfo ((mrcStatusMask)0x00080000L) #define mrcStatusMask_Mode ((mrcStatusMask)0x00040000L) #define mrcStatusMask_ReferMemory ((mrcStatusMask)0x00020000L) #define mrcStatusMask_ReferFile ((mrcStatusMask)0x00010000L)
mrcファイルの状態
#define MRC_Status_OK ((mrcStatusType)0x00000000L) #define MRC_FileCanNotOpened ((mrcStatusType)0x00000000L|mrcStatusMask_ReferFile) #define MRC_ImageHeaderCanNotRead ((mrcStatusType)0x00000001L|mrcStatusMask_ReferFile) #define MRC_ImageDataCanNotRead ((mrcStatusType)0x00000002L|mrcStatusMask_ReferFile) #define MRC_ImageHeaderCanNotWrite ((mrcStatusType)0x00000003L|mrcStatusMask_ReferFile) #define MRC_ImageDataCanNotWrite ((mrcStatusType)0x00000004L|mrcStatusMask_ReferFile) #define MRC_FileCanNotClosed ((mrcStatusType)0x00000005L|mrcStatusMask_ReferFile) #define MRC_MemoryNotAllocate ((mrcStatusType)0x00000000L|mrcStatusMask_ReferMemory) #define MRC_ModeNotSupported ((mrcStatusType)0x00000000L|mrcStatusMask_Mode) #define MRC_ModeZNot1FFTMode ((mrcStatusType)0x00000001L|mrcStatusMask_Mode) #define MRC_SizeOfZero ((mrcStatusType)0x00000000L|mrcStatusMask_ReferMRCInfo)
extern char* mrcImageGeneralTypeList[];
ヘッダ情報
Header Size
#define MRC_HEADER (1024) #define MRC_TAILER (1024)
Image Data Type
for compatibility
typedef mrcImageParaTypeInteger mrcImageMode; #define MRC_char_image ((mrcImageMode)0) #define MRC_short_image ((mrcImageMode)1) #define MRC_float_image ((mrcImageMode)2) #define MRC_complex_short_ft ((mrcImageMode)3) #define MRC_complex_float_ft ((mrcImageMode)4) #define MRC_long_image ((mrcImageMode)101) #define MRC_double_image ((mrcImageMode)102) #define MRC_complex_long_ft ((mrcImageMode)103) #define MRC_complex_double_ft ((mrcImageMode)104) #define MRC_complex_short_fft ((mrcImageMode)1003) #define MRC_complex_float_fft ((mrcImageMode)1004) #define MRC_complex_long_fft ((mrcImageMode)1103) #define MRC_complex_double_fft ((mrcImageMode)1104) #define MRC_NotSuppotedMode ((mrcImageMode)0xffffffffL)
#define mrcCharImage ((mrcImageMode)0) #define mrcShortImage ((mrcImageMode)1) #define mrcFloatImage ((mrcImageMode)2) #define mrcComplexShortFT ((mrcImageMode)3) #define mrcComplexFloatFT ((mrcImageMode)4) #define mrcComplexFloatOneDimFT ((mrcImageMode)5)
Change of IMOD
#define mrcUShortImage ((mrcImageMode)6) #define mrcCharImageRGB ((mrcImageMode)16)
Change for IMOD: (6,7) -> (26,27)
#define mrcFloatRadonImage ((mrcImageMode)26) #define mrcFloatRadonFT ((mrcImageMode)27) #define mrcLongImage ((mrcImageMode)101) #define mrcDoubleImage ((mrcImageMode)102) #define mrcComplexLongFT ((mrcImageMode)103) #define mrcComplexDoubleFT ((mrcImageMode)104) #define mrcComplexShortFFT ((mrcImageMode)1003) #define mrcComplexFloatFFT ((mrcImageMode)1004) #define mrcComplexLongFFT ((mrcImageMode)1103) #define mrcComplexDoubleFFT ((mrcImageMode)1104) #define mrcNotSuppotedMode ((mrcImageMode)0xffffffffL)
#define MRC_MAX_LABEL_N (10) #define MRC_MAX_LABEL_LEN (80) #define MRC_MAX_EXTRA (29) #define MRC_MAX_EXTRA_CCP4 (15)
トレイラ情報
#define mrcImageTailerMode2DProjection 0
画像情報
typedef enum mrcImageInformationMode { meanOfAll = 0, meanOfSphereSurface = 1, meanOfCylinderSurface = 2, meanOfEdge = 3, meanOfCentre = 4, meanOfSparse = 5, meanOf2DEdgeY = 6, RMSofAllPixels = 7 } mrcImageInformationMode;
typedef enum mrcPixelDataType { mrcPixelRePart = 0, mrcPixelImPart, mrcPixelMag, mrcPixelPow, mrcPixelLogMag, mrcPixelLogPow, mrcPixelPhase } mrcPixelDataType;
typedef enum mrcPixelDataHowToGet { mrcPixelHowNearest = 0, mrcPixelHowLinear, mrcPixelHowCubicConv, mrcPixelHowPolySig } mrcPixelDataHowToGet;
構造体
マスク処理
typedef union MRC_Status_t { mrcStatusType all; struct Content { #ifdef M68 unsigned char level; unsigned char category; short detail; #else short detail; unsigned char category; unsigned char level; #endif } status; } MRC_Status_t;
#ifdef M68 typedef int mrcImageParaTypeInteger; #else typedef int mrcImageParaTypeInteger; #endif
typedef float mrcImageParaTypeReal; typedef unsigned char mrcImageParaTypeCharacter;
座標情報(int)
typedef struct mrcImageParaTypeIntegerCoord { mrcImageParaTypeInteger x; mrcImageParaTypeInteger y; mrcImageParaTypeInteger z; } mrcImageParaTypeIntegerCoord;
座標情報(実数)
typedef struct mrcImageParaTypeRealCoord { mrcImageParaTypeReal x; mrcImageParaTypeReal y; mrcImageParaTypeReal z; } mrcImageParaTypeRealCoord;
ヘッダ情報
デフォルト
typedef struct _mrcImageHeader { mrcImageParaTypeIntegerCoord N; /* Data Numbers */ mrcImageMode Mode; /* Data Type */ mrcImageParaTypeIntegerCoord StartN; /* First in map */ mrcImageParaTypeIntegerCoord M; /* Number of Intervals */ mrcImageParaTypeRealCoord Length; /* Cell Dimensions : Angstrom */ mrcImageParaTypeReal Alpha; /* Cell Angle : Degree */ mrcImageParaTypeReal Beta; mrcImageParaTypeReal Gamma; mrcImageParaTypeInteger MAPC; /* Which Axis Corresponds to Columns */ mrcImageParaTypeInteger MAPR; /* Which Axis Corresponds to Rows */ mrcImageParaTypeInteger MAPS; /* Which Axis Corresponds to Sections */ mrcImageParaTypeReal AMin; /* Minimum Density Value */ mrcImageParaTypeReal AMax; /* Maximum Density Value */ mrcImageParaTypeReal AMean; /* Mean Density Value */ mrcImageParaTypeInteger ISPG; /* Space Group */ mrcImageParaTypeInteger NSYMBT; /* Number of bytes */ mrcImageParaTypeReal EXTRA[MRC_MAX_EXTRA]; /* Extra, user defined storage space */ mrcImageParaTypeReal OriginX; mrcImageParaTypeReal OriginY; mrcImageParaTypeInteger LabelN; char Label[MRC_MAX_LABEL_N][MRC_MAX_LABEL_LEN]; } _mrcImageHeader;
CCP4 2000
typedef struct _mrcImageHeaderCCP4 { mrcImageParaTypeIntegerCoord N; /* Data Numbers */ mrcImageMode Mode; /* Data Type */ mrcImageParaTypeIntegerCoord StartN; /* First in map */ mrcImageParaTypeIntegerCoord NI; /* Number of Intervals */ mrcImageParaTypeRealCoord Length; /* Cell Dimensions : Angstrom */ mrcImageParaTypeReal Alpha; /* Cell Angle : Degree */ mrcImageParaTypeReal Beta; mrcImageParaTypeReal Gamma; mrcImageParaTypeInteger MAPC; /* Which Axis Corresponds to Columns */ mrcImageParaTypeInteger MAPR; /* Which Axis Corresponds to Rows */ mrcImageParaTypeInteger MAPS; /* Which Axis Corresponds to Sections */ mrcImageParaTypeReal AMin; /* Minimum Density Value */ mrcImageParaTypeReal AMax; /* Maximum Density Value */ mrcImageParaTypeReal AMean; /* Mean Density Value */ mrcImageParaTypeInteger ISPG; /* Space Group */ mrcImageParaTypeInteger NSYMBT; /* Number of bytes used for storing symmetry operators */ mrcImageParaTypeInteger LSKFLG; /* Flag for Skew transformation */ mrcImageParaTypeReal SKWMAT[9]; /* Skew matrix S */ mrcImageParaTypeReal SKWTRN[3]; /* Skew translation */ mrcImageParaTypeReal EXTRA[MRC_MAX_EXTRA_CCP4]; /* Extra, user defined storage space */ mrcImageParaTypeCharacter MAP[4]; mrcImageParaTypeInteger MARCHST; mrcImageParaTypeReal ARMS; mrcImageParaTypeInteger LabelN; char Label[MRC_MAX_LABEL_N][MRC_MAX_LABEL_LEN]; } _mrcImageHeaderCCP4;
IMOD
typedef struct _mrcImageHeaderIMOD { mrcImageParaTypeIntegerCoord N; /* Data Numbers */ mrcImageMode Mode; /* Data Type */ mrcImageParaTypeIntegerCoord StartN; /* First in map */ mrcImageParaTypeIntegerCoord M; /* Number of Intervals */ mrcImageParaTypeRealCoord Length; /* PixelSize=Len/M: Angstrom */ mrcImageParaTypeReal Alpha; /* Cell Angle : Degree */ mrcImageParaTypeReal Beta; mrcImageParaTypeReal Gamma; mrcImageParaTypeInteger MAPC; /* Which Axis Corresponds to Columns */ mrcImageParaTypeInteger MAPR; /* Which Axis Corresponds to Rows */ mrcImageParaTypeInteger MAPS; /* Which Axis Corresponds to Sections */ mrcImageParaTypeReal AMin; /* Minimum Density Value */ mrcImageParaTypeReal AMax; /* Maximum Density Value */ mrcImageParaTypeReal AMean; /* Mean Density Value */ mrcImageParaTypeInteger ISPG; /* Space Group */ mrcImageParaTypeInteger NSYMBT; /* Number of bytes used for extended header */ short creatid; /* Flag for Skew transformation */ short extra1[15]; /* Flag for Skew transformation */ short nint; short nreal; /*l Skew matrix S */ mrcImageParaTypeReal extra2[5]; /* Extra, user defined storage space */ mrcImageParaTypeInteger imodStamp; mrcImageParaTypeInteger imodFlags; short idtype; short lnes; short nd1; short nd2; short vd1; short vd2; mrcImageParaTypeReal tiltangles[6]; /* Extra, user defined storage space */ mrcImageParaTypeReal xorg; /* Extra, user defined storage space */ mrcImageParaTypeReal yorg; /* Extra, user defined storage space */ mrcImageParaTypeReal zorg; /* Extra, user defined storage space */ mrcImageParaTypeCharacter MAP[4]; mrcImageParaTypeInteger MARCHST; mrcImageParaTypeReal ARMS; mrcImageParaTypeInteger LabelN; char Label[MRC_MAX_LABEL_N][MRC_MAX_LABEL_LEN]; } _mrcImageHeaderIMOD;
共用体
typedef union mrcImageHeader { unsigned char All[MRC_HEADER]; mrcImageParaTypeInteger Int[MRC_HEADER/4]; mrcImageParaTypeReal Real[MRC_HEADER/4]; _mrcImageHeader Cont; _mrcImageHeaderCCP4 ContCCP4; _mrcImageHeaderIMOD ContIMOD; } mrcImageHeader; ====FEI==== #define FEI_EXTENDED_HEADER_EACH (128) typedef struct _FEIextendedHeader { mrcImageParaTypeReal aTilt; // Alpha tilt, in degrees mrcImageParaTypeReal bTilt; // Beta tilt, in degrees mrcImageParaTypeReal xStage; // Stage x position. // Normally in SI units (meters), // but some older files may be in micrometers. // Check by looking at values for x,y,z. // If one of these exceeds 1, it will be micrometers. mrcImageParaTypeReal yStage; // Stage y position. For testing of units see x_stage. mrcImageParaTypeReal zStage; // Stage z position. For testing of units see x_stage. mrcImageParaTypeReal xShift; // Image shift x. For testing of units see x_stage. mrcImageParaTypeReal yShift; // Image shift y. For testing of units see x_stage. mrcImageParaTypeReal defocus; mrcImageParaTypeReal expTime; mrcImageParaTypeReal mean; mrcImageParaTypeReal tiltAxis; mrcImageParaTypeReal pixelSize; mrcImageParaTypeReal magnification; mrcImageParaTypeReal ht; mrcImageParaTypeReal binning; mrcImageParaTypeReal appliedDefocus; mrcImageParaTypeReal reminder[FEI_EXTENDED_HEADER_EACH/4-16]; } _FEIextendedHeader; typedef union FEIextendedHeader { mrcImageParaTypeReal Real[FEI_EXTENDED_HEADER_EACH/4]; _FEIextendedHeader Cont; } FEIextendedHeader;
トレイラ情報
typedef struct _mrcImageTailer { char Code[4]; mrcImageParaTypeInteger Mode; char EulerAngleMode[4]; mrcImageParaTypeReal Rot1; mrcImageParaTypeReal Rot2; mrcImageParaTypeReal Rot3; mrcImageParaTypeInteger Mirror; } _mrcImageTailer;
typedef union mrcImageTailer { unsigned char All[MRC_TAILER]; mrcImageParaTypeInteger Int[MRC_TAILER/4]; mrcImageParaTypeReal Real[MRC_TAILER/4]; _mrcImageTailer Cont; } mrcImageTailer;
typedef char mrcImageSymmetryOperator;
mrcImageの構造
typedef struct mrcImage { mrcImageHeader Header; size_t BytePerImage; mrcImageParaTypeInteger BytePerBand; mrcImageParaTypeInteger BandPerPixel; size_t PixelPerImage; mrcStatusType status; char* Image; unsigned char* ImageCharImage; unsigned short* ImageShortImage; unsigned long* ImageLongImage; float* ImageFloatImage; double* ImageDoubleImage; /* Tailer */ mrcImageTailer* Tailer; mrcImageParaTypeInteger numTailer; /* FEIextendedHeader */ FEIextendedHeader* FEIHeader; mrcImageParaTypeInteger numFEIHeader; /* Dummy */ char* dummyHeader; mrcImageParaTypeInteger dummyHeaderByte; /* Swap Bytes */ mrcImageParaTypeInteger flagByteSwap; /* flag CCP4-2000 */ mrcImageParaTypeInteger flagCCP4; mrcImageParaTypeInteger ImageOffset; mrcImageSymmetryOperator* SymmetryOperator; } mrcImage;
mrcImageの情報
typedef struct mrcImageInformation { double mean; double rootMeanSquare; double meanOfSphereSurface; double sdOfSphereSurface; double seOfSphereSurface; double meanOfCylinderSurface; double sdOfCylinderSurface; double seOfCylinderSurface; double meanOfEdge; double sdOfEdge; double seOfEdge; double max; mrcImageParaTypeRealCoord maxCoord; double min; mrcImageParaTypeRealCoord minCoord; double sd; double se; double rmax; long flagrmax; double RMS; /* Root mean square */ int mode; /* For Search Area */ int flagXrange; int flagYrange; int flagZrange; double XrangeMin; double XrangeMax; double YrangeMin; double YrangeMax; double ZrangeMin; double ZrangeMax; int flagCentre; double CentreX; double CentreY; double CentreZ; } mrcImageInformation;
メンバー呼び出しの省略表記を定義
#define HeaderN Header.Cont.N #define HeaderMode Header.Cont.Mode #define HeaderStartN Header.Cont.StartN #define HeaderM Header.Cont.M #define HeaderLength Header.Cont.Length #define HeaderAlpha Header.Cont.Alpha #define HeaderBeta Header.Cont.Beta #define HeaderGamma Header.Cont.Gamma #define HeaderMAPC Header.Cont.MAPC #define HeaderMAPR Header.Cont.MAPR #define HeaderMAPS Header.Cont.MAPS #define HeaderAMin Header.Cont.AMin #define HeaderAMax Header.Cont.AMax #define HeaderAMean Header.Cont.AMean #define HeaderISPG Header.Cont.ISPG #define HeaderNSYMBT Header.Cont.NSYMBT #define HeaderLSKFLG Header.ContCCP4.LSKFLG #define HeaderEXTRA Header.Cont.EXTRA #define HeaderOriginX Header.Cont.OriginX #define HeaderOriginY Header.Cont.OriginY #define HeaderLabelN Header.Cont.LabelN #define HeaderLabel Header.Cont.Label
画像情報
typedef struct lmrcImageSmoothingInfo { long mode; mrcImageParaTypeReal sx; /* Filter kernel size */ mrcImageParaTypeReal sy; mrcImageParaTypeReal sz; double sigma; /* for Lee-Sigma filter */ } lmrcImageSmoothingInfo; typedef struct lmrcImageCVEInfo { long sx; /* CVE domain size. */ long sy; /* Apr.30,1996 */ long sz; } lmrcImageCVEInfo; typedef struct lmrcImageTfunctionInfo { long sx; /* Tfunction domain size. */ long sy; /* June 11,1996 */ long sz; } lmrcImageTfunctionInfo; typedef struct lmrcImageHighlightInfo { long mode; float Bias; /* May 1,1996 */ float Grad; } lmrcImageHighlightInfo; typedef struct lmrcImageLowPassFilterInfo { long mode; /* June 4,1996 */ float hvp, width; } lmrcImageLowPassFilterInfo; typedef struct lmrcImageHighPassFilterInfo { long mode; /* June 5,1996 */ float hvp, width; float pathWindowLow; float pathWindowHigh; } lmrcImageHighPassFilterInfo; typedef struct lmrcImageBandPassFilterInfo { long mode; /* June 5,1996 */ float hvl, hvh, wl, wh; int flagLowWeight; double lowWeight; } lmrcImageBandPassFilterInfo;
mrcImageDataSet
typedef enum mrcImageDataSetMode { mrcImageDataSetModeCentre = 0, mrcImageDataSetModeOrigin = 1 } mrcImageDataSetMode;
lmrcImageMultiCTFCompensation
typedef enum lmrcImageMultiCTFCompensationSolventFlatteningMode { lmrcImageMultiCTFCompensationSolventFlatteningModeDensityPeakIsSolvent=0, lmrcImageMultiCTFCompensationSolventFlatteningModeDefinedSTDIsSolvent=1 } lmrcImageMultiCTFCompensationSolventFlatteningMode; typedef struct lmrcImageMultiCTFCompensationInfo { long maxiteration; double okrms; long diversefactor; double z; long nhist; FILE* logfile; lmrcImageMultiCTFCompensationSolventFlatteningMode solventMode; float definedSTD; } lmrcImageMultiCTFCompensationInfo;
lmrcImageCTFSN
typedef struct lmrcImageCTFSNInfo { long mode, deg1, deg2; ctfInfo ctf; } lmrcImageCTFSNInfo;
lmrcImageFourierPowerSpectrum
typedef struct mrcImageFourierPowerSpectrum { long n; double d, *h; } mrcImageFourierPowerSpectrum; typedef struct lmrcImageFourierPowerSpectrumInfo { double dX, dY, dZ; } lmrcImageFourierPowerSpectrumInfo;
API
初期化
extern mrcStatusType mrcInit(mrcImage* mrc, char* filaname); extern mrcStatusType mrcDefaultHeaderValueSet(mrcImage* mrc); extern mrcStatusType mrcHiddenDataSet(mrcImage* mrc, long mode); extern void mrcImageFree(mrcImage* mrc, char* message); extern void mrcTailerInit(mrcImage* mrc, long mode); extern void mrcTailerCopy(mrcImage* dst, mrcImage* src, long mode);
ファイルの読み込み
extern mrcStatusType mrcFileRead (mrcImage* mrc, char* filename, char* message, long mode); extern mrcStatusType mrcImageRead (mrcImage* mrc, char* filename, char* message, long mode); extern mrcStatusType mrcHeaderRead(mrcImage* mrc, char* filename, char* message, long mode); extern mrcStatusType mrcTailerRead(mrcImage* mrc, char* filename, char* message, long mode);
extern mrcStatusType mrcFileReadGeneral (mrcImage* mrc, char* filename, char* type, char* message, long mode); extern mrcStatusType mrcImageReadGeneral (mrcImage* mrc, char* filename, char* type, char* message, long mode); extern mrcStatusType mrcHeaderReadGeneral(mrcImage* mrc, char* filename, char* type, char* message, long mode); extern mrcStatusType mrcTailerReadGeneral(mrcImage* mrc, char* filename, char* type, char* message, long mode); extern mrcStatusType mrcFileGeneralListPrint(FILE* fpt);
extern mrcStatusType mrcFileReadFEIextendedMRC (mrcImage* mrc, char* filename, char* message, long mode); extern mrcStatusType mrcImageReadFEIextendedMRC (mrcImage* mrc, char* filename, char* message, long mode); extern mrcStatusType mrcHeaderReadFEIextendedMRC(mrcImage* mrc, char* filename, char* message, long mode); extern mrcStatusType mrcTailerReadFEIextendedMRC(mrcImage* mrc, char* filename, char* message, long mode);
extern mrcStatusType mrcFileReadIMODextendedMRC (mrcImage* mrc, char* filename, char* message, long mode); extern mrcStatusType mrcImageReadIMODextendedMRC (mrcImage* mrc, char* filename, char* message, long mode); extern mrcStatusType mrcHeaderReadIMODextendedMRC(mrcImage* mrc, char* filename, char* message, long mode); extern mrcStatusType mrcTailerReadIMODextendedMRC(mrcImage* mrc, char* filename, char* message, long mode);
ファイルの書き出し
extern mrcStatusType mrcFileWrite (mrcImage* mrc, char* filename, char* message, long mode); extern mrcStatusType mrcImageWrite (mrcImage* mrc, char* filename, char* message, long mode); extern mrcStatusType mrcHeaderWrite(mrcImage* mrc, char* filename, char* message, long mode); extern mrcStatusType mrcTailerWrite(mrcImage* mrc, char* filename, char* message, long mode); #define mrcFileWriteModeGet(x) BYTE4GETBYTE(x,0) #define mrcImageWriteModeGet(x) BYTE4GETBYTE(x,1) #define mrcHeaderWriteModeGet(x) BYTE4GETBYTE(x,2) #define mrcTailerWriteModeGet(x) BYTE4GETBYTE(x,3) #define mrcFileWriteMode_InfoPrint 1 #define mrcFileWriteMode_NoCallForStatData 2 extern mrcStatusType mrcFileWriteLowerResolution(mrcImage* mrc, char* filename, char* message, long sample , long mode);
画素のデータの取得
ピクセル単位での画素データの取得
extern double mrcPixelDataGet(mrcImage* mrc, mrcImageParaTypeReal x /* Pixel Unit */, mrcImageParaTypeReal y /* Pixel Unit */, mrcImageParaTypeReal z /* Pixel Unit */, double* data, mrcPixelDataType mode, mrcPixelDataHowToGet how);
オングストローム単位での画素データの取得
extern double mrcImageDataGetbyAU(mrcImage* mrc, mrcImageParaTypeReal x /* Angstrom */, mrcImageParaTypeReal y /* Angstrom */, mrcImageParaTypeReal z /* Angstrom */, double* data, mrcPixelDataType mode, mrcPixelDataHowToGet how);
インライン型関数の画素データの取得
extern inline double mrcPixelDataGetFloatImage(mrcImage* img, mrcImageParaTypeInteger x, mrcImageParaTypeInteger y, mrcImageParaTypeInteger z)
extern inline double mrcPixelDataGetFloatFTRe(mrcImage* fft, mrcImageParaTypeInteger x, mrcImageParaTypeInteger y, mrcImageParaTypeInteger z)
extern inline double mrcPixelDataGetFloatFTIm(mrcImage* fft, mrcImageParaTypeInteger x, mrcImageParaTypeInteger y, mrcImageParaTypeInteger z)
画素のデータの設定
画像情報の設定
extern void mrcInfoSet(mrcImage* mrc, FILE* fptIn, FILE* fptOut, mrcImageParaTypeInteger mode);
ピクセル単位の画素データの設定
extern void mrcImageDataSet(mrcImage* out, mrcImage* in, mrcImageParaTypeReal x, mrcImageParaTypeReal y, mrcImageParaTypeReal z, float phi, mrcImageDataSetMode mode);
オングストローム単位での画素データの設定
extern double mrcImageDataSetbyAU(mrcImage* mrc, mrcImageParaTypeReal x /* Angstrom */, mrcImageParaTypeReal y /* Angstrom */, mrcImageParaTypeReal z /* Angstrom */, double data, mrcPixelDataType mode);
インライン型関数の画素データの設定
extern double mrcPixelDataSet(mrcImage* mrc, mrcImageParaTypeReal x /* Pixel Unit */, mrcImageParaTypeReal y /* Pixel Unit */, mrcImageParaTypeReal z /* Pixel Unit */, double data, mrcPixelDataType mode); static inline void mrcPixelDataSetFloatImage(mrcImage* img, mrcImageParaTypeInteger x, mrcImageParaTypeInteger y, mrcImageParaTypeInteger z, double data) { img->ImageFloatImage[x + y*img->HeaderN.x + z*img->HeaderN.x*img->HeaderN.y] = (float)data; } static inline void mrcPixelDataSetFloatFTRe(mrcImage* fft, mrcImageParaTypeInteger x, mrcImageParaTypeInteger y, mrcImageParaTypeInteger z, double data) { fft->ImageFloatImage[2*(x + y*(fft->HeaderN.x/2+1) + z*(fft->HeaderN.x/2+1)*fft->HeaderN.y)] = (float)data; } static inline void mrcPixelDataSetFloatFTIm(mrcImage* fft, mrcImageParaTypeInteger x, mrcImageParaTypeInteger y, mrcImageParaTypeInteger z, double data) { fft->ImageFloatImage[2*(x + y*(fft->HeaderN.x/2+1) + z*(fft->HeaderN.x/2+1)*fft->HeaderN.y)+1] = (float)data; } extern void mrcStatDataSet(mrcImage* mrc, mrcImageParaTypeInteger mode);
画像処理一般のライブラリ
画像の演算(mrcImageOperation)
単項演算
加算
void lmrcImageAddedByReal(mrcImage* img, double d); img := img+d
減算
void lmrcImageMeanSubtraction(mrcImage* out, mrcImage* in) out := in - mean void lmrcImageMeanSubtraction(mrcImage* out, mrcImage* in) out := in - mean void lmrcImageSubtraction(mrcImage* out, mrcImage* in, double value) out := in - value
除算
void lmrcImageDevidedByReal(mrcImage* img, double d); img := img/d void lmrcImageDevidedByRealForOnlyPositive(mrcImage* img, double d); img := img/d, if img>=0, else img:=img; void lmrcImageDevidedByImage(mrcImage* img, mrcImage* i1, mrcImage* i2); img := i1/i2, 0 := 0/0, 0 := x/0
絶対値
void ImageAbs(mrcImage* out ,mrcImage* in); out := fabs(in);
それ以外
void lmrcImageSquare(mrcImage* in); in := in*in void lmrcImageRoot(mrcImage* in ,double eps); in := sqrt(in), if in>eps else in:=0
2項演算
連続加算
void lmrcImageAdd(mrcImage* img2, mrcImage* img1, long* number); img2 := img2 + img1, number++
2乗加算
void lmrcSQRImageAdd(mrcImage* img2, mrcImage* img1, long* number); img2 := img2 + img1*img1, number++
減算
void lmrcImageSubtractionCalc(mrcImage* out, mrcImage* in1, mrcImage* in2) out := in1 - in2 void llmrcImageSub(mrcImage*out ,mrcImage* in1 ,mrcImage* in2) out : = in1 - in2(上記と同じ機能)
畳み込み
extern void lmrcImageConvolution(mrcImage* out, mrcImage* in, mrcImage* kernel, int mode);
統計
void lmrcImagetMapForZero(mrcImage* t, mrcImage* avg, mrcImage* se); t := avg/se, at present, avg, se must be given. void lmrcImageSN(mrcImage* sn, mrcImage* avg, mrcImage* sig); sn := avg/sig, at present, avg, sig must be given. void mrcImageStdDev(mrcImage* sig, mrcImage* avg, mrcImage* avg2, long number); sig:= sqrt(avg2 - avg*avg), at present, avg, avg2 must be given. void mrcImageStdErr(mrcImage* se, mrcImage* avg, mrcImage* avg2, long number); se := sqrt(avg2 - avg*avg)/number, at present, avg, avg2 must be given.
エンディアン関連
extern int mrcImageByteSwap(mrcImage* mrc, long mode);
投影
extern void lmrcImageXSection(mrcImage* proj, mrcImage* img, double Y); extern void lmrcImageYSection(mrcImage* proj, mrcImage* img, double X); extern void lmrcImageXProjection(mrcImage* proj, mrcImage* img); extern void lmrcImageYProjection(mrcImage* proj, mrcImage* img); extern void lmrcImageZProjection(mrcImage* proj, mrcImage* img); extern void lmrcImageXYProjection(mrcImage* proj, mrcImage* img); extern void lmrcImageYZProjection(mrcImage* proj, mrcImage* img); extern void lmrcImageZXProjection(mrcImage* proj, mrcImage* img);
ヒストグラム
extern void lmrcImageHistgram(double** hist, unsigned long nlevel, mrcImage* img); extern void lmrcImageHistgram2(double** hist, double delta, mrcImage* img); extern void lmrcImageHistgram5(double** hist, unsigned long nlevel, double Low, double High, mrcImage* img); extern void lmrcImageHistgram3(double** hist, unsigned long nlevel, mrcImage* img); extern void lmrcImageHistgram4(double** hist, double delta, mrcImage* img); extern void lmrcImageHistgram6(double** hist, unsigned long nlevel, double Low, double High, mrcImage* img); extern void lmrcImageInformation(mrcImageInformation* info, mrcImage* img);
データコピー(mrcImageCopy)
extern void lmrcImageCopy(mrcImage* dst, mrcImage* src, mrcImageParaTypeRealCoord to);
画像モードのチェック
extern int IsFloatImage(mrcImage* mrc, char* message, mrcImageParaTypeInteger mode); extern int IsFloatFT(mrcImage* mrc, char* message, mrcImageParaTypeInteger mode); extern int IsImage(mrcImage* mrc, char* message, mrcImageParaTypeInteger mode); extern int IsFT(mrcImage* mrc, char* message, mrcImageParaTypeInteger mode); extern int IsFFT(mrcImage* mrc, char* message, mrcImageParaTypeInteger mode);
画像情報の設定
extern void mrcImageInfoSet(mrcImage* mrc, FILE* fptIn, FILE* fptOut, long mode);
mrcImageToIntImage
extern void lmrcImageToIntImage(mrcImage* dst, mrcImage* src, mrcImageParaTypeInteger mode);
FFT関連
extern void lfft2d(mrcImage* fft, mrcImage* img); extern void lmrcImageFFT(mrcImage* fft, mrcImage* img, long mode); extern void lmrcImageFFTFirstDimensionOnly(mrcImage* fft, mrcImage* img, long mode); #ifdef FFTW extern void lmrcImageFFTWFlagSet(unsigned int flag); #endif extern void lmrcFFTFGconj(mrcImage* out, mrcImage* in1, mrcImage* in2); extern void lmrcFFTFxG(mrcImage* out, mrcImage* in1, mrcImage* in2); extern void lmrcImageCorrelation(mrcImage* out, mrcImage* in, mrcImage* ref, long mode); extern void lmrcImageCorrelationModePrint(FILE* fpt); extern void lmrcImageCorrelationWithCTFCompensation(mrcImage* out, mrcImage* in, ctfInfo* inctf, mrcImage* ref, ctfInfo* refctf, long mode);
mrcError
extern mrcStatusType mrcError(char* message, char* inRtnName, char* callRtnName, mrcStatusType status); extern mrcStatusType mrcErrorMsg(FILE* fpt, char* message, char* inRtnName, char* callRtnName, mrcStatusType status);
lmrcFFTInfo
extern floatVector* lmrcFSInfoXAxisMag(mrcImage* fft); extern floatVector* lmrcFSInfoXAxisPhase(mrcImage* fft); extern floatVector* lmrcFSInfoYAxisMag(mrcImage* fft); extern floatVector* lmrcFSInfoYAxisPhase(mrcImage* fft);
extern floatVector* lmrcFSInfoScatteringAngularDistribution(mrcImage* fft); extern floatVector* lmrcFSInfoScatteringAngularDistributionAverage(mrcImage* fft); extern floatVector* lmrcFSInfoScatteringAngularDistributionSD(mrcImage* fft); extern floatVector* lmrcFSInfoSpacing(mrcImage* fft);
mrcImageCheckSameSize
extern void lmrcImageCheckSameSize(mrcImage* in, mrcImage* ref);
mrcImageCheckFFT
extern void lmrcImageCheckFFT(mrcImage* in, mrcImage* ref);
lmrcImageSmoothing
extern void lmrcImageSmoothingModeInfo(FILE* fpt); extern void lmrcImageSmoothing(mrcImage* dst, mrcImage* src, lmrcImageSmoothingInfo* info, long mode);
lmrcImageCVE
extern void lmrcImageCVE(mrcImage* dst, mrcImage* src, lmrcImageCVEInfo* info);
lmrcImageTfunction
extern void lmrcImageTfunction(mrcImage* dst, mrcImage* src, lmrcImageTfunctionInfo* info, long mode);
lmrcImageHighlighting
extern void lmrcImageHighlightingInfoPrint(FILE* fpt); extern void lmrcImageHighlighting(mrcImage* dst, mrcImage* src, lmrcImageHighlightInfo* info);
lmrcImageLowPassFilter
extern void lmrcImageLowPassFilter(mrcImage* dst, mrcImage* src, lmrcImageLowPassFilterInfo* info, long mode);
lmrcImageHighPassFilter
extern void lmrcImageHighPassFilter(mrcImage* dst, mrcImage* src, lmrcImageHighPassFilterInfo* info, long mode);
lmrcImageBandPassFilter
extern void lmrcImageBandPassFilter(mrcImage* dst, mrcImage* src, lmrcImageBandPassFilterInfo* info, long mode);
lmrcImageShrink
extern void lmrcImageShrink(mrcImage* dst, mrcImage* src, mrcImageParaTypeIntegerCoord S, long mode);
lmrcImageShift
extern void lmrcImageShift(mrcImage* out, mrcImage* in, mrcImageParaTypeRealCoord shift, mrcPixelDataHowToGet mode); extern void lmrcImageShiftFollowingGC(mrcImage* out, mrcImage* in, int mode);
mrcImagePrint
extern void mrcInfoPrint(FILE* fpt, mrcImage* mrc, mrcImageParaTypeInteger mode);
lmrcImageMultiCTFCompensation
extern void lmrcImageMultiCTFCompensationInfoInit(lmrcImageMultiCTFCompensationInfo* info); extern void lmrcImageMultiCTFCompensationInfoSet(lmrcImageMultiCTFCompensationInfo* info); extern void lmrcImageMultiCTFCompensation(mrcImage* dst, mrcImage* src, ctfInfo* ctf, long nfile, lmrcImageMultiCTFCompensationInfo info, long mode);
lmrcImageCTFObservation
extern void lmrcImageCTFObservation(mrcImage* dst, mrcImage* src, ctfInfo* ctf, long mode); extern void lmrcImagesCTFObservation(mrcImage* dst, mrcImage* src, ctfInfo* ctf, long mode);
lmrcImageCTFCompensation
extern void lmrcImageCTFCompensation(mrcImage* src, ctfInfo* ctf, long mode);
lmrcImageCTFSN
extern void lmrcImageCTFSN(mrcImage *outMrc, const mrcImage *inMrc, lmrcImageCTFSNInfo *info, long mode);
lmrcImageFourierPowerSpectrum
extern mrcImageFourierPowerSpectrum lmrcImageFourierPowerSpectrumMakeSpectrum(const mrcImage *inFFT, lmrcImageFourierPowerSpectrumInfo *info); extern void lmrcImageFourierPowerSpectrumDeleteSpectrum(mrcImageFourierPowerSpectrum theSpectrum); extern mrcImageFourierPowerSpectrum lmrcImageFourierPowerSpectrum(const mrcImage *inMrc, lmrcImageFourierPowerSpectrumInfo *info, long mode); extern void lmrcImageRotation3DModePrint(FILE* fpt); extern void lmrcImageRotation3DFollowingEulerAngle(mrcImage* dst, mrcImage* src, const char* m, mrcImageParaTypeReal rot1, mrcImageParaTypeReal rot2, mrcImageParaTypeReal rot3, mrcPixelDataHowToGet mode, int mode2); extern void lmrcImageRotation3DFollowingMatrix3D(mrcImage* dst, mrcImage* src, Matrix3D mat, mrcPixelDataHowToGet mode, int mode2); extern void lmrcImageRotation3DZXY(mrcImage* dst, mrcImage* src, mrcImageParaTypeReal rotx, mrcImageParaTypeReal roty, mrcImageParaTypeReal rotz, mrcPixelDataHowToGet mode, int mode2); extern void lmrcImageRotation2D(mrcImage* dst, mrcImage* src, double angle, mrcPixelDataHowToGet mode); extern void lmrcImageRotation2DPeriodicBoundary(mrcImage* dst, mrcImage* src, double angle, mrcPixelDataHowToGet mode);
mrcImageSectionGet
extern char* mrcImageSectionGetModeDescription(long mode); extern void mrcImageSectionGet(mrcImage* dst, mrcImage* src, mrcImageParaTypeReal z, long mode);
mrcImageSectionSet
extern void mrcImageSectionSet(mrcImage* dst, mrcImage* src, mrcImageParaTypeReal z, long mode); extern double lmrcImageMaxDataGet(mrcImage* src, long mode);