Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  538 1541 (N41E38 L= 65)   347     64        1 ( 519/1521- 564/1572) -31.0 E39-E34 N39-N42
      2 1049 1559 (N39W06 L=108)   566    111        1 ( 987/1537-1111/1572)  -6.3 W01-W10 N37-N40
      3 1239 1463 (N31W19 L=121)   345    121        0 (1196/1413-1273/1504) -63.4 W16-W20 N27-N34
      4  521 1301 (N24E34 L= 68)   610    120        0 ( 500/1231- 538/1355) -76.0 E35-E31 N19-N27
      5 1832 1269 (N12W65 L=167)  1232    172        1 (1801/1194-1862/1335)  63.4 W61-W71 N07-N15
      6 1788 1243 (N11W58 L=160)   514     68        0 (1766/1208-1808/1279) -63.4 W56-W60 N09-N13
      7  299 1223 (N20E53 L= 49)   368    105        0 ( 286/1170- 320/1272)  80.5 E53-E51 N16-N22
      8 1741 1131 (N05W51 L=153)   338     57        2 (1718/1107-1769/1151) -26.6 W48-W53 N03-N06
      9  190 1013 (N06E63 L= 39)   498     77        2 ( 158/ 997- 227/1032)  12.5 E68-E58 N05-N07
     10 1857  986 (S05W62 L=164)   644    129        0 (1812/ 937-1883/1029) -36.9 W56-W66 S08-S02
     11  202  886 (S01E63 L= 39)  1444    248        1 ( 180/ 773- 233/ 991)  65.2 E66-E57 S08-N05
     12 1783  775 (S17W55 L=157)  4617    612       11 (1706/ 602-1860/ 930) -84.8 W46-W65 S29-S07
     13 1402  797 (S12W23 L=125)   218     87        0 (1373/ 766-1432/ 830)  36.9 W20-W24 S13-S10
     14  413  651 (S16E46 L= 56)   388     83        0 ( 391/ 612- 431/ 690) -90.0 E47-E44 S18-S13
     15  892  465 (S31E13 L= 89)  1824    354        4 ( 769/ 383-1006/ 604)  36.9 E24-E03 S36-S21
     16 1144  521 (S28W05 L=107)   362    111        1 (1121/ 463-1159/ 567)  81.9 W03-W05 S32-S25
     17  439  496 (S27E49 L= 53)   534    128        1 ( 392/ 480- 503/ 512)   0.0 E54-E42 S28-S25
     18  727  422 (S33E27 L= 75)   324     74        2 ( 689/ 402- 768/ 436)  -9.5 E29-E23 S34-S32
  

    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