Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1533 1461 (N19W44 L=135)   953    231        0 (1427/1430-1628/1497) -36.9 W38-W51 N15-N24
      2  239 1415 (N46E57 L= 33)   577    163        0 ( 210/1398- 273/1438) -14.0 E61-E53 N44-N47
      3 1797 1254 (N00W59 L=149)  2272    257        2 (1703/1129-1844/1366) -76.0 W53-W63 S07-N09
      4 1132 1214 (N16W11 L=101)   348     97        0 (1086/1184-1165/1259)  18.4 W07-W13 N14-N17
      5  209 1103 (N26E57 L= 33)   600    105        0 ( 180/1071- 245/1136) -29.7 E60-E53 N23-N28
      6 1825  956 (S17W53 L=143)   857    123        2 (1776/ 929-1855/ 990) -20.6 W48-W56 S19-S15
      7  733  570 (S14E28 L= 62)  1069    271        0 ( 688/ 464- 784/ 680)  42.7 E34-E21 S20-S07
      8 1582  522 (S39W28 L=118)  1682    355        5 (1522/ 469-1671/ 590) -37.6 W22-W36 S43-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