Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  420 1568 (N11E55 L=298)   469     69        0 ( 398/1533- 438/1595) -26.6 E58-E52 N09-N12
      2 1544 1518 (N37W21 L= 14)   371    106        0 (1504/1502-1590/1531)  35.5 W17-W24 N34-N39
      3 1761 1367 (N35W41 L= 34)   922    143        3 (1718/1317-1806/1411) -51.3 W37-W46 N31-N38
      4  300 1203 (S13E52 L=301)   444     80        2 ( 278/1164- 323/1240) -53.1 E53-E49 S14-S10
      5  172 1120 (S19E62 L=291)  1322    374        4 ( 100/1077- 252/1172)  -6.0 E73-E54 S21-S16
      6 1563 1110 (N12W29 L= 22)   286    102        0 (1538/1054-1576/1160) -56.3 W26-W30 N08-N14
      7  318 1069 (S20E47 L=306)  1205    363        0 ( 202/ 988- 433/1102)  22.8 E56-E37 S26-S17
      8 1834  923 (N11W56 L= 49)  1378    243        5 (1800/ 851-1893/ 988) -69.4 W52-W62 N06-N14
      9  646  915 (S22E20 L=333)   232     81        0 ( 597/ 904- 684/ 927)  21.8 E23-E17 S23-S20
     10  235  834 (S36E52 L=301)   801    292        1 ( 161/ 810- 313/ 850)  11.3 E59-E44 S38-S33
     11 1445  595 (S19W38 L= 31)   917    245        1 (1397/ 514-1464/ 708) -61.4 W32-W41 S24-S13
     12 1122  549 (S31W20 L= 13)  1906    535        3 ( 870/ 488-1380/ 608)  15.9 W02-W37 S35-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