Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  578 1635 (N17E48 L=284)   775    202        2 ( 545/1563- 602/1724) -45.0 E55-E43 N12-N22
      2 1745 1478 (N45W42 L= 14)  1891    251        1 (1676/1383-1831/1563) -20.2 W33-W53 N41-N49
      3 1174 1414 (N19E03 L=330)   207     88        1 (1136/1380-1206/1435)  68.2 E03-W00 N16-N21
      4 1661 1352 (N32W32 L=  4)   860    156        0 (1613/1286-1698/1411) -29.1 W26-W35 N28-N34
      5 1591 1236 (N22W28 L=  0)   421    111        0 (1560/1183-1609/1287) -38.7 W24-W30 N18-N24
      6  956  931 (S14E01 L=331)   197     69        0 ( 926/ 904- 987/ 957)  63.4 E02-W00 S15-S11
      7  789  822 (S24E08 L=324)   286    151        0 ( 718/ 799- 858/ 840)  40.6 E11-E04 S27-S21
      8 1268  571 (S24W29 L=  1)   762    238        0 (1163/ 530-1383/ 604)  42.3 W24-W35 S29-S19
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y