Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

How to save region of interest from Band 3 Himawari Data as png image #272

Closed
xuesongle opened this issue Apr 27, 2018 · 1 comment
Closed

Comments

@xuesongle
Copy link

xuesongle commented Apr 27, 2018

Code Sample, a minimal, complete, and verifiable piece of code

# Your code here
from satpy.scene import Scene
from satpy.resample import get_area_def
from datetime import datetime
from glob import glob
from pyresample.utils import get_area_def

FILES = glob("/home/griffith/GIS/testData/*DAT")
scn = Scene(sensor='ahi',start_time=datetime(2018, 4, 14, 2, 00),end_time=datetime(2018, 4, 14, 2, 10),filenames=FILES,reader='ahi_hsd')
print(scn.available_dataset_names())
band = "B03"
scn.load([band])

area_id = 'japan_vol'
x_size = 392
y_size = 449
area_extent =  (-196090.1123778188, -222697.93829216747, 196090.11237782123, 226337.63024405806)
projection = '+proj=laea +lat_0=31.914 +lon_0=130.89 +ellps=WGS84'
description = "japan_vol"
proj_id = 'laea_130.89_31.914'


areadef = get_area_def(area_id, description, proj_id, projection,x_size, y_size, area_extent)
local_scene = scn.resample(areadef)
local_scene.save_datasets(writer='simple_image')

#### Problem description
If I use other bands, there is no error. Error only occurred in Band 3
I know error occurred at local_scene = scn.resample(areadef)
but not sure about how to fix it, could someone shed a light on this issue?

#### Expected Output

#### Actual Result, Traceback if applicable
/home/griffith/.conda/envs/satPy/lib/python3.6/site-packages/satpy/scene.py:128: UserWarning: Deprecated: Use 'filter_parameters' to filter loaded files by 'start_time', 'end_time', or 'area'.
  "'end_time', or 'area'.")
['B03']
Traceback (most recent call last):
  File "demo5_ch3.py", line 23, in <module>
    local_scene = scn.resample(areadef)
  File "/home/griffith/.conda/envs/satPy/lib/python3.6/site-packages/satpy/scene.py", line 596, in resample
    **resample_kwargs)
  File "/home/griffith/.conda/envs/satPy/lib/python3.6/site-packages/satpy/dataset.py", line 572, in resample
    new_data = resample(source_area, data, destination_area, **kwargs)
  File "/home/griffith/.conda/envs/satPy/lib/python3.6/site-packages/satpy/resample.py", line 521, in resample
    return resampler.resample(data, **kwargs)
  File "/home/griffith/.conda/envs/satPy/lib/python3.6/site-packages/satpy/resample.py", line 173, in resample
    cache_id = self.precompute(cache_dir=cache_dir, **kwargs)
  File "/home/griffith/.conda/envs/satPy/lib/python3.6/site-packages/satpy/resample.py", line 267, in precompute
    segments=segments)
  File "/home/griffith/.conda/envs/satPy/lib/python3.6/site-packages/pyresample/kd_tree.py", line 337, in get_neighbour_info
    nprocs=nprocs)
  File "/home/griffith/.conda/envs/satPy/lib/python3.6/site-packages/pyresample/kd_tree.py", line 409, in _get_valid_input_index
    source_lons, source_lats = source_geo_def.get_lonlats(nprocs=nprocs)
  File "/home/griffith/.conda/envs/satPy/lib/python3.6/site-packages/pyresample/geometry.py", line 1294, in get_lonlats
    data_slice=data_slice, dtype=dtype)
  File "/home/griffith/.conda/envs/satPy/lib/python3.6/site-packages/pyresample/geometry.py", line 1192, in get_proj_coords
    target_x, target_y = np.meshgrid(target_x, target_y)
  File "/home/griffith/.conda/envs/satPy/lib/python3.6/site-packages/numpy/lib/function_base.py", line 4698, in meshgrid
    output = [x.copy() for x in output]
  File "/home/griffith/.conda/envs/satPy/lib/python3.6/site-packages/numpy/lib/function_base.py", line 4698, in <listcomp>
    output = [x.copy() for x in output]
MemoryError
#### Versions of Python, package at hand and relevant dependencies
Python 3.6
Thank you for reporting an issue !
@mraspaud
Copy link
Member

I believe this is fixed, so I'm closing it for now. Feel free to reopen this if necessary.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants