Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  885 1656 (N34E26 L=182)  2237    518        4 ( 676/1551-1068/1757)  -4.9 E45-E10 N30-N37
      2 1157 1543 (N33E03 L=204)   378     81        1 (1124/1513-1198/1564)  59.0 E05-W00 N29-N35
      3  487 1492 (N15E46 L=161)   429     63        1 ( 467/1458- 505/1529) -68.2 E48-E44 N12-N17
      4  624 1412 (N13E33 L=174)   232     69        0 ( 596/1385- 646/1449)  78.7 E34-E32 N10-N15
      5  444 1383 (N08E46 L=162)   814    329        5 ( 334/1352- 584/1420)  24.0 E54-E36 N04-N12
      6  273 1306 (N00E59 L=148)   502    108        1 ( 242/1287- 304/1324)   9.5 E62-E55 S01-N01
      7 1799  969 (N11W54 L=261)  1432    229        3 (1761/ 877-1828/1056) -45.0 W48-W58 N06-N15
      8  288  998 (S18E49 L=158)   336     90        0 ( 252/ 979- 314/1014)  18.4 E52-E46 S19-S16
      9  757  981 (S10E15 L=192)   352    118        2 ( 693/ 957- 816/ 999)  23.2 E18-E11 S11-S07
     10  540  707 (S31E26 L=182)   254     49        0 ( 521/ 683- 561/ 734) -56.3 E27-E23 S32-S29
     11  994  620 (S27W07 L=214)   362     73        2 ( 951/ 598-1028/ 642)   0.0 W04-W09 S27-S24
  

    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