在Angular 6应用程序中使用Leaflet-geotiff
我面临的问题无法在地图上显示任何.tif图像.我正在使用这是我使用它的样本: import * as geotiff from 'leaflet-geotiff/leaflet-geotiff'; import * as plotty from 'leaflet-geotiff/leaflet-geotiff-plotty'; export class MapComponent { ngOnInit() { this.map.on('load', () => { const options = { band: 0, name: 'Custom Mask', opacity: 1, renderer: new plotty.Plotty({ colorScale: '
0 2024-02-21
编程技术问答社区
在构建中包含 GeoTiff.js 和 Web Worker 时出现 Rollup/Vite 构建错误
我正在使用Web Worker在Geotiff文件上进行计算.我从CalculateRidge.ts>调用函数getRidgePoints,此功能创建一个工人,将消息发布给工人,然后等待工人的结果,最后在获得结果后终止该工人并返回结果. .工人正在用"等待"的功能包裹(我听说过comlink,但我找到了这个解决方案,而且效果很好). 工人具有addEventListener执行将进行计算的函数CreateRidge_Init.此函数加载一个geotiff文件,读取并根据从此文件读取的值进行计算. 我的问题是: 此代码无法构建.尝试构建代码时,我会遇到以下错误: Error [RollupError]: Invalid value "iife" for option "output.format" - UMD and IIFE output formats are not supported for code-splitting builds. at err
0 2023-12-03
编程技术问答社区
查找GeoTiff图像中每个像素的经/纬度坐标
我目前从Geotiff文件中有171 x 171的图像(尽管在其他情况下,我可能会有更大的图像).我的目标是将图像中的每个像素取转换为纬度/经度对. 我已经能够基于此堆叠式帖子将图像的角度转换为纬度/经度对:从Geotiff文件中获得纬度和经度.这篇文章很有帮助,因为我的原始坐标在UTM 15区. . 但是,我现在想将图像的所有像素转换为纬度,经度对,并将结果存储在相同维度的数字阵列中.因此,输出将是一个Numpy阵列,为171 x 171 x 2,numpy阵列的每个元素都是(经度,纬度)对的元组. 我在此上看到的最相关的帖子是 https://scriptndebug.wordpress.com/2014/11/11/11/24/latitudelongitude-each-pixel-pixel-using-python-pixel-using-python-and-gdal/.但是,该帖子建议基本上在每个像素上创建一个用于循环,然后转换为纬度,经度.有没有一种更有
2 2023-10-21
编程技术问答社区
用QPixmap::load( )读取GeoTiff时的TIFFReadDirectory警告
我有一个用qpixmap :: load()加载到qpixmap中的geotiff文件.我收到以下警告在控制台上打印了很多次.但是,直接使用libtiff可以打开没有警告的情况. 关于如何减轻QT中这些难看警告的任何想法? TIFFReadDirectory: Warning, foo: unknown field with tag 33550 (0x830e) encountered. TIFFReadDirectory: Warning, foo: unknown field with tag 33922 (0x8482) encountered. TIFFReadDirectory: Warning, foo: unknown field with tag 34735 (0x87af) encountered. TIFFReadDirectory: Warning, foo: unknown field with tag 34736 (0x87b0) encounter
46 2023-10-21
编程技术问答社区
-3.399999999999996e+38 NODATA 值和光栅 R
nodata 某些geotiffs的值( -3.3999999999999999996e+38 )被识别为 -Inf in r. 我使用栅格软件包读取文件并生成栅格堆栈. 使用 calc 计算栅格后, nodata 值将转换为 0 .. 是否有任何方法可以保留 nodata 标志值. 我认为这是一个浮点精度限制. 数据来自 v1.4/cmip5_30s.html 这是我的R代码. library(raster) library(rasterVis) tiffs
8 2023-08-13
编程技术问答社区
如何在python中写入/创建一个GeoTIFF RGB图像文件?
我有5个numpy阵列的形状NX,纽约 lons.shape = (nx,ny) lats.shape = (nx,ny) reds.shape = (nx,ny) greens.shape = (nx,ny) blues.shape = (nx,ny) 红色,绿色和蓝色阵列包含的值范围为0-255,LAT/LON阵列为纬度/经度像素坐标. 我的问题是如何将这些数据写给Geotiff? 我最终想使用基础图绘制图像. 这是我到目前为止的代码,但是我得到了一个庞大的Geotiff文件(〜500MB),它是空白的(只是一个黑色图像).另请注意,NX,NY = 8120,5416. from osgeo import gdal from osgeo import osr import numpy as np import h5py import os os.environ['GDAL_DATA'] = "/Users/andyprata/Library/Ent
10 2023-08-11
编程技术问答社区
R:直接从网页 url 读取 geotiff 数据(httr::GET 原始内容)
我想从服务器提供的Geotiff数据中创建一个rasterlayer.我将使用httr ::获取呼叫查询服务器此数据(按需提供数据,因此在应用程序中不会有.tif中的URL,而是查询URL). 在将此呼叫磁盘的结果作为geotiff文件编写后,很容易从磁盘上的结果geotiff文件创建rasterlayer: library(httr) library(raster) url
10 2023-07-27
编程技术问答社区
将 geotiff 转换为 datarame,并在列中保留纬度/纵深
我从此处下载了一个Geotiff: https usda.gov/research_and_science/crop_progress_gridded_layers/index.php (也可用文件:查看2021年的一个星期,我想将Geotiff转换为数据框架,因此我在Geotiff中的每个LAT/LON对有关联的值. 我尝试了: import rioxarray fl = 'data/cpc2021/corn/cpccorn2021/condition/cornCond21w24.tif' da = rioxarray.open_rasterio(fl, masked=True) df = da[0].to_pandas() df['y'] = df.index pd.melt(df, id_vars='y') 但是,这返回带有X和Y的数据框,似乎与LAT/LON不符.转换时如何添加(或保留)此信息? 期望lat/lon点处于连续的美国 编辑: 我找到
2 2023-07-25
编程技术问答社区
在R中把多个输入文件名作为输出文件名
我正在尝试使用栅格包装中的简单计算操作(从开尔文转换为摄氏度)的R中的数千张Geotiff图像. 计算工作起作用,但我需要原始输入文件名作为输出才能分配图像. 使用我的当前代码,仅列出列表中的第一个文件名,并且输出文件不断覆盖. 使用"计算"​​路径的文件列表: library(raster) files
12 2023-07-25
编程技术问答社区
为什么我用R-package Raster的 "image(x,y,z) "和 "plot (raster) "来绘制NetCDF图层时,得到的结果不同?
这可能很容易解决. 我想使用函数图(栅格)从NetCDF文件上绘制映射数据层.我不知道为什么我会得到栅格偏斜/偏移(我猜测问题是在转换中,分辨率?),如下图所示. 不正确的地图 如果我使用lat,lng,value矩阵使用函数映像(x,y,z ...),我将获得正确的显示,如下所示: 正确地图 这是我正在使用的R中的代码: library(ncdf) library(raster) # This is opening the NetCDF layer as a raster varRaster
10 2023-07-25
编程技术问答社区
使用R将光栅图像从等高线重新投影到纬线
嗨,我一直在尝试从equirectangular到EPSG:4326(latlon)重新投射栅格图像,问题是每次我在r上运行代码时,我都会在新图像上遇到错误的坐标;我不知道代码中的错误在哪里,我也与QGIS进行了相同的过程,而且我得到了相同的结果,这很奇怪,我有机会在Envi中进行相同的重新投影过程,结果是成功的,请帮助! ! a a class : RasterLayer dimensions : 1800, 1800, 3240000 (nrow, ncol, ncell) resolution : 1100, 1100 (x, y) extent : -988900, 991100, 1677577, 3657577 (xmin, xmax, ymin, ymax) coord. ref. : +proj=eqc +lat
10 2023-07-25
编程技术问答社区
R:编写RasterStack并保留图层名称
我有一个栅格堆栈,stk,由R中的三个栅格图像组成.这是一个简单的示例 # set up a raster stack with three layers > library(raster) > r r[] stk names(stk) [1] "layer.1" "layer.2" "layer.3" 我将名称分配给栅格层: # set layer names to "one", "two" and "three" > names(stk) names(stk) [1] "one" "two" "three" 当我将栅格堆栈写入Geotiff(多层)时,请使用: writeRaster(stk,"myStac
98 2023-07-25
编程技术问答社区
将hdf文件读入R,并将其转换为geoTIFF光栅
我正在尝试将MODIS 17数据文件读取到R中,操纵它们(裁剪等),然后将其保存为Geotiff.数据文件采用.hdf格式,似乎没有一种简单的方法将它们读入R. 与其他主题相比,那里没有很多建议,其中大部分已有几年历史了.其中一些也建议使用其他程序,但我想坚持使用R. 人们在r? 中使用哪种软件包来处理.hdf文件 解决方案 好,所以我的MODIS HDF文件是HDF4而不是HDF5格式.令人惊讶的是,很难发现这一点,MODIS在他们的网站上没有提及它,但是各种博客和堆栈交换帖子中有一些提示.最后,我不得不下载 hdfview 确定. r不做HDF4文件,几乎所有软件包(例如rgdal)仅支持HDF5文件.有一些有关下载驱动程序和从源头编译rgdal的文章,但这似乎很复杂,这些帖子适用于Mac或Unix,我正在使用Windows. . 基本上gdal_translate gdalUtils软件包是为想要在R中使用HDF4文件的任何人的保存宽限期所有例如通过
184 2023-07-18
编程技术问答社区
geotiff坐标转换的错误
在使用GDAL的应用程序中有错误,并且 crash 在OpenSuse中运行它的Geotiff图像中提取纬度和经度,而在Ubuntu中为我的同事们效果很好.对于GDAL-V3和GDAL-V2版本而言,错误是不同的,但是似乎问题在OGRCreateCoordinateTransformation对象创建中:在两种情况下返回NULL.请参阅下面的详细信息: 代码: QGeoCoordinate toGeoCoordinate(double* adGeotransform, OGRSpatialReference& srcRef, int x, int y) { double worldX = adGeotransform[0] + x * adGeotransform[1] + y * adGeotransform[2]; double worldY = adGeotransform[3] + x * adGeotransform[4] + y * adGeotr
14 2023-07-17
编程技术问答社区
怎么把png转换成geotiff?
我有一个png图像文件.我想将其转换为Geotiff.我安装了QGIS软件,但我无法使用它,也不知道如何查看图像.请帮我.有没有在线软件? 解决方案 这是一个很好的教程,带有有关如何在GUI中查看图像的逐步说明. /a> 其他解决方案 如果您想立即进行操作,则最快的选项将通过以下链接: https://mygeodata.cloud/cloud/converter/png-to-gtiff 不幸的是,该工具不会让您每月转换超过3个文件或总计超过5MB文件. 如果您想免费进行转换,这里有一个工具: geoconverter.hsr.ch/raster ,但质量比mygeodata云最差.但是,如果您只需要在边界上需要它,那么就可以了.
10 2023-07-16
编程技术问答社区
如何将png文件转换为geotiff文件
我有一组PNG文件和EPSG:每个文件边缘的3006个坐标.如何使用Python将这些PNG文件转换为Geotiff文件,以便TIFF文件包含Geo Metadata. 我想可以用rasterio lib来完成,但我不确定如何确切. 解决方案 找到了解决方案: dataset = rasterio.open(input_file_path, 'r') bands = [1, 2, 3] data = dataset.read(bands) transform = rasterio.transform.from_bounds(west, south, east, north, data.shape[1], data.shape[2]) crs = {'init': 'epsg:3006'} _, height, width = data.shape with rasterio.open(output_file_path, 'w', driver='GTiff',
6 2023-07-13
编程技术问答社区
Python 3:如何改变GDAL中的图像数据?
i有一个geotiff映像,其中包含一个带有8位表键的彩色表和单个栅格频带,并且使用lzw压缩,我加载了 gdal.Open .我还有一个包含24位RGB值的Numpy阵列(对于图像的模糊版本),对应于三个8位栅格频段.我需要将这三个栅格频段替换为当前在图像中的栅格频段,然后保存图像(如果可能的话,最好用作新文件).我该怎么做? 我想以RGB形式将数据保存在numpy数组中,因此我想以三个栅格频段而不是一个.我看到有一个 解决方案 我这样做的方式,只需创建一个带有新值的模板栅格的新副本...如果您想避免不惜一切代价获得副本,也可以覆盖.但是创建副本的错误较小,您可以保留原始版本. 此函数假设您有一个数组arr,该数组在第三维中有三个频段(So第二轴). import gdal def createRGB(template,arr,filename): '''Creates a copy of a 3-band raster with values from
38 2023-07-04
编程技术问答社区
如何在python中用gdal打开geotiff图片?
我正在尝试运行以下代码: from osgeo import gdal import sys # This allows GDAL to throw Python exceptions src_ds = gdal.Open("fused.tif") src_ds.show() 但是我收到以下错误: Traceback (most recent call last): File ".../gdalopen1.py", line 5, in module src_ds.show() AttributeError: 'Dataset' object has no attribute 'show' 为什么会发生这种情况? 解决方案 以下代码打开栅格文件,并将栅格读数读取到 numpy 数组. from osgeo import gdal ds = gdal.Open('input.tif', gdal.GA_ReadOnly) rb = ds.GetR
30 2023-07-03
编程技术问答社区