Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  994 1845 (N59W38 L=290)   618     77        2 ( 964/1808-1023/1877)  12.1 W31-W45 N56-N61
      2  859 1779 (N60W18 L=270)   308     79        1 ( 840/1747- 896/1798)  12.5 W14-W23 N57-N61
      3 1385 1726 (N36W53 L=305)   620     66        2 (1341/1705-1419/1752)  -6.3 W48-W57 N34-N37
      4 1644 1374 (N09W49 L=301)   341     62        1 (1607/1362-1683/1386) -26.6 W46-W52 N07-N10
      5 1800 1224 (S04W57 L=309)   931    128        0 (1775/1156-1826/1285)  32.0 W52-W61 S06-N00
      6  863 1240 (N23E04 L=248)   296     91        1 ( 820/1213- 901/1276) -51.3 E06-E02 N20-N25
      7  839 1152 (N19E08 L=245)   586    244        0 ( 739/1116- 948/1181) -13.0 E14-E01 N15-N20
      8  869  982 (N08E10 L=243)   223    160        0 ( 838/ 913- 905/1049)  40.6 E13-E06 N05-N11
      9 1403 1002 (S04W21 L=273)   285     67        0 (1377/ 965-1421/1028)  56.3 W19-W22 S05-S02
     10  278  875 (N14E54 L=199)   740    158        1 ( 243/ 794- 310/ 938) -73.3 E56-E51 N08-N18
     11  169  854 (N14E68 L=185)   458     66        0 ( 153/ 820- 188/ 896) -78.7 E69-E65 N11-N16
     12 1743  803 (S26W42 L=294)   607     92        1 (1701/ 760-1780/ 838)  23.2 W37-W45 S27-S23
     13  153  771 (N09E74 L=179)   614     74        1 ( 136/ 736- 171/ 798) -45.0 E76-E70 N06-N10
     14  319  505 (S09E65 L=187)   778     94        3 ( 292/ 457- 340/ 549)  60.3 E68-E62 S12-S05
     15  464  426 (S16E56 L=197)  3039    305        8 ( 383/ 334- 553/ 528) -48.0 E65-E49 S22-S09
     16  982  467 (S24E17 L=235)   654    148        1 ( 924/ 419-1041/ 521) -73.3 E19-E15 S29-S19
     17  638  430 (S19E41 L=211)  1130    220        3 ( 579/ 349- 689/ 504)  39.3 E49-E35 S24-S14
     18 1260  344 (S42E04 L=248)  2685    597        7 (1045/ 208-1497/ 503)  10.5 E26-W17 S47-S36
  

    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