Questions tagged «floating-point»

5
是否将栅格完全加载到numpy数组中?
我一直在尝试检查我在DEM栅格上的过滤器以进行模式识别,并且它总是导致缺少最后一行和最后一列(如..20)。我已经尝试了PIL库,图像加载。然后用numpy。输出是相同的。 我以为我的循环出了点问题,当检查数组中的值时(只是用ArcCatalog中的Identification选择像素),我意识到像素值未加载到数组中。 因此,只需打开,放入数组并保存数组中的图像即可: a=numpy.array(Image.open(inraster)) #raster is .tif Float32, size 561x253 newIm=Image.new(Im.mode, Im.size) Image.fromarray(a).save(outraster) 导致删除最后的行和列。抱歉,无法发布图片 任何人都可以帮助理解为什么?并建议一些解决方案? 编辑: 因此,在伙计们的帮助下,我成功地将小型栅格加载到numpy数组中,但是当图像较大时,我开始出现错误。我想这是关于numpy数组的限制,因此数组会像这样自动重塑或变形……所以ex: Traceback (most recent call last): File "<pyshell#36>", line 1, in <module> ima=numpy.array(inDs.GetRasterBand(1).ReadAsArray()) File "C:\Python25\lib\site-packages\osgeo\gdal.py", line 835, in ReadAsArray buf_xsize, buf_ysize, buf_obj ) File "C:\Python25\lib\site-packages\osgeo\gdal_array.py", line 140, in BandReadAsArray ar = numpy.reshape(ar, [buf_ysize,buf_xsize]) File …

2
在GIS中处理空间数据时的精度
是否有基础研究或入门研究在工作时检查和比较空间数据的精度 具有不同的数据输入精度,例如具有1、2,...小数位? 有不同的浮点实现(float,double)? 与赤道附近的数据相比,赤道附近的数据? 用隧道距离,大圆距,文森特,鲍林,兰伯特计算的地理距离? 到目前为止,我发现的所有工作都指出,这些都是错误的来源,但没有给出人们可以预期的确切错误界限。

1
在ArcMap定义查询中引用NaN或±∞(特殊浮点值)?
我偶然发现了ArcMap如何向用户显示特殊的浮点值。 +∞(正无穷大)显示为 1.#INF –∞(负无穷大)应该显示为-1.#INF—我尚未验证这一点。 NaN(不是数字)显示为右对齐<Null>-请勿与left-aligned混淆<Null>,后者表示NULL(缺少值): (顺便说一下,在字段计算器中获取唯一值根本不会列出NaN。) 但是我还没有发现如何编写图层定义查询以基于这些特殊值选择行: ColumnName IS NULL 只会选择常规的NULL值,而不会选择NaN。 ColumnName = 1.#INF 因语法无效而被拒绝。 有谁知道如何做到这一点? 用于将1.#INF值存储到表字段的C#ArcObjects代码片段(基本概念): 按照要求。由于我不再工作了,因此以下代码不是我使用的真实代码,现在无法测试,但是它会产生如上屏幕截图所示的效果: ITable table = …; int doubleFieldIndex = table.FindField(…); IRow row = table.CreateRow(); row.Value[doubleFieldIndex] = double.PositiveInfinity; row.Store();
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.