Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1531 1703 (N28W63 L=260)   598     77        0 (1503/1672-1558/1737)  59.0 W60-W66 N25-N30
      2 1175 1678 (N39W34 L=231)   554    110        1 (1125/1661-1233/1696) -20.6 W31-W39 N36-N41
      3 1388 1603 (N27W43 L=239)   853    194        1 (1294/1585-1477/1626) -38.7 W37-W50 N23-N31
      4  997 1604 (N40W18 L=214)   316     71        1 ( 957/1582-1030/1619) -21.8 W15-W20 N38-N41
      5 1460 1458 (N17W40 L=236)   754    310        1 (1310/1435-1591/1484) -22.8 W29-W48 N12-N21
      6 1657 1379 (N07W50 L=247)   315     64        0 (1634/1351-1674/1405) -76.0 W49-W51 N04-N08
      7 1078 1344 (N22W12 L=208)   584    154        1 (1004/1322-1148/1371) -11.3 W07-W17 N20-N24
      8  584 1212 (N28E22 L=174)   561    238        0 ( 490/1149- 670/1274) -45.0 E28-E17 N22-N33
      9 1865 1029 (S19W58 L=254)   429     64        1 (1842/ 991-1883/1063)  31.0 W54-W61 S20-S17
     10  471 1008 (N19E34 L=162)   672    235        2 ( 400/ 976- 591/1042) -23.2 E39-E25 N16-N22
     11  358  953 (N18E45 L=152)   333     47        0 ( 327/ 936- 383/ 969) -31.0 E47-E42 N16-N19
     12 1769  921 (S22W45 L=242)   330     74        0 (1727/ 909-1809/ 935) -14.0 W41-W49 S23-S19
     13  332  813 (N10E50 L=146)   294     41        0 ( 318/ 785- 346/ 841)  56.3 E51-E48 N08-N12
     14 1437  762 (S20W17 L=213)   546    252        0 (1368/ 677-1515/ 839)  24.8 W10-W23 S23-S17
     15  355  749 (N06E50 L=146)   574      0       -1 ( 309/ 730- 407/ 771) -20.6 E54-E46 N04-N08
     16 1140  653 (S18E03 L=193)   279     82        1 (1123/ 625-1169/ 688) -78.7 E04-E02 S20-S15
     17 1574  430 (S49W23 L=219)   383     74        0 (1542/ 404-1602/ 460)  18.4 W19-W26 S50-S47
     18  494  303 (S24E66 L=130)   758     61        2 ( 472/ 268- 516/ 338)  45.0 E70-E62 S27-S21
     19 1345  298 (S51E02 L=194)   783     91        0 (1302/ 272-1393/ 335)   0.0 E06-W03 S53-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