Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  709 1789 (N39E51 L=138)  1421    170        3 ( 640/1730- 766/1842) -12.1 E68-E40 N35-N42
      2  740 1630 (N28E35 L=154)  1121    226        1 ( 629/1591- 825/1673)  45.0 E41-E31 N23-N33
      3 1585 1511 (N42W31 L=221)  1798    223        4 (1522/1412-1633/1598) -59.7 W25-W36 N35-N48
      4  425 1485 (N13E52 L=137)   747    121        0 ( 392/1434- 464/1532)  82.9 E54-E50 N08-N16
      5 1539 1300 (N26W27 L=217)  2096    646        2 (1375/1030-1650/1506) -40.4 W12-W39 N11-N35
      6 1129 1409 (N23E02 L=187)   448    159        1 (1057/1363-1184/1464) -18.4 E08-W01 N21-N25
      7  994 1317 (N14E08 L=181)   197     54        2 ( 976/1288-1015/1343) -56.3 E09-E06 N12-N15
      8 1286 1243 (N16W11 L=200)   301    139        1 (1224/1223-1349/1263)  23.2 W06-W14 N14-N18
      9 1825 1170 (N25W56 L=245)  1330    148        1 (1803/1087-1844/1236)  84.3 W53-W58 N19-N29
     10  279 1028 (S17E52 L=138)  2034    394        6 ( 176/ 873- 365/1137)  42.0 E65-E45 S28-S09
     11  373  973 (S19E42 L=147)   337     87        0 ( 354/ 936- 403/ 998)  63.4 E43-E39 S20-S16
     12  526  897 (S20E29 L=161)   384    118        0 ( 473/ 857- 564/ 940)  49.4 E32-E26 S23-S16
     13 1293  838 (S07W20 L=209)   267     94        0 (1246/ 820-1336/ 857)  31.0 W17-W22 S08-S05
     14 1508  791 (S05W35 L=225)   279     97        0 (1470/ 761-1558/ 811)  38.7 W33-W38 S07-S02
     15 1458  721 (S10W34 L=223)   306     63        1 (1435/ 688-1481/ 753) -45.0 W31-W35 S11-S08
     16 1445  651 (S14W35 L=225)   401    107        3 (1412/ 616-1475/ 688) -31.0 W32-W38 S16-S12
     17  368  613 (S43E44 L=146)  1185    155        2 ( 291/ 592- 442/ 644)   8.1 E55-E34 S45-S40
     18 1337  620 (S19W29 L=219)   341     77        0 (1290/ 608-1374/ 642)  38.7 W26-W31 S20-S16
     19 1527  509 (S21W48 L=238)   295     89        0 (1484/ 484-1565/ 528)  45.0 W45-W51 S23-S18
  

    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