Zoe_左-CSDN博客

原创 python_MODIS HDF数据转为tif并拼接图像

首先是拼接，原始hdf数据是这样# -*- coding:utf-8 -*-import osimport arcpyfrom arcpy import envsourceDir=u'E:\' #中文路径输入targetDir=u'E:\' #英文路径输出arcpy.CheckOutExtension("Spatial")env.workspace = sourceDirarcpy.env.scratchWorkspace = sourceDirhdfList = arcpy

2020-08-07 19:09:15 7341 15

原创 python_拼接遥感影像Mosaic to New Raster

在使用Python GDAL将MODIS的HDF数据转换为tif之后，可使用Python中的ArcGISMosaic to New Raster工具来拼接多幅影像。或者也可继续使用GDAL 来拼接。数据文件名如：MOD11A2.A2000049.h23v04.006.2015058135050.tif# -*- coding: utf-8 -*-"""Created on Thu Au...

2019-08-10 09:58:45 4614 10

原创 python_批量将txt/csv坐标文件转为shp文件

# -*- coding: utf-8 -*-import arcpyfrom arcpy import envimport globimport osimport arcgisscriptinggp=arcgisscripting.create()spatial_ref=r'C:\Users\zoe\Documents\ArcGIS\unionCopy2.prj' #转换点文件需...

2019-05-11 09:41:02 8462 14

原创 python_计算格点数据降水变化斜率

import numpy as npimport pandas as pdimport codecspath=r'C:\try.csv' #读取格点数据CSVdata=pd.read_csv(path,encoding='gbk')outpath=r'C:\rate.csv'#新建存放斜率的CSVfile_csv = codecs.open(outpath,'w+','utf-8'...

2019-01-10 20:22:51 4231 8

原创 python_合并多年气象站点数据至新文件

# -*- coding: utf-8 -*-"""Created on Sat Oct 13 14:08:19 2018@author: Administrator"""#import reimport osimport pandas as pdinws=r'G:\QX_data\2002_2017\PRE_CSVaddTtile' #输入路径outws=r'G:\QX_d...

2018-10-13 16:46:28 2090

原创 python_MOD11A2 地表温度数据剔除缺省值和无效值

输入数据为使用MRT 拼接并转换投影后的tif文件# -*- coding: utf-8 -*-import arcpyfrom arcpy import envfrom arcpy.sa import *arcpy.CheckOutExtension("ImageAnalyst") #检查许可arcpy.CheckOutExtension("spatial") #检查许可...

2018-09-19 22:09:11 8030 30

原创 python_MOD13A3 NDVI数据剔除无效值

无效值为 -3000，使用setNull函数来执行# -*- coding: utf-8 -*-import arcpyfrom arcpy import envfrom arcpy.sa import *arcpy.CheckOutExtension("ImageAnalyst") #检查许可arcpy.CheckOutExtension("spatial") #检查许可wh...

2018-09-18 17:20:02 4955 10

原创 python_MOD13A3_NDVI 最大合成

如图，输入的文件路径下存放NDVI数据和 EVI数据，因此需要用到正则表达式匹配文件名。import arcpyfrom arcpy import envfrom arcpy.sa import *import rearcpy.CheckOutExtension("spatial") #检查模块许可# Set environment settingsenv.workspa...

2018-09-18 13:36:42 4984 3

原创 python_爬虫获取Modis数据下载页面的数据链接

# -*- coding: utf-8 -*-"""Created on Sat Jul 21 14:35:57 2018@author: Administrator"""import reimport urllibfrom bs4 import BeautifulSoupimport csvrawurl=r'https://e4ftl01.cr.usgs.gov/MOLT/...

2018-08-22 10:58:06 1873 4

原创 python_爬虫下载GLDAS数据

利用python requests以及正则表达式匹配来下载GLDAS数据，也可略微修改后用于Modis数据下载。# -*- coding: utf-8 -*-"""Created on Thu Aug 9 11:36:33 2018@author: Administrator"""# -*- coding: utf-8 -*-from bs4 import Beautifu...

2018-08-22 10:50:04 2034

原创 python_批量新建空白文件夹

# -*- coding: utf-8 -*-"""Created on Fri Oct 27 20:19:11 2017@author: Administrator"""import osbase = r'F:\QX_data\extract_process\TEM_persite\TEM_' #新建文件夹的路径i = 1981while i < 2011: #新建1...

2018-07-20 12:24:55 1534

原创 python_进行栅格计算_加法_批量处理

import arcpy from arcpy.sa import * arcpy.CheckOutExtension("spatial") arcpy.gp.overwriteOutput=1 #定义输出路径 outpath="e:/ccc/" outfilename="mean"#更改工作空间将，将"e:/bbb"更改为处理cdf文件时的输出路径arcpy.en...

2018-07-20 12:20:04 10599 11

原创 python 新建txt，并写入第一行title，将数据写入新txt

import pandas as pdimport os s = [] title=r'区站号纬度经度观测场拔海高度年月日平均气温日最高气温日最低气温平均气温质量控制码日最高气温质量控制码日最低气温质量控制码'+'\n'path = r'C:\Users\Administrator\

2017-10-21 12:34:46 6860 1

原创 python_使用需要的气象台站提取气象数据

# -*- coding: utf-8 -*-"""Created on Fri Oct 20 15:55:54 2017使用气象台站提取数据 @author: Administrator"""import pandas as pdimport os site=['58027','58040','58102','58221','58238','58251','58314','5...

2017-10-20 21:26:15 4368 8

转载 pandas 读取txt

转载于：http://blog.csdn.net/u011077672/article/details/509605801.原始文本格式：空格分隔的txt,例如1 2016-03-22 00:06:24.4463094 中文测试字符 2 2016-03-22 00:06:32.4565680 需要编辑encoding 3 2016-03-22 00:06:32.6835

2017-10-20 21:16:28 15119

转载 python—DataFrame的存储

转载于：http://blog.sina.com.cn/s/blog_4ddef8f80102vu4e.htmlPandas中的DataFrame数据既可以存储在SQL数据库中，也可以直接存储在CSV文件中。数据库1. dataframe.to_sql()函数将DataFrame数据存储到数据库中。name :数据表的名称， string, Name of

2017-10-20 20:59:04 25636 1

Python2.7_根据位置选择shp_备份.py

根据shp面要素的范围，批量选取在其范围内的shp点。只用Python2.7，批量实现该功能。代码中以及写了具体的注释，能够帮助大家理解代码含义。

2020-08-03

python2.7_批量TRMM3B43_单位mm/h转为mm/month

python2.7_批量TRMM3B43_单位mm/h转为mm/month 分闰年2月，不是2月，和不是闰年不是2月来转换单位

2018-09-16

python2.7栅格数据批量转换投影

python2.7栅格数据批量转换投影:ProjectRaster_management (in_raster, out_raster, out_coor_system, {resampling_type}, {cell_size}, {geographic_transform}, {Registration_Point}, {in_coor_system}) in_raster 输入栅格数据集。Mosaic Layer; Raster Layer out_raster 要创建的输出栅格数据集。以文件格式存储栅格数据集时，需要指定文件扩展名，具体如下：.bil - Esri BIL, .bip - Esri BIP, .bmp - BMP, .bsq - Esri BSQ, .dat - ENVI DAT,.gif - GIF,.img - ERDAS IMAGINE,.jpg - JPEG,.jp2 - JPEG 2000,.png - PNG,.tif - TIFF,无扩展名 - Esri Grid,以地理数据库形式存储栅格数据集时，不应向栅格数据集的名称添加文件扩展名。将栅格数据集存储到 JPEG 文件、JPEG 2000 文件、TIFF 文件或地理数据库时，可以指定压缩类型和压缩质量。 Raster Dataset out_coor_system 输入栅格待投影到的目标坐标系。默认值将基于“输出坐标系”环境设置进行设定。该参数的有效值是扩展名为 .prj 的文件。现有要素类、要素数据集、栅格目录（基本上包含了与坐标系相关的所有内容）。坐标系的字符串表示。要生成此类较长的字符串，可向模型构建器添加一个坐标系变量，并根据需要设置该变量的值，然后将模型导出到 Python 脚本。 Coordinate System resampling_type (可选) 要使用的重采样算法。默认设置为 NEAREST。 NEAREST —最邻近分配法 BILINEAR —双线性插值法 CUBIC —三次卷积插值法 MAJORITY —众数重采样法 NEAREST 和 MAJORITY 选项用于分类数据，如土地利用分类。NEAREST 选项是默认设置，因为它是最快的插值法，同时也因为它不会更改像元值。请勿对连续数据（如高程表面）使用 NEAREST 或 MAJORITY。BILINEAR 选项和 CUBIC 选项最适用于连续数据。不推荐对分类数据使用 BILINEAR 或者 CUBIC，因为像元值可能被更改。 cell_size (可选) 新栅格数据集的像元大小。默认像元大小为所选栅格数据集的像元大小。 Cell Size XY geographic_transform (可选) 在两个地理坐标系或基准面之间实现变换的方法。当输入和输出坐标系的基准面相同时，地理（坐标）变换为可选参数。如果输入和输出基准面不同，则必须指定地理（坐标）变换。有关各个受支持的地理（基准面）变换的详细信息，请参阅位于 ArcGIS 安装目录的 \Documentation 文件夹下的 geographic_transformations.pdf 文件。 Registration_Point(可选)用于对齐像素的 x 和 y 坐标（位于输出空间中）。配准点的工作原理与捕捉栅格的概念类似。通过配准点可指定用于定位输出像元的原点，而不是仅将输出捕捉到现有栅格像元。所有输出像元与该点之间必须间隔一个像元。该点的坐标不必位于一角，也不必落入栅格数据集中。捕捉栅格环境设置参数将优先于 Registration_Point 参数。因此，如果您要设置配准点，请确保尚未设置捕捉栅格。 in_coor_system (可选) 输入栅格数据集的坐标系。Coordinate System

2018-09-15

空空如也

TA创建的收藏夹 TA关注的收藏夹

TA关注的人

夏天的西瓜

原创 python_MODIS HDF数据转为tif并拼接图像

原创 python_拼接遥感影像Mosaic to New Raster

原创 python_批量将txt/csv坐标文件转为shp文件

原创 python_计算格点数据降水变化斜率

原创 python_合并多年气象站点数据至新文件

原创 python_MOD11A2 地表温度数据剔除缺省值和无效值

原创 python_MOD13A3 NDVI数据剔除无效值

原创 python_MOD13A3_NDVI 最大合成

原创 python_爬虫获取Modis数据下载页面的数据链接

原创 python_爬虫下载GLDAS数据

原创 python_批量新建空白文件夹

原创 python_进行栅格计算_加法_批量处理

原创 python 新建txt，并写入第一行title，将数据写入新txt

原创 python_使用需要的气象台站提取气象数据

转载 pandas 读取txt

转载 python—DataFrame的存储

Python2.7_根据位置选择shp_备份.py

python2.7_批量TRMM3B43_单位mm/h转为mm/month

python2.7栅格数据批量转换投影

python 转dbf文件为dataframe（CSV）

python_批量重采样

python2.7_批量掩膜提取数据

python批量定义投影

计量地理习题PPT

主成分分析 matlab

matlab 判定多边形顶点凹凸性

南京市地图 shp文件

python程序设计基础第八章习题代码

STELLA上机指导

全国地图 SHP

全国地图 SHP 带坐标系

计量地理课程

pythonpythonpython

空空如也

Python2.7_根据位置选择shp_备份.py

python2.7_批量TRMM3B43_单位mm/h转为mm/month

python2.7栅格数据批量转换投影

python 转dbf文件为dataframe（CSV）

python_批量重采样

python2.7_批量掩膜提取数据

python批量定义投影

计量地理习题PPT

主成分分析 matlab

matlab 判定多边形顶点凹凸性

南京市地图 shp文件

python程序设计基础第八章习题 代码

STELLA上机指导

全国地图 SHP

全国地图 SHP 带坐标系

计量地理课程

pythonpythonpython

空空如也

python程序设计基础第八章习题代码