Magnetic cutoff filtering

[imogenagle]

Figure out why the hmi data is filtered the way it is.

[samaloney]

Can you link to the code here?

[imogenagle]

CHIMERApy/chimerapy/chimera.py

Lines 292 to 325 in 910d3c4

	pos = np.zeros((len(poslin[0]), 2), dtype=np.uint)
	pos[:, 0] = np.array((poslin[0] - (s[0] / 2)) * convermul + (s[1] / 2), dtype=np.uint)
	pos[:, 1] = np.array((poslin[1] - (s[0] / 2)) * convermul + (s[1] / 2), dtype=np.uint)
	npix = list(
	np.histogram(
	datm[pos[:, 0], pos[:, 1]],
	bins=np.arange(
	np.round(np.min(datm[pos[:, 0], pos[:, 1]])) - 0.5,
	np.round(np.max(datm[pos[:, 0], pos[:, 1]])) + 0.6,
	1,
	),
	)
	)
	npix[0][np.where(npix[0] == 0)] = 1
	npix[1] = npix[1][:-1] + 0.5
	wh1 = np.where(npix[1] > 0)
	wh2 = np.where(npix[1] < 0)
	# =====magnetic cut offs dependent on area=========
	if (
	np.absolute((np.sum(npix[0][wh1]) - np.sum(npix[0][wh2])) / np.sqrt(np.sum(npix[0])))
	<= 10
	and arcar < 9000
	):
	continue
	if (
	np.absolute(np.mean(datm[pos[:, 0], pos[:, 1]])) < garr[int(cent[0]), int(cent[1])]
	and arcar < 40000
	):
	continue
	iarr[poslin] = ident