Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  503 1713 (N27E57 L=291)  2324    244        7 ( 439/1622- 568/1812) -74.1 E65-E52 N20-N36
      2  445 1710 (N27E64 L=284)   706     93        0 ( 419/1671- 473/1752) -66.8 E66-E61 N23-N30
      3 1336 1667 (N41W08 L=355)  1013    249        4 (1237/1611-1443/1706)  -7.1 W00-W16 N38-N43
      4  472 1559 (N16E49 L=299)   840    205        0 ( 417/1486- 536/1635)  80.5 E50-E47 N10-N22
      5  580 1512 (N15E38 L=309)   278     75        1 ( 564/1482- 604/1543)  76.0 E39-E36 N12-N16
      6  715 1365 (N08E25 L=322)   206     56        1 ( 693/1341- 736/1390)  76.0 E25-E23 N05-N09
      7 1362 1297 (N16W15 L=  2)   396     93        0 (1313/1263-1403/1318)  -9.5 W11-W17 N14-N16
      8 1724 1268 (N23W42 L= 29)   282     54        0 (1701/1241-1745/1294) -45.0 W39-W44 N21-N24
      9 1184 1254 (N09W05 L=352)   264     51        1 (1164/1226-1209/1282) -71.6 W03-W06 N07-N11
     10  326 1112 (S13E48 L=300)   613    178        1 ( 264/1041- 387/1155)  41.6 E52-E42 S17-S09
     11 1416 1117 (N06W21 L=  9)   208     49        1 (1392/1092-1431/1149)  90.0 W19-W22 N04-N08
     12 1640 1075 (N09W37 L= 24)   412    106        0 (1626/1024-1672/1117)  80.5 W35-W38 N06-N12
     13  497  998 (S17E33 L=315)   254     98        0 ( 486/ 943- 509/1043)  90.0 E33-E31 S20-S14
     14  381  749 (S34E42 L=306)   853    242        1 ( 282/ 702- 463/ 825)  25.3 E53-E33 S38-S28
     15 1769  722 (S07W58 L= 45)   399     56        1 (1751/ 690-1794/ 746)  90.0 W55-W60 S09-S05
     16 1755  595 (S14W63 L= 51)  2137    252       11 (1663/ 543-1832/ 660)  28.1 W55-W73 S19-S09
     17 1471  542 (S25W41 L= 29)   950    228        2 (1428/ 444-1519/ 654) -70.3 W38-W44 S31-S17
     18 1777  535 (S16W71 L= 59)   828     71        2 (1756/ 502-1798/ 565) -76.0 W67-W75 S18-S14
     19 1564  371 (S32W62 L= 49)   561     59        1 (1536/ 346-1599/ 389)  56.3 W57-W65 S33-S29
  

    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