pyhton 离散点数据转栅格img

最新推荐文章于 2022-12-04 12:35:58 发布

天道酬勤_业道酬信

最新推荐文章于 2022-12-04 12:35:58 发布

阅读量3.4k

点赞数 3

CC 4.0 BY-SA版权

文章标签： python 转img arcgis

本文链接：https://blog.youkuaiyun.com/qq_16645423/article/details/80359207

这个Python类实现了将离散点数据或格点数据转换为ArcGIS的栅格图像（img）。它包括两个主要方法：station_to_img用于处理站点数据，grid_to_img用于处理格点数据。转换过程中涉及的空间范围可以通过设置上下左右经纬度来定义，并且提供了重采样的选项。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

# -*-coding:utf-8 -*-

import arcpy
import os
from dbcommon.views import OperateDB
from arcpy import env
import shutil
import numpy as np

class DataToImg(object):
    """
    :站点和格点数据转img图层公用算法
    """
    def __init__(self):
        #img文件临时存储路径
        self.OUTPATH = 'F:/arcgis_temp/raster/'
        #img文件存储路径
        self.IMGPATH = 'F:/arcgis_temp/img/'
        #环境变量设置   上下左右经纬度
        self.x_min = 104.250
        self.y_min = 35.200
        self.x_max = 107.700
        self.y_max = 39.400

    
    def station_to_img(self,data,out_path,file_name):
        """
        :站点数据转换为img图层
        :param  data  :数据    二维数据[[经度，纬度，数据],………]
                out_path :输出路径 例如：FGFP/
                file_name:文件名称
        """
        try:
            env.overwriteOutput = True
            #工作空间
            env.workspace = self.OUTPATH + out_path
            #环境变量设置   上下左右经纬度
            env.extent = arcpy.Extent( self.x_min, self.y_min, 
                                       self.x_max, self.y_max
                                      )
            #环境变量设置   上下左右经纬度
            img_path = self.discretePoint2raster(data, out_path, file_name)
            #如果路径不存在，新建目录
            if not os.path.exists(self.IMGPATH+out_path):
                os.makedirs(self.IMGPATH+out_path)
            file_name = self.IMGPATH+out_path+"/"+file_name+".img"
            #将文件永久保存
            shutil.copyfile(img_path,file_name)
            #删除所有的临时文件
            self.remove_all_file(self.OUTPATH+out_path)
            return file_name
        except:
            raise
        
    def grid_to_img(self,data,out_path,file_name,usecols=None):
        """
        :站点数据转换为img图层
        :param  data     :数据    二维数据[[经度，纬度，数据],………]
                out_path :输出路径 例如：FGFP/
                file_name:文件名称
                usecols  :文件参与绘图的列序号，判断传入数据是否为文件  例如（0,1,2）