小编Kur*_*ehr的帖子

在ipython(0.10.2)中,我想运行bash shell样式(或glob.glob)扩展.我想做类似这样的事情,并且如果我在bash中那样扩展出来的话.这可能吗？

filenames = !ls *.sbet
run sbet.py $filenames
# fails because filenames is a list

# or the simpler case:
run sbet.py *.sbet
# failes with "*.sbet" being directly passed in without being expanded.

建议？谢谢!

我试图用jqgrid制作一个静态json文件的漂亮视图,并且是javascript的初学者.我创建了一个简化文件,但是当我查看日志时,jqGrid正在尝试在服务器上进行查询.是否可以将jqGrid与静态json一起使用？

我的测试文件:

http://schwehr.org/NGDC_DEMs.json

我在日志中看到了这一点:

http://schwehr.org/NGDC_DEMs.json?_search=false&nd=1344908237004&rows=20&page=1&sidx=&sord=asc

哪个不行.是否可以在浏览器中轻松完成所有这些操作？

该文件真的很无聊:

[
    {
        "title": "Elfin Cove MHHW",
        "demid": 4230,
    },
    {
        "title": "Astoria V2",
        "demid": 4090,
    },
    {
        "title": "South Padre Island",
        "demid": 3984,
    },
]

谢谢!-kurt

我有很好的numpy 2D浮动网格(例如> 20k x 10k单元格).我想计算网格中nodata值的连续组的数量.我可以在python中实现一个简单的pain fill方法,但似乎必须有方法在numpy或scipy中有效地做到这一点,但我没有找到一个明显的方法与ndimage.我当时认为必须有一种计算补丁的方法,我可以创建一个二进制网格并计算与nodata相对应的补丁.用scipy的fcluster可以做到这一点吗？

我试图在具有大的无数据值 (1e6) 的 2D 网格中找到所有整体的边界多边形。我已经得到了使用 scipy 标签工作的漏洞列表。如果不深入 gdal 的多边形化，是否有一种简单的方法来生成边界多边形？我看到有 matplotlib.pylab.contour，但这试图绘制一个我真的不想要的情节。关于如何为每个标签获取边界多边形的任何建议（如果可能的话，最好使用一种稍微简化多边形的方法）？我敢肯定我可以写一些东西来遍历每个标记孔的边界，但是有没有已经存在的东西？

from osgeo import gdal
from scipy import ndimage

dem_file = gdal.Open('dem.tif')
dem = dem.file.GetRasterBand(1).ReadAsArray()

# Get a binary image of the no-data regions.  The no-data value is large
bin = dem > 9e5

# Find all the wholes.  Anything with a label > 0.
labels, num_labels = ndimage.measurements.label(bin)
num_labels
1063

# The hole's label and size. Skip 0 as that label has all the valid data.
holes = [(label, sum(labels==label)) …