Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  614 1616 (N34E32 L=262)  1986    378        4 ( 487/1471- 731/1708) -29.5 E50-E20 N23-N42
      2  322 1432 (N22E53 L=241)   506     78        0 ( 303/1392- 353/1462)  63.4 E55-E50 N19-N24
      3 1805 1309 (N16W57 L=352)  1274    254        4 (1765/1209-1833/1411)  65.2 W52-W62 N10-N23
      4  381 1130 (N03E42 L=252)   550     90        3 ( 342/1101- 422/1163)  -8.1 E45-E38 N01-N04
      5 1765 1118 (N04W51 L=345)   533     50        2 (1745/1080-1786/1152)  78.7 W48-W52 N01-N06
      6 1688 1070 (N01W44 L=338)   471     83        3 (1660/1036-1729/1095)  31.0 W41-W47 S01-N02
      7  185  750 (S19E67 L=227)  4721    186        8 ( 135/ 667- 238/ 840)  70.0 E81-E60 S24-S13
      8  454  577 (S31E43 L=251)   927    164        1 ( 409/ 522- 486/ 643)  35.0 E49-E38 S34-S26
      9  296  553 (S32E63 L=232)   855    105        0 ( 255/ 509- 337/ 584)  11.9 E74-E55 S34-S30
     10 1153  225 (S60W17 L=312)  2820    414        6 ( 959/ 195-1334/ 267)   6.2 E06-W41 S63-S55
  

    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