Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1118 1840 (N53E26 L=155)  1665    179        2 (1033/1799-1188/1870)  21.0 E33-E18 N48-N56
      2  817 1790 (N38E43 L=138)   416     41        1 ( 794/1764- 839/1808) -18.4 E45-E39 N35-N39
      3 1306 1710 (N48E02 L=178)   315     38        0 (1286/1689-1332/1726)  53.1 E04-W00 N45-N49
      4  451 1630 (N17E60 L=121)  1105    129        1 ( 422/1591- 489/1691) -49.4 E64-E54 N14-N21
      5  492 1507 (N11E48 L=133)   387     66        0 ( 467/1476- 520/1534)  90.0 E49-E46 N08-N12
      6 1718 1414 (N40W41 L=222)   750    139        1 (1649/1397-1776/1428)  20.6 W33-W49 N36-N43
      7 1125 1404 (N20E04 L=176)   450    155        1 (1043/1369-1192/1424)  12.5 E09-W00 N18-N21
      8 1118 1240 (N10E00 L=180)   213     12        0 (1099/1223-1140/1258)  90.0 E01-W00 N08-N11
      9 1876 1106 (N25W61 L=241)   492     64        1 (1852/1081-1900/1131)  -8.1 W57-W64 N23-N26
     10 1738 1073 (N18W45 L=226)   591     59        1 (1711/1034-1758/1107)  90.0 W42-W46 N15-N20
     11  895  938 (S12E05 L=176)   530    167        0 ( 863/ 852- 921/1017) -84.3 E06-E04 S17-S07
     12  298  980 (S24E48 L=132)   376    139        1 ( 247/ 956- 352/1009)  15.3 E54-E43 S26-S21
     13  503  921 (S23E30 L=151)   290     55        1 ( 474/ 898- 521/ 947)  90.0 E32-E28 S24-S21
     14 1167  824 (S12W13 L=194)   266     45        0 (1137/ 805-1194/ 843)   0.0 W11-W15 S12-S10
     15  425  739 (S36E34 L=146)   557    107        1 ( 386/ 699- 462/ 781)  45.0 E38-E31 S39-S32
     16 1653  494 (S16W61 L=242)   466     57        0 (1624/ 468-1683/ 517)  63.4 W59-W63 S17-S13
     17 1485  432 (S24W51 L=232)   658    151        0 (1429/ 378-1536/ 482)  83.7 W49-W52 S29-S20
  

    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