如何使用 libtiff.net 获取 GeoTiff 世界定位元数据答案

【问题标题】：How to get GeoTiff world-positioning metadata using libtiff.net如何使用 libtiff.net 获取 GeoTiff 世界定位元数据
【发布时间】：2021-02-02 09:24:43
【问题描述】：

我正在尝试从 ALOS 数据集 (JAXA) 中获取 GeoTiff 文件的边界框：

Tiff terrainTiff = Tiff.Open(@"Assets/Project/Heightmaps/" + "N046E007" + "/ALPSMLC30_" + "N046E007" + "_DSM.tif", "r");
                FieldValue[] modelPointTags = terrainTiff.GetField(TiffTag.GEOTIFF_MODELTIEPOINTTAG);
                foreach (FieldValue modelPointTag in modelPointTags)
                {
                    System.Type valueType = modelPointTag.Value.GetType();
                    Debug.Log(modelPointTag.Value.ToString());
                    Debug.Log(valueType);
                }

modelPointTags 中有 2 个值：System.Int32 和 System.Byte[]。如何继续阅读世界定位元数据？

【问题讨论】：

标签： geotiff libtiff libtiff.net

【解决方案1】：

用libtiff怎么样

 using (Tiff tiff = Tiff.Open(fileName, "r"))
        {
            //Image size
            int nWidth = tiff.GetField(TiffTag.IMAGEWIDTH)[0].ToInt();
            int nHeight = tiff.GetField(TiffTag.IMAGELENGTH)[0].ToInt();
            float[,] heightMap = new float[nWidth, nHeight];
            FieldValue[] modelPixelScaleTag = tiff.GetField(TiffTag.GEOTIFF_MODELPIXELSCALETAG);
            FieldValue[] modelTiePointTag = tiff.GetField(TiffTag.GEOTIFF_MODELTIEPOINTTAG);

            byte[] modelPixelScale = modelPixelScaleTag[1].GetBytes();
            double dW = BitConverter.ToDouble(modelPixelScale, 0);
            double dH = BitConverter.ToDouble(modelPixelScale, 8) * -1;

            byte[] modelTransformation = modelTiePointTag[1].GetBytes();
            double originLon = BitConverter.ToDouble(modelTransformation, 24);
            double originLat = BitConverter.ToDouble(modelTransformation, 32);

            double startW = originLon + dW / 2.0;
            double startH = originLat + dH / 2.0;

            FieldValue[] tileByteCountsTag = tiff.GetField(TiffTag.TILEBYTECOUNTS);
            long[] tileByteCounts = tileByteCountsTag[0].TolongArray();

            FieldValue[] bitsPerSampleTag = tiff.GetField(TiffTag.BITSPERSAMPLE);
            int bytesPerSample = bitsPerSampleTag[0].ToInt() / 8;

            FieldValue[] tilewtag = tiff.GetField(TiffTag.TILEWIDTH);
            FieldValue[] tilehtag = tiff.GetField(TiffTag.TILELENGTH);
            int tilew = tilewtag[0].ToInt();
            int tileh = tilehtag[0].ToInt();

            int tileWidthCount = nWidth / tilew;
            int remainingWidth = nWidth - tileWidthCount * tilew;
            if (remainingWidth > 0)
            {
                tileWidthCount++;
            }

            int tileHeightCount = nHeight / tileh;
            int remainingHeight = nHeight - tileHeightCount * tileh;
            if (remainingHeight > 0)
            {
                tileHeightCount++;
            }

            int tileSize = tiff.TileSize();
            for (int iw = 0; iw < nWidth; iw += tilew)
            {
                for (int ih = 0; ih < nHeight; ih += tileh)
                {
                    byte[] buffer = new byte[tileSize];
                    tiff.ReadTile(buffer, 0, iw, ih, 0, 0);
                    for (int itw = 0; itw < tilew; itw++)
                    {
                        int iwhm = ih + itw;
                        if (iwhm > nWidth - 1)
                        {
                            break;
                        }
                        for (int ith = 0; ith < tileh; ith++)
                        {
                            int iyhm = iw + ith;
                            if (iyhm > nHeight - 1)
                            {
                                break;
                            }
                            heightMap[iwhm, iyhm] =
                              BitConverter.ToSingle(buffer, (itw * tileh + ith) * 4);

                            Console.WriteLine(heightMap[itw, ith]);
                        }
                    }
                }
            }
        }

【讨论】：

【解决方案2】：

我能够使用来自https://build-failed.blogspot.com/2014/12/processing-geotiff-files-in-net-without.html 的说明读取边界框

【讨论】：