Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1125 1794 (N51W42 L=288)   484     69        2 (1095/1762-1145/1832)  26.6 W37-W48 N48-N53
      2  666 1424 (N40E14 L=232)   858    200        2 ( 590/1380- 750/1454) -26.6 E20-E08 N34-N42
      3 1182 1417 (N25W20 L=267)   735    171        2 (1093/1395-1260/1445) -20.0 W14-W25 N23-N27
      4 1083 1315 (N22W11 L=258)   294     92        1 (1057/1273-1115/1356) -80.5 W10-W11 N19-N25
      5  980 1225 (N20W02 L=249)   244    126        0 ( 923/1206-1034/1256)   0.0 E01-W06 N18-N20
      6 1859 1194 (S07W64 L=310)   474     36        1 (1837/1166-1875/1216)  56.3 W60-W66 S08-S05
      7  734 1150 (N21E15 L=232)   739    221        1 ( 641/1080- 831/1203) -45.0 E19-E09 N15-N25
      8  300 1013 (N21E50 L=196)   548     98        0 ( 259/ 976- 334/1048) -40.6 E54-E47 N18-N24
      9 1893  921 (S24W63 L=309)  1993    185        2 (1845/ 853-1928/1017)  38.7 W55-W70 S28-S19
     10 1733  583 (S41W43 L=290)  3857    443        7 (1599/ 440-1862/ 723)  26.6 W28-W62 S49-S35
     11 1398  668 (S23W14 L=260)   306    130        0 (1343/ 631-1448/ 707)   7.1 W09-W17 S24-S21
     12  524  500 (S13E47 L=200)   651    110        0 ( 489/ 453- 569/ 550) -81.9 E47-E44 S16-S09
     13 1420  336 (S50W09 L=256)   907    145        2 (1362/ 294-1483/ 380)   8.7 W03-W16 S52-S46
     14 1079  282 (S42E20 L=227)   595     90        0 (1045/ 259-1117/ 308)   8.1 E22-E15 S43-S39
     15  604  277 (S30E57 L=190)   641     77        0 ( 562/ 265- 639/ 295)  -6.3 E62-E52 S31-S29
     16 1262  251 (S52E08 L=238)   523     62        2 (1234/ 224-1295/ 272)   8.1 E11-E04 S54-S49
  

    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