Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  768 1870 (N47E44 L=197)   788    162        1 ( 750/1835- 788/1907) -35.0 E48-E38 N43-N51
      2 1025 1539 (N25E09 L=231)   607    164        0 ( 982/1459-1076/1599)  79.7 E10-E06 N18-N29
      3 1696 1489 (N37W41 L=282)  1142    152        2 (1627/1455-1750/1521)   9.5 W34-W46 N34-N39
      4 1351 1434 (N24W13 L=254)   791    177        0 (1304/1348-1399/1519) -48.4 W09-W17 N19-N28
      5  294 1458 (N10E60 L=181)   882    130        4 ( 265/1413- 333/1499)  71.6 E62-E55 N06-N12
      6  360 1411 (N07E51 L=189)   708    140        2 ( 306/1346- 399/1477)  66.0 E54-E48 N02-N11
      7  420 1189 (S06E41 L=200)   595    117        0 ( 399/1128- 451/1255)  82.9 E42-E39 S09-S01
      8 1173 1212 (N07W05 L=246)   337     77        0 (1134/1168-1202/1248) -38.7 W02-W07 N04-N08
      9 1664 1200 (N16W38 L=279)   628    100        1 (1631/1136-1686/1243) -63.4 W34-W39 N12-N18
     10 1886 1056 (N13W60 L=301)   445     82        1 (1868/1032-1909/1084)  45.0 W58-W63 N11-N15
     11  197  856 (S27E58 L=183)   454    130        2 ( 177/ 820- 221/ 901)  45.0 E60-E55 S29-S24
     12 1304  804 (S14W21 L=262)  2399    597        4 (1007/ 740-1570/ 845)  10.7 W02-W39 S19-S09
     13 1615  789 (S09W42 L=283)   268     63        1 (1588/ 759-1640/ 822)  76.0 W40-W43 S10-S06
     14  672  654 (S34E17 L=223)  1159    409        3 ( 514/ 593- 846/ 694)  23.7 E30-E05 S41-S29
     15  464  571 (S43E35 L=206)   562     68        2 ( 430/ 542- 495/ 602)  -9.5 E38-E32 S45-S40
     16 1709  494 (S22W62 L=303)  1668    245        1 (1617/ 460-1778/ 538)  28.3 W53-W67 S25-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