Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  492 1656 (N36E52 L=162)   609    104        0 ( 454/1637- 538/1674)  26.6 E55-E48 N34-N37
      2  596 1653 (N37E43 L=171)   469     90        1 ( 570/1621- 623/1674) -38.7 E46-E40 N34-N38
      3  763 1302 (N16E20 L=195)   314     97        0 ( 727/1267- 815/1323)  31.0 E21-E16 N13-N17
      4  207 1227 (N05E64 L=150)   599     65        0 ( 190/1190- 230/1263)  68.2 E66-E62 N02-N07
      5  313 1148 (N01E51 L=163)   391     68        2 ( 292/1117- 338/1176) -33.7 E53-E48 N00-N02
      6  647 1135 (N04E25 L=189)   301     59        1 ( 623/1101- 669/1168)  68.2 E26-E23 N01-N06
      7 1088  894 (S06W05 L=219)   763    234        0 (1036/ 771-1133/ 989) -58.4 W01-W09 S13-N00
      8  723  871 (S11E17 L=197)   429    142        3 ( 647/ 856- 789/ 892)   6.3 E22-E13 S12-S09
      9 1178  685 (S18W13 L=227)   234     53        0 (1149/ 667-1210/ 701)  26.6 W11-W15 S19-S17
  

    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