Previous topic

asciiread

Next topic

binwrite

This Page

binread

mipylib.dataset.midata.binread(fn, dim, datatype=None, skip=0, byteorder='little_endian')

Read data array from a binary file.

Parameters:
  • fn – (string) The binary file name for data reading.
  • dim – (list) Dimensions.
  • datatype – (string) Data type string.
  • skip – (int) Skip bytes number.
  • byteorder – (string) Byte order. little_endian or big_endian.
Returns:

(MIArray) Data array

DMSP SSM/I-SSMIS Pathfinder Daily EASE-Grid Brightness Temperatures, Version 2:

#Data source: https://nsidc.org/data/docs/daac/nsidc0032_ssmi_ease_tbs.gd.html
#DMSP SSM/I-SSMIS Pathfinder Daily EASE-Grid Brightness Temperatures, Version 2

fn = 'D:/Temp/binary/ID2-F17-SH2015035D-V2.91V'

#Set projection
proj = projinfo(proj='laea', lat_0=-90, lon_0=0)
xn = 1441             #Column number
yn = 1441             #Row number
dx = 12.5 * 1000      #X resolution (m)
dy = 12.5 * 1000      #Y resolution (m)
x0 = -9000 * 1000 + dx / 2     #Lower left x corner
y0 = -9000 * 1000 + dy / 2    #Lower left y corner
x = arange1(x0, xn, dx)
y = arange1(y0, yn, dy)

#Read data
data = binread(fn, [yn, xn], 'short')
data = data[::-1,:].astype('float')
data = data / 10
data[data==0] = nan   #0 is missing value

#Plot
axesm(projinfo=proj, gridline=True, griddx=30, griddy=30)
lworld = shaperead('D:/Temp/Map/country1.shp')
geoshow(lworld, edgecolor='k')
layer = imshowm(x, y, data, 20, cmap='BlAqGrYeOrRe', proj=proj)
colorbar(layer)
title('Brightness Temperatures')
../../../../_images/ssmis_bt_laea1.png