Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  757 1700 (N34E37 L=185)  1804    447        4 ( 633/1626- 886/1784) -10.5 E52-E25 N30-N37
      2  968 1605 (N32E18 L=204)   540    100        1 ( 923/1570-1029/1624)  39.8 E20-E14 N29-N34
      3  278 1357 (N03E60 L=162)   427     86        0 ( 266/1320- 296/1394) -38.7 E63-E57 N00-N05
      4 1276 1277 (N18W10 L=232)   244     82        0 (1245/1236-1311/1314) -31.0 W07-W12 N16-N19
      5 1628 1021 (N11W38 L=260)  1046    252        2 (1574/ 895-1659/1129) -61.7 W32-W41 N02-N16
      6  535 1055 (S10E31 L=191)   323    103        0 ( 483/1039- 579/1070)   8.1 E34-E27 S11-S08
      7 1360  801 (S08W25 L=247)   327    139        0 (1343/ 730-1382/ 871) -58.0 W22-W28 S11-S03
      8  583  716 (S30E22 L=200)   243     48        1 ( 556/ 700- 609/ 737)  14.0 E24-E20 S30-S28
  

    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