Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1411 1745 (N48W51 L= 48)   890    108        1 (1361/1708-1457/1785) -56.3 W48-W54 N44-N51
      2  497 1383 (N32E36 L=321)  1161    324        4 ( 382/1267- 592/1443) -29.7 E48-E27 N24-N36
      3 1717 1105 (N02W48 L= 45)   295     52        2 (1695/1078-1737/1134)  45.0 W45-W50 N00-N03
      4  640  875 (S00E25 L=331)   276    141        3 ( 582/ 835- 689/ 932) -33.7 E28-E22 S03-N03
      5 1665  850 (S13W42 L= 38)   291     68        0 (1633/ 831-1692/ 878)  21.8 W38-W43 S13-S11
      6 1827  814 (S18W59 L= 55)   678    134        0 (1794/ 777-1851/ 856)  53.1 W54-W62 S20-S15
      7  612  750 (S08E29 L=328)   380     71        1 ( 582/ 722- 655/ 771)   9.5 E31-E25 S09-S06
      8  533  666 (S13E36 L=320)   402    116        1 ( 512/ 618- 555/ 714)  74.1 E37-E34 S16-S09
      9 1234  665 (S19W10 L=  6)   384    117        2 (1203/ 628-1272/ 710) -78.7 W07-W12 S21-S16
     10 1432  579 (S27W24 L= 20)  1445    410        2 (1301/ 480-1598/ 701)  25.6 W13-W36 S33-S21
     11  713  647 (S15E24 L=333)   452     95        0 ( 674/ 619- 743/ 693)  45.0 E26-E21 S16-S12
     12  884  638 (S17E13 L=344)   332    110        0 ( 831/ 614- 930/ 658)  15.9 E16-E09 S18-S16
     13 1731  626 (S29W51 L= 48)   647     64        1 (1688/ 609-1772/ 646) -11.3 W46-W56 S30-S27
     14 1000  571 (S23E06 L=350)   571    207        0 ( 953/ 498-1075/ 638)  38.7 E10-W00 S27-S19
     15  598  591 (S18E33 L=323)   259     59        1 ( 567/ 573- 628/ 607) -14.0 E35-E31 S19-S17
     16 1521  485 (S36W32 L= 29)   598    148        0 (1448/ 451-1576/ 527)  20.0 W26-W37 S37-S33
     17 1010  370 (S38E09 L=348)   461    134        0 ( 957/ 354-1074/ 384) -12.5 E12-E03 S39-S36
     18 1064  299 (S45E06 L=350)  1848    270        3 ( 993/ 210-1156/ 377)  45.0 E13-W02 S54-S39
  

    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