Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  501 1714 (N70E37 L=102)  1277    156        2 ( 434/1677- 572/1744)  -2.9 E62-E21 N66-N72
      2 1103 1645 (N41W27 L=166)   376     80        1 (1066/1622-1144/1659)  -8.1 W23-W30 N39-N42
      3 1525 1518 (N20W48 L=187)   758    196        2 (1450/1463-1591/1583) -73.3 W46-W52 N15-N25
      4 1656 1504 (N15W58 L=197)   741     79        3 (1631/1460-1679/1540)  29.7 W54-W62 N12-N16
      5 1256 1376 (N20W24 L=163)   577    333        1 (1109/1321-1387/1410) -38.2 W14-W29 N13-N25
      6  719 1270 (N29E12 L=127)   222    115        1 ( 682/1244- 765/1312) -45.0 E14-E09 N26-N32
      7  240 1283 (N40E58 L= 81)   512     42        1 ( 213/1262- 266/1302) -18.4 E63-E53 N38-N41
      8  593 1186 (N27E23 L=116)   356    106        2 ( 559/1152- 638/1226) -51.3 E25-E20 N24-N29
      9  231 1175 (N33E58 L= 82)   376     56        0 ( 219/1147- 245/1206) -76.0 E59-E55 N31-N35
     10  219  957 (N20E60 L= 80)   706    120        5 ( 166/ 924- 252/ 999)  18.4 E67-E55 N18-N22
     11  418  880 (N12E42 L= 97)   849    264        3 ( 315/ 854- 544/ 919)  -6.0 E50-E31 N10-N14
     12  238  880 (N15E59 L= 80)   375     49        2 ( 221/ 853- 257/ 909)  76.0 E61-E56 N13-N17
     13  572  834 (N07E32 L=108)   403    123        5 ( 543/ 770- 594/ 873) -90.0 E33-E29 N02-N09
     14  876  808 (S01E14 L=125)   301     73        0 ( 844/ 774- 908/ 832) -63.4 E15-E12 S04-N00
     15  551  526 (S11E44 L= 96)   430     74        1 ( 534/ 487- 568/ 566)  59.0 E45-E41 S13-S08
     16 1358  503 (S33W07 L=146)   592     73        3 (1325/ 472-1392/ 538) -90.0 W04-W09 S35-S30
     17  385  463 (S12E62 L= 77)  1989    201        3 ( 319/ 379- 437/ 518)  84.3 E68-E57 S18-S08
  

    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