Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  504 1676 (N23E52 L= 41)   705    140        1 ( 473/1626- 529/1729) -53.1 E55-E49 N19-N27
      2  907 1513 (N19E18 L= 75)   465    197        0 ( 864/1425- 938/1586) -45.0 E22-E13 N14-N23
      3  258 1321 (S02E59 L= 33)  3150    483        8 ( 178/1099- 332/1505) -80.9 E69-E55 S16-N09
      4 1233 1337 (N15W06 L= 98)   202     85        1 (1219/1296-1247/1382) -68.2 W04-W07 N13-N18
      5 1725 1304 (N26W40 L=133)   367     70        0 (1701/1265-1741/1335) -45.0 W38-W41 N23-N28
      6  772 1168 (S04E17 L= 75)   322    141        0 ( 716/1107- 816/1232)  69.4 E19-E15 S08-N00
      7 1814 1075 (N15W51 L=144)   356     68        0 (1793/1039-1833/1105) -45.0 W48-W53 N13-N16
      8  297 1054 (S18E49 L= 44)   311     60        0 ( 286/1021- 308/1091)  90.0 E50-E47 S19-S15
      9  742  835 (S23E13 L= 79)   206     90        0 ( 699/ 823- 795/ 847)  15.9 E16-E09 S24-S21
  

    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