Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  455 1646 (N17E59 L=117)  1208    198        3 ( 428/1556- 487/1721) -59.7 E62-E53 N10-N22
      2  579 1463 (N08E39 L=137)   578     80        1 ( 535/1429- 620/1498)  56.3 E41-E37 N05-N11
      3 1604 1251 (N25W29 L=206)   340    109        0 (1590/1202-1618/1304)  90.0 W28-W30 N21-N28
      4  216 1037 (S24E56 L=120)  1097    271        5 ( 163/ 954- 265/1139)  80.5 E60-E52 S29-S17
      5 1778  871 (N09W53 L=229)  1108    302        2 (1724/ 739-1838/1012) -86.8 W49-W55 N00-N18
      6  404  911 (S27E37 L=140)   907    259        4 ( 319/ 871- 500/ 939)  29.7 E43-E29 S31-S23
      7  946  881 (S15E01 L=176)   184    112        0 ( 884/ 866- 995/ 891)  31.0 E03-W02 S17-S13
      8  978  772 (S20W04 L=181)   210    134        0 ( 907/ 762-1038/ 786)  29.7 W00-W07 S21-S17
      9 1833  747 (N04W64 L=241)   450     67        1 (1813/ 723-1852/ 771)  71.6 W62-W67 N02-N05
     10  671  581 (S41E10 L=166)   233     72        1 ( 655/ 553- 691/ 618) -76.0 E11-E08 S42-S38
     11 1620  443 (S19W60 L=236)  2460    426       10 (1483/ 353-1739/ 563)  56.3 W53-W67 S27-S09
     12 1431  413 (S27W47 L=224)   342    123        0 (1419/ 366-1447/ 461) -50.2 W44-W51 S29-S23
     13 1387  360 (S31W48 L=224)   410     70        2 (1364/ 334-1409/ 397) -59.0 W45-W49 S33-S28
  

    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