遥感影像的保存需要无压缩形式的保存，普通JPG, TIFF等格式会压缩数据，改变原始数据值
所以应该采取二进制流保存的方式，用fopen, fwrite
后缀写成“.raw”（PS：单反相机采用的raw格式也是一种无压缩存储的形式）
例子：
str='january.raw';

fid=fopen(str,'wb');

fwrite(fid,temp1,'double');  %temp1是需要保存的矩阵，‘double’是数据精度

fclose(fid);



然鹅我经过试验还存在一个问题！
假设我们经过计算得到的temp1是一个3600行*7200列的数据，同时想要保存成这样的大小，如果直接输入temp1，竟然保存的的是temp1的转置！即7200行*3600列（这可能是matlab自身矩阵是按列优先，然后和写入文件的时候有差别造成的，猜的），所以在保存之前需要将矩阵转置一下：
temp1=temp1';


在envi中打开图像时，编辑头文件，输入行列号，数据精度，排列方式，就可以正常显示啦！


哈哈第一次写博客，请多多指教！

一、背景

在日吃昂的程序书写中，测试数据的来源分很多种。笔者一般喜欢用C原因呢的rand指令随机产生一些数据，但是在图像处理的程序中，C语言通常是通过读写图像文件。图像文件的产生可以用matlab，一般的流程为用matlab读入图像，在保存为二进制文件。二进制文件的保存通常不用添加文件后缀名，为了区分，也可以自行添加。



二、matlab写二进制文件



2.1 步骤

1）定义文件指针；2）写入



clear all;close all;clc;

image=imread('lena.bmp');
[rows,cols]=size(image);
vector(1,:)=reshape(image',[1,rows*cols]);

fid = fopen('data','wb');
fwrite(fid,vector(1,:),'single');



三、C语言读二进制文件



3.1知道文件中数据的size



void readDataLength(float *vector, int data_length)
{
    char flieName[50] = "data";  //文件路径
    FILE *fileIn;            //fp文件指针
    errno_t errIn;          //标记是否成功打开文件

    errIn = fopen_s(&fileIn, flieName, "rb"); // 成功ff=0,否则，ff!=0
    if (errIn != 0)
    {
        std::cout << "不能打开文件" << std::endl;
        return;
    }
    // 读入数据到内存中
    for (int i = 0; i < data_length; i++)
    {
        if (fread(&vector[i], sizeof(float), 1, fileIn) != 1){
            if (feof(fileIn)){
                fclose(fileIn);
                return;
            }
            printf("读出文件失败\n");
        }
        //fread(&vector[i], sizeof(float), 1, fileIn);
    }
    fclose(fileIn);
    std::cout << "读出文件成功！" << std::endl;
}



3.2不知道文件中数据的size

对于很多文件，出于一些原因，我们不能提前知道文件中数据的size，此时可以使用c语言中的文件结束标识符feof读取文件.



void readData(float* vector)
{
    char flieName[50] = "b";  //文件路径
    FILE *fileIn;            //fp文件指针
    errno_t errIn;          //标记是否成功打开文件
    int i = 0;
    errIn = fopen_s(&fileIn, flieName, "rb"); // 成功ff=0,否则，ff!=0
    if (errIn != 0)
    {
        std::cout << "不能打开文件" << std::endl;
        return;
    }

    // 读入数据到内存中
    //for (int i = 0; i < 100000; ++i)
    //{
    //  while (!feof(fileIn))
    //  {
    //      fread(&vector[i], sizeof(float), 1, fileIn);
    //      //if (i % 10 == 0){
    //      //  printf("%d", i);
    //      //}
    //  }
    //}
    while (!feof(fileIn)){
        fread(&vector[i], sizeof(float), 1, fileIn);
        i++;
        if (i % 1000 == 0){
            printf("%d\n", i);
        }
    }
    fclose(fileIn);
    std::cout << "读出文件成功！" << std::endl;
}



c语言二进制写文件



//二进制写文件
void saveData(float *vector, int data_length)
{
    char sXName[50] = "result";  // 文件路径

    FILE *fX;            //fp文件指针
    errno_t efX;            //标记是否成功打开文件

    efX = fopen_s(&fX, sXName, "wb");//成功ff=0,否则，ff!=0 
    if (efX != 0)
    {
        std::cout << "不能打开文件" << std::endl;
        return;
    }

    int eff = 0;
    for (int i = 0; i < data_length; i++)
    {
        eff = fwrite(&vector[i], sizeof(float), 1, fX);
        if (eff != 1)
        {
            std::cout << "文件写入错误" << std::endl;
            return;
        }
    }
    fclose(fX);

    std::cout << "已保存处理后的文件" << std::endl;
}

引言
PHP几乎很少处理二进制文件。但是便宜也完整的保留了这个功能。当你需要的时候，PHP自带的pack() & unpack()能能够极大地提供便利。下面我们从一个编程问题开始，讨论二进制文件的操作。
下文讨论gif文件，我们会编写一个函数，处理的内容跟GIF图像后缀无关。当然，我们也不打算尝试PHP的GD库。
gif文件头
不使用任何与图像处理相关的函数，为了解决这个问题，我们得从GIF文件本身获取数据。
与HTML、XML或其他文本格式文件不同，GIF文件和大多数其他图像格式是以二进制格式存储的。
大多数二进制文件的顶部都有一个头文件，它提供关于特定文件的元信息。我们可以使用这些信息来查找文件的类型和其他信息，比如GIF文件的高度和宽度。
下面显示了一个典型的原始GIF头文件，使用的是十六进制编辑器。
标题的详细描述如下。
因此，要检查图像文件是否是有效的GIF，我们需要检查文件的头3个字节，它有“GIF”标记，然后3个字节，它给出了版本号;“87a”或“89a”。
对于这样的需求，unpack()函数是必不可少的。在查看解决方案之前，我们快速查看一下unpack()函数本身。
使用unpack()函数
unpack()是pack()的补充——它根据指定的格式将二进制数据转换为关联数组。
它有点类似于sprintf，根据给定的格式转换字符串数据。
这两个函数（pack & unpack），允许我们根据指定的格式字符串读写二进制数据缓冲区。这使编程人员能够轻松地与用其他语言，或其他格式编写的程序交换数据。
举个小例子：
$data = unpack('C*', 'codediesel');
var_dump($data);
这会打印以下十进制代码的“codediesel”:
array
  1 => int 99
  2 => int 111
  3 => int 100
  4 => int 101
  5 => int 100
  6 => int 105
  7 => int 101
  8 => int 115
  9 => int 101
  10 => int 108
在上面的例子中，第一个参数是格式字符串，第二个参数是实际数据。
格式字符串指定应该如何解析数据参数。在本例中，格式“C”的第一部分指定我们应该将数据的第一个字符视为无符号字节。下一部分' * '告诉函数将前面指定的格式代码应用于所有剩余的字符。
乍一看上去挺混乱，但是慢慢梳理你会发现其中的规律。下一节我们提供一个具体的例子。
抓取头部数据
下面是使用unpack()函数解决GIF问题的方法。如果给定的文件是GIF格式的，则is_gif()函数将返回true。
function is_gif($image_file)
{
 
    /* Open the image file in binary mode */
    if(!$fp = fopen ($image_file, 'rb')) return 0;
 
    /* Read 20 bytes from the top of the file */
    if(!$data = fread ($fp, 20)) return 0;
 
    /* Create a format specifier */
    $header_format = 'A6version';  # Get the first 6 bytes

    /* Unpack the header data */
    $header = unpack ($header_format, $data);
 
    $ver = $header['version'];
 
    return ($ver == 'GIF87a' || $ver == 'GIF89a')? true : false;
 
}
 
/* Run our example */
echo is_gif("aboutus.gif");
需要注意的重要行是格式说明符。
' A6 '字符指定unpack()函数获取数据的前6个字节并将其解释为字符串。然后将检索到的数据存储在一个关联数组中，该数组的键名为“version”。
下面给出了另一个例子。
这将返回GIF文件的一些附加头数据，包括图像的宽度和高度。
function get_gif_header($image_file)
{
 
    /* Open the image file in binary mode */
    if(!$fp = fopen ($image_file, 'rb')) return 0;
 
    /* Read 20 bytes from the top of the file */
    if(!$data = fread ($fp, 20)) return 0;
 
    /* Create a format specifier */
    $header_format = 
            'A6Version/' . # Get the first 6 bytes
            'C2Width/' .   # Get the next 2 bytes
            'C2Height/' .  # Get the next 2 bytes
            'C1Flag/' .    # Get the next 1 byte
            '@11/' .       # Jump to the 12th byte
            'C1Aspect';    # Get the next 1 byte

    /* Unpack the header data */
    $header = unpack ($header_format, $data);
 
    $ver = $header['Version'];
 
    if($ver == 'GIF87a' || $ver == 'GIF89a') {
        return $header;
    } else {
        return 0;
    }
}
 
/* Run our example */
print_r(get_gif_header("aboutus.gif"));
上面的示例运行后打印以下内容。
Array
(
    [Version] => GIF89a
    [Width1] => 97
    [Width2] => 0
    [Height1] => 33
    [Height2] => 0
    [Flag] => 247
    [Aspect] => 0
)
下面我们将详细讨论格式说明符的工作方式。我将分解格式，给出每个字符的详细信息。
$header_format = 'A6Version/C2Width/C2Height/C1Flag/@11/C1Aspect';
更多格式选项可以在上图找到。
写在最后
我们展示的只是一个小小的例子，按照上图所揭示的规律，您可以任意组装成强大的解析函数。

引言
PHP几乎很少处理二进制文件。但是便宜也完整的保留了这个功能。当你需要的时候，PHP自带的pack() & unpack()能能够极大地提供便利。下面我们从一个编程问题开始，讨论二进制文件的操作。
下文讨论gif文件，我们会编写一个函数，处理的内容跟GIF图像后缀无关。当然，我们也不打算尝试PHP的GD库。
gif文件头
不使用任何与图像处理相关的函数，为了解决这个问题，我们得从GIF文件本身获取数据。
与HTML、XML或其他文本格式文件不同，GIF文件和大多数其他图像格式是以二进制格式存储的。
大多数二进制文件的顶部都有一个头文件，它提供关于特定文件的元信息。我们可以使用这些信息来查找文件的类型和其他信息，比如GIF文件的高度和宽度。
下面显示了一个典型的原始GIF头文件，使用的是十六进制编辑器。
标题的详细描述如下。
因此，要检查图像文件是否是有效的GIF，我们需要检查文件的头3个字节，它有“GIF”标记，然后3个字节，它给出了版本号;“87a”或“89a”。
对于这样的需求，unpack()函数是必不可少的。在查看解决方案之前，我们快速查看一下unpack()函数本身。
使用unpack()函数
unpack()是pack()的补充——它根据指定的格式将二进制数据转换为关联数组。
它有点类似于sprintf，根据给定的格式转换字符串数据。
这两个函数（pack & unpack），允许我们根据指定的格式字符串读写二进制数据缓冲区。这使编程人员能够轻松地与用其他语言，或其他格式编写的程序交换数据。
举个小例子：
$data = unpack('C*', 'codediesel');
var_dump($data);
这会打印以下十进制代码的“codediesel”:
array
  1 => int 99
  2 => int 111
  3 => int 100
  4 => int 101
  5 => int 100
  6 => int 105
  7 => int 101
  8 => int 115
  9 => int 101
  10 => int 108
在上面的例子中，第一个参数是格式字符串，第二个参数是实际数据。
格式字符串指定应该如何解析数据参数。在本例中，格式“C”的第一部分指定我们应该将数据的第一个字符视为无符号字节。下一部分' * '告诉函数将前面指定的格式代码应用于所有剩余的字符。
乍一看上去挺混乱，但是慢慢梳理你会发现其中的规律。下一节我们提供一个具体的例子。
抓取头部数据
下面是使用unpack()函数解决GIF问题的方法。如果给定的文件是GIF格式的，则is_gif()函数将返回true。
function is_gif($image_file)
{
 
    /* Open the image file in binary mode */
    if(!$fp = fopen ($image_file, 'rb')) return 0;
 
    /* Read 20 bytes from the top of the file */
    if(!$data = fread ($fp, 20)) return 0;
 
    /* Create a format specifier */
    $header_format = 'A6version';  # Get the first 6 bytes

    /* Unpack the header data */
    $header = unpack ($header_format, $data);
 
    $ver = $header['version'];
 
    return ($ver == 'GIF87a' || $ver == 'GIF89a')? true : false;
 
}
 
/* Run our example */
echo is_gif("aboutus.gif");
需要注意的重要行是格式说明符。
' A6 '字符指定unpack()函数获取数据的前6个字节并将其解释为字符串。然后将检索到的数据存储在一个关联数组中，该数组的键名为“version”。
下面给出了另一个例子。
这将返回GIF文件的一些附加头数据，包括图像的宽度和高度。
function get_gif_header($image_file)
{
 
    /* Open the image file in binary mode */
    if(!$fp = fopen ($image_file, 'rb')) return 0;
 
    /* Read 20 bytes from the top of the file */
    if(!$data = fread ($fp, 20)) return 0;
 
    /* Create a format specifier */
    $header_format = 
            'A6Version/' . # Get the first 6 bytes
            'C2Width/' .   # Get the next 2 bytes
            'C2Height/' .  # Get the next 2 bytes
            'C1Flag/' .    # Get the next 1 byte
            '@11/' .       # Jump to the 12th byte
            'C1Aspect';    # Get the next 1 byte

    /* Unpack the header data */
    $header = unpack ($header_format, $data);
 
    $ver = $header['Version'];
 
    if($ver == 'GIF87a' || $ver == 'GIF89a') {
        return $header;
    } else {
        return 0;
    }
}
 
/* Run our example */
print_r(get_gif_header("aboutus.gif"));
上面的示例运行后打印以下内容。
Array
(
    [Version] => GIF89a
    [Width1] => 97
    [Width2] => 0
    [Height1] => 33
    [Height2] => 0
    [Flag] => 247
    [Aspect] => 0
)
下面我们将详细讨论格式说明符的工作方式。我将分解格式，给出每个字符的详细信息。
$header_format = 'A6Version/C2Width/C2Height/C1Flag/@11/C1Aspect';
更多格式选项可以在上图找到。
写在最后
我们展示的只是一个小小的例子，按照上图所揭示的规律，您可以任意组装成强大的解析函数。