Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  438 1707 (N20E65 L=131)  2098    320        4 ( 410/1651- 489/1758) -45.0 E71-E61 N15-N24
      2  654 1461 (N09E34 L=162)   230     38        0 ( 637/1436- 671/1488)  90.0 E35-E33 N06-N10
      3  685 1020 (S15E19 L=177)   652    131        1 ( 620/ 978- 745/1067)  53.1 E22-E16 S19-S11
      4 1827  836 (N09W59 L=255)  1481    222        3 (1791/ 756-1871/ 905) -26.6 W53-W68 N05-N12
      5  810  839 (S22E07 L=190)   310     81        0 ( 778/ 812- 851/ 877) -11.3 E09-E04 S23-S20
      6 1405  765 (S09W29 L=226)   715    162        1 (1384/ 682-1425/ 842) -48.8 W25-W32 S13-S05
      7  989  764 (S21W06 L=202)   331    131        1 ( 927/ 744-1060/ 790)  20.6 W01-W09 S22-S19
      8 1347  710 (S14W28 L=224)   303    122        2 (1278/ 675-1384/ 738)   7.1 W22-W30 S14-S12
      9  579  643 (S41E18 L=178)   834    211        2 ( 481/ 624- 679/ 671)  14.0 E26-E10 S43-S37
     10 1416  461 (S25W44 L=240)  2543    899        1 (1096/ 407-1716/ 538)  15.5 W20-W74 S32-S15
     11 1217  258 (S43W46 L=242)   421     72        0 (1184/ 237-1243/ 275)   0.0 W42-W48 S44-S42
  

    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