Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1165 1825 (N51E07 L=101)  1876    191        5 (1126/1727-1208/1884)  86.8 E11-E03 N41-N59
      2  984 1826 (N47E22 L= 86)   407     88        0 ( 964/1782-1003/1868) -45.0 E26-E18 N43-N51
      3  790 1594 (N25E25 L= 82)   848    232        3 ( 702/1513- 866/1660) -28.1 E33-E18 N20-N28
      4 1807 1372 (N31W54 L=162)   591    119        0 (1770/1326-1852/1405) -16.7 W49-W59 N28-N32
      5  232 1149 (S08E57 L= 50)  4182    350       12 ( 137/1036- 321/1308)  43.2 E69-E47 S15-N01
      6 1367 1230 (N11W17 L=125)   216     54        2 (1336/1208-1394/1251)  45.0 W15-W18 N09-N13
      7 1819 1204 (N19W53 L=161)   335     49        1 (1799/1180-1836/1227) -90.0 W50-W55 N18-N21
      8 1387 1161 (N08W19 L=127)   207     81        0 (1371/1121-1402/1207) -78.7 W18-W20 N05-N10
      9 1678 1074 (N08W41 L=148)   605    191        2 (1641/ 989-1704/1155) -60.9 W37-W43 N03-N12
     10  626 1047 (S11E24 L= 83)   191     92        0 ( 601/1011- 660/1087) -38.7 E26-E21 S12-S08
     11  185 1046 (S15E62 L= 46)   456     57        2 ( 166/1017- 207/1074)  31.0 E64-E59 S16-S13
     12  255  889 (S24E54 L= 54)  1675    328        7 ( 209/ 741- 293/1015)  54.8 E63-E49 S33-S16
     13  390  803 (S28E41 L= 66)  1524    443        3 ( 368/ 606- 441/1002)  72.9 E46-E37 S41-S15
     14 1029  969 (S09W01 L=109)   205     76        0 ( 987/ 954-1066/ 987)  26.6 E01-W03 S10-S07
     15 1756  571 (S18W65 L=172)  1170    132        0 (1710/ 526-1806/ 630)  58.0 W61-W68 S21-S13
     16 1200  458 (S37W24 L=132)   285     77        0 (1183/ 417-1217/ 490) -71.6 W22-W26 S40-S34
     17 1584  380 (S33W62 L=169)  1583    189        3 (1501/ 323-1666/ 417)  42.0 W56-W70 S38-S29
  

    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