Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1134 1611 (N42W08 L=342)   713    149        2 (1111/1544-1158/1673)  78.7 W06-W10 N37-N47
      2 1673 1433 (N29W50 L= 23)   543    107        0 (1647/1404-1711/1456)   8.1 W46-W53 N26-N30
      3 1801 1423 (N27W63 L= 36)   837    181        2 (1765/1403-1843/1438) -11.3 W58-W68 N26-N28
      4  476 1384 (N25E40 L=293)   273     80        0 ( 450/1367- 506/1408) -11.3 E42-E37 N23-N26
      5 1207 1305 (N21W12 L=345)   656    119        1 (1168/1243-1248/1358) -50.2 W09-W14 N17-N24
      6  205 1142 (N08E62 L=271)  1928    340        3 ( 182/1050- 230/1241) -65.6 E65-E59 N02-N13
      7 1865 1182 (N12W65 L= 39)  1079    151        0 (1823/1129-1895/1234) -37.9 W60-W69 N08-N15
      8 1663  988 (N01W43 L= 17)  1136    193        2 (1590/ 943-1728/1041)  20.0 W37-W48 S01-N03
      9 1487  459 (S34W37 L= 10)   353     80        1 (1462/ 431-1507/ 498)  90.0 W34-W38 S36-S31
  

    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