Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  987 1922 (N64W48 L= 37)   501     42        0 ( 966/1896-1002/1942)  15.9 W41-W55 N61-N66
      2  836 1853 (N64W21 L= 10)  1399    186        3 ( 755/1814- 908/1903)   3.9 W08-W37 N60-N67
      3  717 1780 (N62W02 L=350)   472     27        0 ( 699/1751- 737/1806)  59.0 E01-W04 N59-N65
      4  942 1628 (N42W13 L=  2)  1439    332        0 ( 800/1597-1102/1661) -23.2 W01-W22 N35-N45
      5 1692 1559 (N16W63 L= 52)   680     59        0 (1655/1540-1728/1576) -12.5 W59-W68 N14-N17
      6 1241 1517 (N26W27 L= 16)   390     65        0 (1224/1475-1262/1551)  31.0 W24-W29 N23-N28
      7  403 1426 (N43E36 L=313)   432     43        0 ( 376/1404- 438/1444)   0.0 E39-E32 N41-N44
      8 1537 1314 (N07W38 L= 27)   703    181        0 (1498/1233-1571/1397)  84.3 W36-W39 N02-N12
      9 1827 1271 (S03W60 L= 49)   716    114        1 (1791/1224-1856/1329)  20.0 W55-W66 S04-N00
     10 1210 1111 (N04W12 L=  1)   373    138        0 (1170/1057-1254/1166) -76.0 W11-W13 N00-N08
     11  161 1025 (N22E63 L=286)  6232    255        5 (  94/ 923- 239/1133)  76.0 E76-E53 N16-N28
     12  390  967 (N15E41 L=308)  1362    352        6 ( 268/ 880- 511/1113) -48.4 E49-E33 N07-N25
     13 1208  810 (S12W05 L=354)   529    164        0 (1182/ 739-1238/ 896)  53.1 W02-W08 S15-S07
     14  154  742 (N06E73 L=276)  3294    169        5 ( 120/ 655- 196/ 810)  78.7 E84-E66 N00-N10
     15 1761  604 (S42W46 L= 34)  1076    172        1 (1683/ 582-1841/ 626) -16.9 W36-W59 S45-S37
     16 1593  520 (S42W27 L= 16)   484     80        0 (1559/ 491-1630/ 552)  15.9 W23-W30 S43-S40
  

    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