Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  481 1421 (N37E37 L=130)   372     60        3 ( 448/1400- 519/1435)   0.0 E40-E33 N35-N38
      2  841 1362 (N29E08 L=159)   268     63        2 ( 812/1339- 878/1385) -14.0 E09-E05 N27-N30
      3  911 1209 (N19E05 L=162)   214     91        0 ( 901/1161- 923/1253)  80.5 E05-E03 N15-N21
      4 1643 1209 (N08W44 L=211)   307     64        1 (1623/1168-1661/1243)  53.1 W42-W46 N05-N09
      5  457 1113 (N17E37 L=130)  1112    397        1 ( 396/1000- 539/1189) -76.0 E42-E30 N10-N22
      6  288  906 (N06E53 L=114)   987    119        2 ( 249/ 843- 331/ 955)  35.0 E58-E48 N01-N08
      7 1582  393 (S47W41 L=209)  1432    179        1 (1515/ 358-1657/ 442)   7.1 W33-W49 S49-S44
      8  495  352 (S33E57 L=110)   565    106        1 ( 479/ 304- 509/ 403)  31.0 E66-E51 S37-S28
      9  604  289 (S39E51 L=116)  1357    148        3 ( 551/ 244- 657/ 332)  18.4 E65-E41 S43-S35
  

    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