Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1756 1492 (N17W69 L=159)  1078    117        3 (1711/1460-1807/1524) -33.7 W63-W75 N14-N20
      2 1548 1438 (N20W45 L=135)   572     86        2 (1520/1391-1572/1472)  33.7 W42-W48 N17-N22
      3 1726 1420 (N15W60 L=150)   812     88        1 (1687/1380-1761/1459) -45.0 W56-W64 N11-N16
      4  624 1402 (N37E22 L= 69)   231     86        0 ( 605/1362- 647/1441)  45.0 E23-E19 N34-N39
      5 1088 1170 (N14W06 L= 97)   218    103        1 (1048/1141-1130/1209) -59.0 W04-W08 N11-N17
      6  155 1143 (N25E70 L= 20)  1218    133        2 ( 123/1087- 189/1203) -26.6 E78-E63 N21-N28
      7 1596  941 (S09W34 L=125)   238     89        0 (1583/ 897-1610/ 988) -90.0 W33-W35 S12-S06
      8 1221  886 (S05W09 L= 99)   491    210        1 (1164/ 836-1286/ 979) -60.9 W06-W12 S08-N01
      9 1837  890 (S19W57 L=147)   901    115        1 (1805/ 847-1886/ 928)   5.2 W52-W63 S21-S16
     10 1513  895 (S10W28 L=118)   268     68        0 (1491/ 868-1537/ 926) -68.2 W26-W29 S12-S07
     11  350  838 (N05E48 L= 42)   412    110        0 ( 311/ 803- 397/ 874)  18.4 E52-E43 N03-N06
     12 1285  775 (S13W11 L=102)   409     97        3 (1244/ 740-1321/ 806)  36.9 W08-W13 S14-S11
     13  325  597 (S10E59 L= 31)  1676    194        3 ( 282/ 527- 389/ 673) -68.2 E62-E55 S14-S04
     14  714  548 (S17E29 L= 61)   362     94        1 ( 695/ 496- 731/ 596)  71.6 E30-E27 S20-S14
     15  293  538 (S13E67 L= 23)   970     49        0 ( 272/ 500- 311/ 572)  59.0 E70-E63 S16-S11
     16 1441  530 (S33W19 L=109)   711    117        0 (1382/ 510-1494/ 566)   0.0 W14-W23 S34-S30
     17  485  462 (S20E50 L= 40)   515     97        0 ( 459/ 413- 510/ 514)  37.9 E54-E45 S23-S16
     18 1007  455 (S29E12 L= 78)   460     92        0 ( 988/ 408-1029/ 505)  56.3 E14-E10 S32-S25
     19  594  328 (S31E50 L= 41)  1145    177        1 ( 517/ 300- 681/ 349)  -3.4 E57-E40 S33-S29
  

    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