Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  913 1905 (N48E52 L=219)   776     54        1 ( 882/1873- 947/1924)  53.1 E56-E46 N44-N50
      2  800 1845 (N38E51 L=221)  1295    138        3 ( 767/1781- 848/1892) -50.7 E58-E44 N32-N43
      3 1096 1631 (N32E15 L=256)   666     96        1 (1064/1571-1121/1681) -63.4 E17-E12 N27-N35
      4 1357 1525 (N33W05 L=276)  1305    334        2 (1276/1373-1429/1657) -43.2 E04-W13 N24-N40
      5 1473 1112 (N12W23 L=294)  1537    413        1 (1410/ 940-1524/1278) -74.1 W17-W26 N00-N21
      6  492 1154 (S13E35 L=236)   488    100        1 ( 471/1099- 520/1208) -56.3 E37-E32 S15-S09
      7  745 1008 (S14E15 L=256)   254    133        0 ( 715/ 958- 791/1046)  71.6 E16-E13 S17-S11
      8  642  961 (S20E20 L=251)   240     67        2 ( 615/ 930- 670/ 987) -26.6 E22-E18 S20-S18
      9  989  413 (S42W21 L=292)   660    211        0 ( 885/ 393-1083/ 433)  18.4 W12-W27 S45-S39
     10 1308  404 (S31W41 L=313)   669    125        1 (1249/ 386-1370/ 431)  18.4 W36-W45 S33-S29
     11 1557  297 (S27W74 L=346)   867     98        0 (1510/ 281-1611/ 314)  26.6 W68-W80 S29-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