C# 基于GDAL读取影像数据并用四叉树存储_c# gdal获取影像有效范围-程序员宅基地

原谅，这是刚开始学C#学的代码，所以比较凌乱，但是是可以运行的

using System;
using System.Collections.Generic;
using System.ComponentModel;
using System.Data;
using System.Drawing;
using System.Linq;
using System.Text;
using System.Windows.Forms;
using OSGeo.GDAL;
using System.IO;
using System.Data.OracleClient;
using System.Data.OleDb;

namespace GDALRead
{
    public struct point
    {
        public float x, y;
        public point(float x, float y)
        {
            this.x = x; this.y = y;
        }
    }
    public partial class Form1 : Form
    {
        int[,] _data = new int[16, 16];
        List<int[]> _QTree = new List<int[]>();
        int s = 0;
        public Form1()
        {
            InitializeComponent();
        }
        private void button1_Click(object sender, EventArgs e)
        {
            Gdal.AllRegister();
            // string fileName = @"E:\img\jiangxueying.img";
            string fileName = @"E:\img\11.tif";

            Dataset DEMdata = Gdal.Open(fileName, Access.GA_ReadOnly);
            if (DEMdata == null)
            {
                MessageBox.Show("不能打开：" + fileName);
                return;
            }

            int size = DEMdata.RasterCount;//影像大小
            int xsize = DEMdata.RasterXSize;//行
            int ysize = DEMdata.RasterYSize;//列

            Band DEMband = DEMdata.GetRasterBand(1);//波段
            double val;
            int hasval;
            DEMband.GetMinimum(out val, out hasval);
            if (hasval != 0)
            {
                string strmin = val.ToString();
            }
            DEMband.GetMaximum(out val, out hasval);
            if (hasval != 0)
            {
                string strmax = val.ToString();
            }

string strProjection = DEMdata.GetProjectionRef();//投影信息

            //输出图像的坐标和分辨率信息
            double[] adfGeoTransform = new double[6];
            DEMdata.GetGeoTransform(adfGeoTransform);

string strOrigin = adfGeoTransform[0] + "," + adfGeoTransform[3];
string strPixelsize = adfGeoTransform[1] + "," + adfGeoTransform[5];

double a;
DEMband.GetNoDataValue(out a, out hasval);

int[] DEM = new int[xsize * ysize];
DEMband.ReadRaster(0, 0, xsize, ysize, DEM, xsize, ysize, 0, 0);

            int t1 = 0;
            AllDEM = new int[xsize, ysize];
            for (int i = 0; i < xsize; i++)
            {
                for (int j = 0; j < ysize; j++)
                {

AllDEM[i, j] = DEM[t1];
t1++;

                }
            }
            //OpenFileDialog opF = new OpenFileDialog();
            //string[] lines = System.IO.File.ReadAllLines(@"E:\新建文本文档.txt", Encoding.Default);//读取文件
            //double[,] data = new double[lines.Length, lines[0].Split(',').Length];//将相应的数据存入数组中并赋值；
            //for (int i = 0; i < lines.Length; i++)
            //{//
            //for (int j = 0; j < lines[0].Split(',').Length; j++)
            //{
            //data[i, j] = double.Parse(lines[i].Split(',')[j]);
            //}
            //}
            int[,] element_data = element_array(AllDEM, 256);
            int[][,] child = new int[element_data.GetLength(0) * element_data.GetLength(1) / 256 / 256][,];
            int count = 0;
            string[][] all = new string[child.GetLength(0) * 256 * 256][];
            Point[][] p = depart(element_data, 256);
            string[][,] result = new string[child.GetLength(0)][,];
            string lll;
            for (int i = 0; i < child.GetLength(0); i++)
            {
                child[i] = index(element_data, p[i]);
                if (i == 132)
                    lll = "";
                result[i] = quadtree(child[i]);

                for (int j = 0; j < result[i].GetLength(0); j++)
                {
                    if (result[i][j, 0] == null && result[i][j, 1] == null || result[i][j, 1] == "0")
                        continue;
                    all[count] = new string[3];
                    all[count][0] = result[i][j, 0]; all[count][1] = result[i][j, 1]; all[count][2] = i.ToString();
                    count++;

}

}
MessageBox.Show("数据获取完毕，四叉分成功");

}

        计算morton值的方法
        //public int GetMorton(int row, int column)
        //{
        //    int i = 0;
        //    string c2 = "";
        //    string r = Convert.ToString(row, 2);
        //    string c = Convert.ToString(column, 2);
        //    if (r.Length > c.Length)
        //    {
        //        i = r.Length - c.Length;
        //        for (int a = 0; a < i; a++)
        //        {
        //            c2 = "0" + c2;
        //        }
        //        c = c2 + c;
        //    }
        //    else
        //    {
        //        i = c.Length - r.Length;
        //        for (int a = 0; a < i; a++)
        //        {
        //            c2 = "0" + c2;
        //        }
        //        r = c2 + r;
        //    }
        //    StringBuilder strBu = new StringBuilder();
        //    char[] rchar = r.ToCharArray();
        //    char[] cchar = c.ToCharArray();
        //    for (int j = 0; j < r.Length; j++)
        //    {
        //        strBu.Append(rchar[j]);
        //        strBu.Append(cchar[j]);
        //    }
        //    string result = strBu.ToString();
        //    int result10 = Convert.ToInt32(result, 2);
        //    return result10;
        //}

        //private void button2_Click(object sender, EventArgs e)
        //{ {
        //    //读取txt文件

本文链接：https://blog.csdn.net/sinat_30805761/article/details/79404904

原作者删帖不实内容删帖广告或垃圾文章投诉

智能推荐

c# 调用c++ lib静态库_c#调用lib-程序员宅基地

文章浏览阅读2w次，点赞7次，收藏51次。四个步骤1.创建C++ Win32项目动态库dll 2.在Win32项目动态库中添加外部依赖项 lib头文件和lib库3.导出C接口4.c#调用c++动态库开始你的表演...①创建一个空白的解决方案，在解决方案中添加 Visual C++ , Win32 项目空白解决方案的创建：添加Visual C++ , Win32 项目这......_c#调用lib

deepin/ubuntu安装苹方字体-程序员宅基地

文章浏览阅读4.6k次。苹方字体是苹果系统上的黑体，挺好看的。注重颜值的网站都会使用，例如知乎：font-family: -apple-system, BlinkMacSystemFont, Helvetica Neue, PingFang SC, Microsoft YaHei, Source Han Sans SC, Noto Sans CJK SC, W..._ubuntu pingfang

html表单常见操作汇总_html表单的处理程序有那些-程序员宅基地

文章浏览阅读159次。表单表单概述表单标签表单域按钮控件demo表单标签表单标签基本语法结构<form action="处理数据程序的url地址“ method=”get|post“ name="表单名称”></form><!--method将表单中的数据传送给服务器处理，get方式直接显示在url地址中，数据可以被缓存，且长度有限制；而post方式数据隐藏传输，_html表单的处理程序有那些

PHP设置谷歌验证器（Google Authenticator）实现操作二步验证_php otp 验证器-程序员宅基地

文章浏览阅读1.2k次。使用说明:开启Google的登陆二步验证（即Google Authenticator服务）后用户登陆时需要输入额外由手机客户端生成的一次性密码。实现Google Authenticator功能需要服务器端和客户端的支持。服务器端负责密钥的生成、验证一次性密码是否正确。客户端记录密钥后生成一次性密码。下载谷歌验证类库文件放到项目合适位置(我这边放在项目Vender下面)https://github.com/PHPGangsta/GoogleAuthenticatorPHP代码示例://引入谷_php otp 验证器

【Python】matplotlib.plot画图横坐标混乱及间隔处理_matplotlib更改横轴间距-程序员宅基地

文章浏览阅读4.3k次，点赞5次，收藏11次。matplotlib.plot画图横坐标混乱及间隔处理_matplotlib更改横轴间距

docker — 容器存储_docker 保存容器-程序员宅基地

文章浏览阅读2.2k次。①Storage driver 处理各镜像层及容器层的处理细节，实现了多层数据的堆叠，为用户提供了多层数据合并后的统一视图②所有 Storage driver 都使用可堆叠图像层和写时复制（CoW）策略③docker info 命令可查看当系统上的 storage driver主要用于测试目的，不建议用于生成环境。_docker 保存容器

随便推点

网络拓扑结构_网络拓扑csdn-程序员宅基地

文章浏览阅读834次，点赞27次，收藏13次。网络拓扑结构是指计算机网络中各组件（如计算机、服务器、打印机、路由器、交换机等设备）及其连接线路在物理布局或逻辑构型上的排列形式。这种布局不仅描述了设备间的实际物理连接方式，也决定了数据在网络中流动的路径和方式。不同的网络拓扑结构影响着网络的性能、可靠性、可扩展性及管理维护的难易程度。_网络拓扑csdn

JS重写Date函数，兼容IOS系统_date.prototype 将所有 ios-程序员宅基地

文章浏览阅读1.8k次，点赞5次，收藏8次。IOS系统Date的坑要创建一个指定时间的new Date对象时，通常的做法是：new Date("2020-09-21 11:11:00")这行代码在 PC 端和安卓端都是正常的，而在 iOS 端则会提示 Invalid Date 无效日期。在IOS年月日中间的横岗许换成斜杠，也就是new Date("2020/09/21 11:11:00")通常为了兼容IOS的这个坑，需要做一些额外的特殊处理，笔者在开发的时候经常会忘了兼容IOS系统。所以就想试着重写Date函数，一劳永逸，避免每次ne_date.prototype 将所有 ios