Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  585 1800 (N37E67 L=149)   651     61        1 ( 563/1772- 614/1828) -11.3 E78-E58 N34-N38
      2 1280 1600 (N40W03 L=218)  1050    224        5 (1211/1513-1322/1673) -42.0 E04-W07 N33-N44
      3  762 1553 (N23E30 L=185)   379    168        0 ( 709/1501- 800/1623) -22.6 E37-E25 N20-N26
      4 1244 1327 (N20W06 L=222)   857    463        2 (1179/1153-1325/1543) -54.5 E02-W14 N11-N32
      5  704 1441 (N15E30 L=186)   257     93        2 ( 676/1404- 726/1489) -38.7 E32-E27 N12-N17
      6  856 1385 (N15E18 L=197)   566    260        1 ( 776/1354- 971/1425)  24.4 E23-E11 N12-N18
      7 1110 1359 (N19E03 L=213)   217    129        1 (1091/1300-1125/1416) -81.9 E04-E01 N15-N22
      8 1556 1271 (N25W29 L=244)   718    207        2 (1528/1202-1608/1323) -54.5 W26-W32 N20-N28
      9  240 1120 (S13E58 L=158)  2319    420        5 ( 206/ 948- 285/1263) -84.0 E63-E53 S23-S04
     10 1800 1063 (N18W53 L=269)   910    163        4 (1784/ 983-1816/1143)  90.0 W51-W55 N13-N23
     11  790  948 (S12E12 L=203)   289    111        1 ( 742/ 924- 830/ 978)   9.5 E15-E09 S13-S10
     12 1227  884 (S06W15 L=230)   345    120        0 (1191/ 837-1254/ 934) -39.8 W11-W17 S08-S03
     13  959  759 (S19W02 L=218)   401    190        1 ( 878/ 737-1029/ 777)  36.9 E02-W06 S22-S16
     14  335  696 (S38E46 L=170)  3031    371        6 ( 171/ 654- 478/ 745)   7.7 E82-E30 S43-S32
     15  694  397 (S51E08 L=207)   370     74        1 ( 659/ 377- 729/ 425)   0.0 E12-E04 S52-S48
     16 1633  392 (S24W69 L=285)   747     68        3 (1611/ 363-1663/ 424)  -3.8 W63-W78 S26-S22
     17 1509  328 (S32W62 L=278)   720     51        2 (1482/ 305-1535/ 355) -11.3 W57-W67 S33-S30
  

    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