Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1590 1597 (N31W48 L=152)  1420    228        1 (1495/1535-1668/1671) -45.0 W43-W53 N26-N37
      2 1393 1558 (N30W30 L=133)  1746    351        3 (1250/1458-1519/1668) -46.8 W21-W37 N22-N38
      3 1015 1621 (N36W03 L=107)   477     56        2 ( 983/1596-1053/1644) -26.6 W00-W05 N34-N37
      4  688 1485 (N29E20 L= 83)   453    129        1 ( 628/1449- 750/1519)  20.6 E24-E16 N26-N30
      5 1003 1400 (N21W01 L=104)   197     96        0 ( 986/1355-1025/1427) -78.7 W00-W02 N17-N22
      6  267 1357 (N23E55 L= 48)   927    124        2 ( 227/1315- 310/1405)  35.5 E59-E51 N20-N26
      7 1072 1341 (N17W05 L=108)   340    134        1 (1008/1322-1130/1369)   0.0 W00-W08 N15-N18
      8 1032  868 (S11E00 L=103)   379    113        0 ( 989/ 831-1071/ 918)  45.0 E03-W02 S13-S08
      9 1539  736 (S22W33 L=137)   896    273        0 (1495/ 610-1575/ 864)  82.4 W30-W35 S30-S14
     10 1040  670 (S24E01 L=103)   470    148        0 (1025/ 600-1061/ 737) -90.0 E02-W00 S28-S19
     11 1710  649 (S28W50 L=154)   412     39        2 (1687/ 622-1733/ 676) -45.0 W47-W52 S29-S26
     12 1776  577 (S33W62 L=166)  1854    198        2 (1696/ 497-1845/ 651)  34.7 W54-W69 S38-S28
     13 1536  480 (S39W38 L=142)   911    127        3 (1508/ 422-1559/ 537)  90.0 W34-W43 S43-S35
     14 1653  467 (S41W53 L=156)   528     34        1 (1628/ 443-1683/ 487) -14.0 W48-W57 S42-S39
     15  619  282 (S49E48 L= 56)   738     94        1 ( 590/ 258- 649/ 314)  15.3 E53-E42 S50-S45
     16 1377  248 (S59W35 L=138)   668     77        0 (1342/ 228-1413/ 269)  -9.5 W28-W40 S60-S56
  

    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