Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1035 1820 (N53E25 L=347)   410     45        1 (1006/1799-1058/1843)  36.9 E27-E22 N51-N55
      2 1428 1722 (N56W19 L= 31)  1282    170        3 (1339/1694-1495/1747)  11.3 W08-W30 N52-N61
      3  898 1372 (N17E15 L=357)   291    104        1 ( 859/1339- 930/1394)  -9.5 E17-E11 N15-N18
      4 1482 1210 (N19W25 L= 37)   225     64        0 (1462/1176-1500/1240) -53.1 W23-W26 N16-N20
      5  260 1053 (S14E54 L=317)   490    115        0 ( 244/ 999- 279/1104) -71.6 E55-E51 S16-S10
      6  706  961 (S11E18 L=354)   465    189        1 ( 628/ 903- 779/ 997)  45.0 E22-E13 S15-S07
      7  383  740 (S31E41 L=331)  2306    441        4 ( 249/ 619- 531/ 874) -18.4 E54-E27 S37-S23
      8 1415  358 (S34W48 L= 60)   372     50        1 (1390/ 338-1438/ 391) -26.6 W45-W51 S35-S32
  

    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