Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1290 1806 (N51W52 L=170)  1559    219        2 (1188/1761-1383/1849)  -1.8 W36-W68 N49-N52
      2 1537 1695 (N36W63 L=181)   434     67        0 (1516/1674-1572/1711) -21.8 W58-W66 N34-N37
      3 1633 1324 (N12W47 L=165)   920    282        2 (1524/1286-1749/1377) -24.0 W39-W57 N08-N17
      4  324 1318 (N35E50 L= 68)   438    100        0 ( 299/1272- 354/1361)  50.2 E52-E47 N31-N37
      5  928 1240 (N21E02 L=116)   229    120        0 ( 901/1182- 966/1282)  51.3 E04-W00 N17-N23
      6 1099 1218 (N17W08 L=126)   388    175        0 (1021/1175-1169/1249) -32.0 W03-W11 N13-N19
      7  655 1086 (N16E22 L= 96)   626    208        1 ( 607/1010- 722/1168) -61.4 E24-E18 N10-N21
      8 1800 1009 (S10W53 L=172)  3231    287        9 (1692/ 929-1878/1102)   8.1 W42-W64 S15-S05
      9 1727 1044 (S06W47 L=165)   322     44        0 (1710/1015-1744/1070) -90.0 W45-W48 S07-S04
     10 1188  750 (S13W05 L=123)   536    113        2 (1149/ 699-1218/ 813)  66.8 W02-W06 S15-S08
     11  850  656 (S13E17 L=101)  1368    326        6 ( 753/ 574- 934/ 770) -65.0 E20-E13 S19-S04
     12 1686  714 (S26W41 L=159)   326    151        0 (1633/ 657-1720/ 765)  35.5 W36-W43 S28-S23
     13  384  535 (S15E56 L= 63)  1657    410        1 ( 326/ 399- 472/ 702)  84.8 E60-E50 S25-S03
     14  934  523 (S23E15 L=103)   364     99        1 ( 913/ 476- 956/ 569)  56.3 E17-E12 S25-S19
     15 1060  432 (S32E09 L=109)   627     98        2 (1014/ 406-1108/ 465)   0.0 E12-E05 S32-S29
     16  643  405 (S27E40 L= 79)   802     71        2 ( 618/ 359- 680/ 443)  54.5 E42-E36 S30-S23
     17 1307  272 (S52W07 L=125)   847    114        2 (1260/ 243-1351/ 303)  11.3 W01-W11 S55-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