Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1034 1610 (N43W12 L=345)   266     49        2 (1009/1584-1060/1632)  36.9 W09-W14 N41-N44
      2 1434 1245 (N13W29 L=  3)   659    165        0 (1408/1173-1479/1316) -83.7 W27-W31 N08-N17
      3  299 1165 (N22E52 L=282)   302     42        1 ( 277/1136- 319/1187) -45.0 E54-E49 N20-N23
      4  724 1017 (N09E19 L=315)   206    101        0 ( 693/ 993- 765/1046) -21.8 E20-E15 N07-N11
      5 1660  918 (S10W40 L= 14)   781    241        1 (1618/ 872-1704/ 989)  74.1 W36-W44 S13-S06
      6 1625  726 (S22W37 L= 10)   504    130        0 (1585/ 663-1657/ 776)  45.0 W33-W39 S25-S19
      7 1435  691 (S21W22 L=355)   294    116        0 (1408/ 634-1456/ 743)  56.3 W19-W23 S24-S17
      8 1102  479 (S30E03 L=330)   385     95        0 (1081/ 428-1133/ 517)  56.3 E05-W00 S33-S27
      9  411  417 (S26E61 L=273)  1726    157        3 ( 382/ 339- 440/ 496)  29.5 E75-E52 S33-S20
  

    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