Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1115 1868 (N56W38 L=  1)  1000    108        3 (1082/1829-1155/1899) -80.5 W32-W46 N52-N59
      2 1031 1828 (N55W26 L=349)   870    111        2 ( 977/1784-1077/1874)  18.4 W18-W33 N51-N58
      3 1484 1774 (N38W61 L= 24)  2884    296        5 (1341/1701-1603/1827) -32.3 W49-W73 N31-N45
      4  888 1582 (N38W02 L=325)  1028    268        2 ( 747/1556- 998/1603) -12.5 E08-W10 N35-N41
      5 1237 1581 (N31W26 L=349)   228     66        0 (1203/1567-1269/1596) -26.6 W24-W28 N29-N32
      6  506 1360 (N31E29 L=294)  1798    497        4 ( 294/1297- 756/1401) -16.7 E50-E10 N24-N37
      7 1199 1353 (N17W17 L=340)   731    230        2 (1137/1308-1277/1395) -14.0 W12-W21 N14-N20
      8  281 1243 (N27E50 L=273)   852    189        3 ( 241/1147- 314/1323) -58.4 E55-E45 N20-N33
      9  983 1270 (N16W02 L=325)   527    225        0 ( 882/1244-1085/1291) -15.3 E03-W08 N13-N19
     10  165 1244 (N29E66 L=257)   636     81        0 ( 137/1212- 197/1287) -24.8 E73-E60 N26-N32
     11  585 1229 (N21E24 L=299)   434    179        0 ( 500/1208- 670/1252) -24.4 E29-E18 N18-N23
     12  628 1065 (N10E23 L=300)   387    122        2 ( 581/1019- 668/1111) -56.3 E25-E20 N07-N13
     13 1811  974 (S16W53 L= 16)   606    137        1 (1785/ 913-1825/1046)  63.4 W50-W55 S19-S11
     14  359  823 (N01E47 L=276)   380    112        0 ( 318/ 784- 398/ 867)  14.0 E51-E42 N00-N02
     15  461  673 (S09E43 L=280)   354     51        1 ( 435/ 645- 494/ 696)  21.8 E44-E39 S10-S08
     16  586  609 (S15E36 L=287)   394    150        0 ( 534/ 576- 653/ 656) -54.5 E38-E31 S18-S11
     17  710  415 (S30E35 L=288)   736    199        0 ( 693/ 320- 741/ 494)  53.1 E40-E29 S36-S24
     18  868  449 (S31E22 L=301)   287    108        0 ( 830/ 410- 913/ 480)  20.6 E26-E18 S32-S29
     19  594  270 (S38E57 L=266)   651     61        2 ( 563/ 242- 622/ 303)  14.0 E63-E51 S39-S35
  

    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