Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  815 1855 (N48E46 L=339)  1402    242        2 ( 778/1808- 881/1900) -14.0 E54-E38 N42-N51
      2  597 1465 (N15E37 L=347)  2039    513        1 ( 464/1334- 740/1603) -20.1 E54-E24 N09-N21
      3 1582 1307 (N27W30 L= 55)   670    193        0 (1516/1227-1637/1403) -39.3 W24-W35 N23-N32
      4 1240 1279 (N17W07 L= 32)   433     88        1 (1201/1245-1279/1324) -31.0 W04-W10 N15-N19
      5  618 1142 (S04E27 L=358)   677    138        2 ( 553/1102- 687/1178)  41.2 E30-E22 S07-N00
      6 1258 1121 (N08W12 L= 36)   271     31        0 (1238/1097-1280/1146) -33.7 W10-W13 N07-N09
      7  860  955 (S10E08 L= 17)   731    238        1 ( 744/ 924- 982/ 988)  11.3 E15-W00 S11-S07
      8  236  812 (S31E53 L=332)  1264    285        3 ( 176/ 746- 301/ 879)  39.8 E60-E46 S36-S26
      9 1320  525 (S25W32 L= 56)   440     80        0 (1279/ 509-1364/ 541)  18.4 W28-W34 S26-S23
     10  887  502 (S37W04 L= 29)   469    139        0 ( 825/ 477- 969/ 526)  29.1 W00-W09 S40-S33
  

    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