9527的代码梦 2022-08-16 14:33 采纳率: 100%
浏览 34
已结题

c#如何在libtiff里边拿到高程数据

需要做一个上传tiff文件的功能,那么就需要把tiff的数据拿出来,由于要用c#和libtiff库来完成,网上找了很多资料,导出来的高程数据还是有问题
using (Tiff tiff = Tiff.Open(path, "r"))
            {
                int height = tiff.GetField(TiffTag.IMAGELENGTH)[0].ToInt();
                int width = tiff.GetField(TiffTag.IMAGEWIDTH)[0].ToInt();

                FieldValue[] modelPixelScaleTag = tiff.GetField((TiffTag)33550);
                FieldValue[] modelTiepointTag = tiff.GetField((TiffTag)33922);

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

                byte[] modelTransformation = modelTiepointTag[1].GetBytes();
                double originLon = BitConverter.ToDouble(modelTransformation, 24);//x
                double originLat = BitConverter.ToDouble(modelTransformation, 32);//y

                int bitsPerPixel = tiff.GetField(TiffTag.BITSPERSAMPLE)[0].ToInt();

                double startLat = originLat + (pixelSizeY / 2.0);
                double startLon = originLon + (pixelSizeX / 2.0);

                var scanline = new byte[tiff.ScanlineSize()];

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

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

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

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

heightMap 不知道这个数据到底拿的有没有问题

img

现在基本上都通了就是高程数据拿的有问题
我想把一个tiff文件里的高程数据存到一个float里边,然后用它生成坐标生成地图

img

  • 写回答

1条回答 默认 最新

  • 石岩Maple 2022-08-16 14:47
    关注

    srtm数据?

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论

报告相同问题?

问题事件

  • 系统已结题 8月24日
  • 已采纳回答 8月16日
  • 创建了问题 8月16日

悬赏问题

  • ¥15 winform的chart曲线生成时有凸起
  • ¥15 msix packaging tool打包问题
  • ¥15 finalshell节点的搭建代码和那个端口代码教程
  • ¥15 用hfss做微带贴片阵列天线的时候分析设置有问题
  • ¥15 Centos / PETSc / PETGEM
  • ¥15 centos7.9 IPv6端口telnet和端口监控问题
  • ¥20 完全没有学习过GAN,看了CSDN的一篇文章,里面有代码但是完全不知道如何操作
  • ¥15 使用ue5插件narrative时如何切换关卡也保存叙事任务记录
  • ¥20 海浪数据 南海地区海况数据,波浪数据
  • ¥20 软件测试决策法疑问求解答