Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1655 1548 (N20W57 L=156)   524     49        1 (1632/1520-1680/1575) -90.0 W55-W59 N18-N22
      2  250 1458 (N42E61 L= 38)   474     47        0 ( 231/1438- 268/1480) -29.7 E64-E56 N40-N44
      3  589 1390 (N31E22 L= 77)   583    220        0 ( 530/1308- 666/1432)  29.1 E27-E16 N26-N32
      4 1832 1040 (S12W56 L=155)   349     74        0 (1802/1010-1855/1061) -63.4 W53-W59 S14-S10
      5 1560 1013 (S09W33 L=132)   381     86        1 (1536/ 960-1581/1058)  68.2 W31-W34 S11-S06
      6  327  871 (N04E48 L= 51)   952    170        0 ( 256/ 833- 382/ 937)  17.1 E55-E42 N02-N07
      7 1110  851 (S10W02 L=101)   793    184        2 (1059/ 771-1166/ 925)  45.0 E02-W06 S14-S06
      8  519  761 (S05E36 L= 63)   486     85        1 ( 492/ 717- 558/ 793)  18.4 E38-E32 S07-S03
      9 1835  711 (S33W59 L=158)   875    110        0 (1803/ 660-1860/ 753) -41.6 W53-W65 S37-S29
     10 1729  672 (S33W46 L=145)   818    154        2 (1681/ 620-1778/ 748)  36.9 W41-W50 S36-S28
     11  527  611 (S14E40 L= 59)   520     48        1 ( 501/ 588- 561/ 643)  11.3 E42-E37 S15-S12
  

    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