Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1582 1665 (N24W64 L=275)   588     87        0 (1538/1639-1619/1686) -59.0 W61-W67 N21-N26
      2 1036 1540 (N35W17 L=228)  1860    548        4 ( 814/1498-1290/1586) -25.0 W01-W31 N25-N40
      3  434 1528 (N54E34 L=177)   324     77        0 ( 423/1492- 446/1565) -71.6 E36-E32 N50-N56
      4 1729 1468 (N09W62 L=273)   428     79        0 (1702/1444-1748/1507)   6.3 W58-W67 N07-N10
      5 1398 1331 (N12W31 L=242)  1103    396        1 (1219/1247-1524/1388) -40.9 W20-W36 N04-N18
      6  906 1253 (N22E01 L=210)   786    244        1 ( 810/1220-1030/1289) -13.0 E06-W07 N20-N24
      7  388 1183 (N31E39 L=172)   339     56        0 ( 372/1146- 407/1213) -78.7 E39-E36 N28-N33
      8 1887 1173 (S12W66 L=277)   891    127        0 (1849/1141-1926/1202) -22.6 W61-W73 S14-S09
      9  885 1125 (N15E05 L=206)   370    174        0 ( 843/1064- 962/1166)  14.0 E08-W00 N12-N16
     10  319  952 (N19E48 L=163)   742    221        2 ( 249/ 928- 412/1007) -29.7 E54-E40 N15-N23
     11 1670  862 (S22W36 L=247)   500    158        0 (1609/ 825-1726/ 919)   5.2 W30-W41 S23-S20
     12  222  741 (N08E65 L=146)  1481    127        3 ( 196/ 686- 257/ 803)  74.1 E68-E60 N04-N11
     13  304  709 (N05E57 L=154)   488     73        0 ( 277/ 674- 332/ 749) -68.2 E58-E54 N02-N07
     14 1779  581 (S46W49 L=260)  1128     91        2 (1729/ 554-1823/ 613) -10.0 W40-W57 S48-S43
     15  913  489 (S22E22 L=189)   739    218        1 ( 877/ 390- 934/ 583)  47.7 E28-E17 S27-S16
     16 1610  466 (S48W27 L=238)   515    104        0 (1560/ 442-1653/ 493)   0.0 W21-W32 S49-S46
     17 1250  252 (S51E13 L=198)  1033    199        1 (1170/ 220-1352/ 298) -12.5 E20-E02 S53-S48
  

    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