Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1407 1820 (N48W69 L=326)   917     72        1 (1375/1792-1440/1847)   2.4 W59-W83 N46-N49
      2 1528 1672 (N35W60 L=317)  1321    247        0 (1479/1600-1569/1725)  60.9 W55-W67 N30-N39
      3  697 1513 (N43E16 L=241)   544    124        2 ( 645/1486- 772/1537) -20.0 E20-E09 N39-N45
      4 1817 1385 (N11W70 L=327)   663     78        0 (1798/1352-1839/1416) -63.4 W68-W72 N08-N12
      5 1249 1280 (N18W19 L=276)  1111    386        0 (1067/1225-1405/1326) -31.0 W07-W27 N11-N23
      6  473 1305 (N32E35 L=222)   261     58        0 ( 440/1289- 508/1324) -18.4 E38-E32 N30-N33
      7  211 1203 (N28E63 L=194)   849     81        1 ( 178/1172- 247/1229)  20.6 E68-E58 N25-N29
      8  745 1199 (N21E15 L=243)   215     73        1 ( 710/1180- 778/1229)  21.8 E17-E12 N20-N22
      9  349 1127 (N22E46 L=211)   298     82        0 ( 314/1114- 399/1147) -20.6 E49-E41 N20-N23
     10  204 1102 (N21E63 L=195)  1220    100        3 ( 169/1065- 250/1138) -35.5 E68-E56 N18-N23
     11 1707  999 (S08W45 L=302)   770    137        2 (1679/ 933-1733/1053)  74.1 W41-W47 S11-S04
     12  155  921 (N10E71 L=187)   763     70        1 ( 133/ 883- 173/ 953) -59.0 E75-E67 N07-N12
     13  842  511 (S23E21 L=236)   464    180        0 ( 799/ 447- 872/ 588)  48.4 E26-E18 S26-S17
     14 1688  528 (S40W46 L=303)   411     44        0 (1662/ 511-1714/ 541) -15.9 W41-W49 S41-S38
     15 1585  410 (S47W37 L=294)   700    103        1 (1544/ 382-1629/ 440)  -5.2 W31-W42 S49-S44
     16  588  329 (S33E50 L=207)   499     82        1 ( 562/ 297- 612/ 369)  28.6 E56-E45 S35-S29
     17 1285  253 (S54W05 L=263)   592     71        2 (1251/ 228-1320/ 276)   0.0 W01-W09 S56-S50
     18  957  203 (S51E28 L=229)   440     36        1 ( 937/ 181- 981/ 224)  29.7 E31-E24 S52-S48
  

    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