Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1221 1896 (N60W63 L= 32)   952    121        1 (1168/1870-1277/1921)  -1.9 W50-W82 N58-N61
      2  559 1680 (N59E35 L=295)  1198    199        3 ( 480/1637- 641/1716)   8.1 E45-E23 N55-N61
      3 1124 1660 (N45W23 L=352)   834    194        0 (1056/1610-1207/1687)   0.0 W15-W30 N41-N46
      4 1270 1627 (N39W32 L=  2)   324     65        0 (1239/1603-1296/1649) -63.4 W30-W34 N36-N41
      5 1516 1606 (N32W51 L= 20)   688     62        2 (1480/1575-1552/1631) -15.9 W48-W55 N29-N33
      6 1336 1492 (N29W31 L=360)   290     74        0 (1299/1473-1376/1511) -45.0 W28-W33 N26-N30
      7  824 1375 (N31E07 L=322)  1922    568        4 ( 599/1295-1081/1471) -28.0 E24-W08 N22-N41
      8 1704 1430 (N17W57 L= 27)   676    113        1 (1648/1410-1753/1458) -10.3 W51-W62 N14-N18
      9 1428 1254 (N13W29 L=358)   302     82        0 (1389/1219-1452/1288) -68.2 W27-W30 N10-N15
     10  613 1072 (N15E25 L=304)   265     99        0 ( 578/1032- 657/1100)  26.6 E27-E21 N13-N16
     11  175 1039 (N17E65 L=264)   729    158        2 ( 146/1006- 215/1098)  90.0 E69-E59 N15-N21
     12  398  987 (N13E42 L=287)  1193    314        1 ( 293/ 865- 494/1061) -40.9 E50-E35 N04-N17
     13 1371  835 (S10W18 L=347)   577    172        0 (1336/ 754-1392/ 903)  53.1 W14-W20 S14-S06
     14 1579  778 (S18W32 L=  1)   279     67        1 (1549/ 764-1610/ 795) -21.8 W29-W34 S19-S16
     15  841  644 (S13E18 L=312)   932    360        1 ( 779/ 575- 923/ 707) -26.6 E20-E12 S17-S09
     16 1050  443 (S30E09 L=320)  1248    291        5 ( 954/ 371-1153/ 564)  30.5 E17-W00 S34-S23
     17  833  319 (S36E29 L=300)   309     38        0 ( 813/ 300- 859/ 340) -56.3 E30-E27 S37-S34
     18  651  224 (S42E54 L=275)  1064     85        1 ( 611/ 191- 689/ 259)  13.4 E68-E47 S45-S38
  

    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