- 博客(181)
- 资源 (86)
- 收藏
- 关注
RSS订阅原创 6s 大气校正处理 ( HJ2A2B CCD1 ，CCD2 ) 分块提速
【代码】6s 大气校正处理 ( HJ2A2B CCD1 ，CCD2 ) 分块提速。
2025-04-30 15:03:28
287
原创 gdal 正射校正
transformer_options.append("RPC_HEIGHT=0") # 默认椭球面高度。output_resolution (float): 输出分辨率（单位与 dst_srs 一致）dem_path (str): DEM 文件路径（可选，不提供则使用椭球面高度）input_path (str): 输入影像路径（需含 RPC 元数据）dem_path="srtm_30m.tif", # 使用 SRTM DEM。(2) 无 DEM 校正（仅 RPC）。"""批量正射校正"""
2025-04-30 13:54:39
226
原创 6s 大气校正处理 波谱曲线重采样
点击OK，Enter Containing Output Wavelength界面选择包含输出波长的文件，即自定义的CSV文件，打开，设置SCAII文件相关参数，即单位，波长所在列等（若波长范围重采样前后数量级一致可不选择波长单位参数，默认Unknown；b = np.arange(450, 520+ 2.5, 2.5) # 蓝光波段。g = np.arange(520, 600+ 2.5, 2.5) # 绿光波段。# 获取所有TXT文件（排除临时文件和已处理文件）# 带进度条的批量处理。
2025-04-30 09:52:11
901
原创 可切换对比图组建swipe
<!DOCTYPE html><html lang="zh-CN"><head> <meta charset="UTF-8"> <meta name="viewport" content="width=device-width, initial-scale=1.0"> <title>可切换对比图组件</title> <script src="https://cdn.tailwindcss.com"></script> <style> .compar
2025-04-29 14:41:21
606
原创 6s 大气校正处理 ( HJ2A2B CCD1 ，CCD2 )
#!usr/bin/env python# -*- coding:utf-8 -*-import globimport osimport sysimport tarfileimport jsonimport numpy as npimport gdalimport mathimport timeimport xml.dom.minidomfrom tqdm import tqdmfrom Py6S import *import argparsefrom base import
2025-04-27 15:17:53
902
原创 6s 大气校正处理 (HJ 1A CCD1 ，CCD2 )
outDataset = Driver.Create(OutRCname, cols, rows, 1, gdal.GDT_Int32) # 改为5个波段。map_projection = dom.find('.//mapProjection').text # 如'UTM'earth_model = dom.find('.//earthModel').text # 如'WGS_84'zone = dom.find('.//zone').text # 如'53N'# 生成WKT格式（备选方案）
2025-04-27 12:28:27
970
原创 XYZ tile 瓦片生成原理
if tile_data.ndim == 3 and tile_data.shape == 3: # 检查形状:ml-citation{ref="6" data="citationList"}tile_data = np.ma.getdata(tile_data).copy() # 处理掩码数据:ml-citation{ref="1" data="citationList"}if tile_data.shape == (3, 256, 256): # 检查是否为(C,H,W)
2025-04-15 13:07:45
860
原创 拼接右移算法
raise FileNotFoundError(f"无法打开文件: {envi_file_path}")# 获取第一个波段（通常ENVI文件是多波段的，但这里只读取一个作为示例）# 遍历数组中的每个元素（除了最后一行和最后一列）# 检查右侧是否为1。# 检查下方是否为1。# 初始化一个空列表来存储满足条件的索引。# 检查当前元素是否为1。# 读取波段数据到NumPy数组。# 检查数据集是否成功打开。# 获取数组的行数和列数。# 打印满足条件的索引。# 打开ENVI文件。
2025-02-26 12:03:38
330
原创 文件分片上传 python
console.error(`上传分片 ${chunkIndex} 失败: ${error.message}, 重试次数: ${chunkRetries[chunkIndex]}/${MAX_RETRIES}`);console.error(`分片 ${chunkIndex} 上传失败，已达到最大重试次数`);url: '/upload_filesliceprocess', // 替换为你的上传URL。// 1MB 分片大小。// 上传下一个分片。// 处理上传失败的情况，比如停止上传或通知用户。
2025-02-18 16:31:30
863
原创 python 影像裁切程序
dtypes = [src.dtypes[i] for i in range(src.count)] # 确保数据类型列表与波段数量匹配。# 假设所有波段的数据类型都是相同的，并且你想要保持与输入影像相同的数据类型。# 将rasterio的bounds转换为shapely的box对象。# 将面要素的边界转换为shapely的box对象。# 为裁剪后的影像设置输出路径和文件名。# 检查面要素的边界是否与影像的边界相交。# 创建输出文件，并写入裁剪后的影像数据。# 使用这个数据类型打开输出文件。
2025-01-13 15:05:42
401
原创 Geoserver 自动发布矢量和Tiff栅格
store_name='TIFF1209'#相对于data_dir路径。一份带投影 3857坐标系的GeoTiff文件。file_name='TIFF1209.tif'#文件名。layer_name='TIFF1209'#图层名。准备 一份带投影 ３８５７ 坐标系shp矢量数据。二自动发布Tiff栅格。
2025-01-07 15:49:53
915
原创 Envi 将扩大10000的反射率还原成float
pro bandmathdivide10000 COMPILE_OPT idl2 envi,/restore_base_save_files envi_batch_init ;image_file=DIALOG_PICKFILE(PATH=’d:\’, $ ; TITLE=’Select DICOM Patient File’, FILTER=’*.tif’) ;...
2019-03-06 14:19:19
2556
转载 mongodb 安装和
MongoDB设置访问权限、设置用户MongoDB已经使用很长一段时间了，基于MongoDB的数据存储也一直没有使用到权限访问（MongoDB默认设置为无权限访问限制），今天特地花了一点时间研究了一下，研究成果如下：注：研究成果基于Windows平台MongoDB在本机安装部署好后1. 输入命令：show dbs，你会发现它内置有两个数据库，一个名为admin，一个名为loca
2016-06-19 15:53:42
826
转载 mongodb 编译
考虑到mongodb使用了boost库源码，参考mongodb官方文档后，下载编译boost版本是1.42(时间为2010-2-2)或更新版本: boost版本1.42： http://sourceforge.net/projects/boost/files/boost/1.42.0/boost_1_42_0.zip/download
2016-06-18 10:55:39
687
原创 Spring 依赖配置
使用方法1： //在执行此实例化的时候就会完成所有注入 ClassPathXmlApplicationContext ctx = new ClassPathXmlApplicationContext("beans.xml"); UserService service = (UserService)ctx.getBean("userService"); 使用方法2：
2016-06-16 13:44:09
295
翻译 ssh 定时任务
-- 刷新EHR用户角色 --> "synRoleId" class="org.springframework.scheduling.quartz.MethodInvokingJobDetailFactoryBean"> "targetObject"> "webService"/> "targetMe
2016-06-15 15:23:28
1118
原创 定时 备份mysql 数据库
需求： 1，每天4点备份mysql数据； 2，为节省空间，删除超过3个月的所有备份数据； 3，删除超过7天的备份数据，保留3个月里的 10号 20号 30号的备份数据； #创建shell文件vim backup_mysql.sh mysqldump -uroot -p123456 --all-databases >/data/dbdata
2016-06-15 15:16:23
354
原创 mysql latin1 utf8 转换
在生成环境中，会经常碰到Mysql字符的设置问题，个人也总结下。如果在应用开始阶段没有正确的设置字符集，在运行一段时间以后才发现存在不能满足要求需要调整，那么就需要进行字符集的修改。字符集的修改不能直接通过alter dataabase character set ***; 或者 alter table tablename character set ***; 命令进行，这两个命令都没有更新已
2016-06-02 13:09:29
580
原创 mysqll 字符集合分析
MySQL中涉及的几个字符集character-set-server/default-character-set：服务器字符集，默认情况下所采用的。character-set-database：数据库字符集。character-set-table：数据库表字符集。优先级依次增加。所以一般情况下只需要设置character-set-server，而在创建数据库和表时不特别指定字符集，
2016-06-02 13:04:38
403
原创 comet 应用
pushlet/comet, 比websocket有名， 也应用的广。 比如webqq就是用erlang的mochiweb做的长连接; python社区有个国产的eurasia框架，也是专注comet,也做了很多年;原来那个挺出名的随机聊天网站omegle也是基于长连接的， 开始用的twisted，后来改成了gevent; ape-project是个专注长连接的应用,用js处理逻辑的，早于no
2016-05-23 19:56:14
314
转载 Mysql字符集设置
最近，在项目组使用的mysql数据库中，插入数据出现乱码，关于这个问题做了下总结，我们从最基本的地方说起，到错误产生的深层次原因和解决办法。基本概念• 字符(Character)是指人类语言中最小的表义符号。例如’A'、’B'等；• 给定一系列字符，对每个字符赋予一个数值，用数值来代表对应的字符，这一数值就是字符的编码(Encoding)。例如，我们给字符’A'赋予数值0，给字符’B'
2016-03-11 15:09:14
644
原创 VC modbus 读
static unsigned char auchCRCHi[] = {0x00, 0xC1, 0x81, 0x40, 0x01, 0xC0, 0x80, 0x41, 0x01, 0xC0, 0x80, 0x41, 0x00, 0xC1, 0x81,0x40, 0x01, 0xC0, 0x80, 0x41, 0x00, 0xC1, 0x81, 0x40, 0x00, 0xC1, 0x81,
2016-03-09 08:54:28
2025
原创 VB Modbus Rtu 读
1） 模块文件：modCRC，其中包含了CRC校验的函数。'data 待校验的数组名称'no 数组中元素个数'btLoCRC 算出的CRC高字节'btHiCRC 算出的CRC低字节Public Function CalCRC16Fast(data() As Byte, no As Integer, btLoCRC As Byte,
2016-03-09 08:52:24
2133
原创 Modbus Ascii Rtu 区别与联系
ModBus-ASCII协议和RTU协议的比较通过比较可以看到，ASCII协议和RTU协议相比拥有开始和结束标记，因此在进行程序处理时能更加方便，而且由于传输的都是可见的ASCII字符，所以进行调试时就更加的直观，另外它的LRC校验也比较容易。但是因为它传输的都是可见的ASCII字符，RTU传输的数据每一个字节ASCII都要用两个字节来传输，比如RTU传输一个十六进制数0xF9,ASC
2016-03-09 08:46:27
1790
转载 开源 android
主要介绍那些不错个性化的View，包括ListView、ActionBar、Menu、ViewPager、Gallery、GridView、ImageView、ProgressBar及其他如Dialog、Toast、EditText、TableView、Activity Animation等等。 一、ListView android-pulltorefresh
2015-12-07 09:57:49
546
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人