Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1012 1729 (N36E19 L= 88)   686    133        2 ( 943/1698-1073/1757)  -5.2 E24-E13 N34-N37
      2  381 1560 (N13E57 L= 50)   736     96        4 ( 349/1518- 418/1598)  68.2 E60-E54 N10-N16
      3 1209 1512 (N26W00 L=108)   339    111        0 (1175/1470-1252/1556)  74.1 E01-W02 N22-N29
      4 1553 1514 (N35W25 L=132)   877    173        2 (1472/1467-1627/1546)   4.4 W18-W31 N32-N37
      5 1047 1267 (N07E04 L=104)   258     72        0 (1011/1247-1088/1287)  36.9 E05-E01 N05-N08
      6  171 1182 (S12E65 L= 42)  1333    253        1 ( 117/1082- 231/1273)  45.0 E71-E59 S17-S05
      7  286 1148 (S13E51 L= 56)   328     60        1 ( 251/1130- 316/1169)   9.5 E54-E48 S13-S11
      8 1812 1102 (N16W50 L=158)   803    114        0 (1765/1052-1849/1147) -18.4 W45-W55 N14-N17
      9 1417 1047 (N03W22 L=129)   327     95        0 (1403/ 995-1441/1098) -51.3 W20-W24 N00-N05
     10  882  981 (S12E07 L=100)  1062    528        1 ( 670/ 902-1080/1060)  31.8 E18-W03 S20-S05
     11  290  962 (S23E49 L= 58)   413     70        1 ( 255/ 943- 321/ 987)  26.6 E52-E46 S24-S21
     12 1841  862 (N04W59 L=166)  1270    237        1 (1806/ 752-1860/ 977) -76.0 W55-W61 S02-N10
     13 1278  736 (S18W21 L=129)   246     25        1 (1259/ 714-1305/ 758)  45.0 W19-W22 S19-S16
     14 1119  691 (S24W13 L=121)   266     86        0 (1081/ 675-1162/ 715)   0.0 W10-W16 S24-S22
     15 1787  596 (S11W66 L=173)  1434    195        6 (1727/ 536-1845/ 657) -21.0 W57-W75 S14-S08
     16 1148  611 (S28W17 L=125)   218     40        1 (1122/ 592-1178/ 627)  56.3 W16-W19 S29-S26
     17 1732  561 (S15W62 L=169)   566    112        2 (1709/ 509-1758/ 609) -63.4 W60-W65 S17-S11
  

    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