Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  504 1587 (N14E51 L= 98)  1270    264        1 ( 479/1487- 550/1678) -53.1 E55-E44 N08-N20
      2  371 1335 (S03E50 L= 98)   460     47        0 ( 351/1299- 391/1369)  90.0 E51-E48 S05-S01
      3  947 1287 (N08E11 L=137)   322     80        2 ( 914/1251- 984/1314)  59.0 E12-E09 N04-N09
      4 1326 1208 (N13W12 L=161)   250    114        0 (1280/1148-1364/1267) -26.6 W07-W15 N11-N15
      5  331 1019 (S22E45 L=104)   564    193        0 ( 269/ 971- 394/1073)  45.0 E49-E40 S26-S17
      6  826  768 (S24E05 L=144)  1244    393        3 ( 671/ 730-1029/ 833)  36.4 E14-W05 S29-S15
      7 1326  680 (S15W28 L=176)   479    224        1 (1218/ 633-1429/ 721)  45.0 W22-W32 S20-S10
      8 1371  589 (S19W34 L=182)   535    218        2 (1271/ 571-1473/ 605)  26.6 W27-W40 S21-S15
      9 1654  475 (S16W61 L=210)  3178    449        9 (1501/ 394-1754/ 566)  54.2 W52-W71 S25-S07
     10 1034  365 (S43W25 L=173)   455    102        1 (1000/ 341-1089/ 398)  18.4 W21-W28 S45-S40
     11 1223  257 (S43W45 L=194)   446     66        0 (1192/ 241-1259/ 273)  45.0 W42-W47 S44-S41
  

    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