Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1459 1531 (N32W25 L=233)   279     91        0 (1440/1485-1473/1567) -80.5 W23-W26 N28-N34
      2  724 1431 (N17E23 L=184)  1123    268        1 ( 678/1327- 813/1527) -47.5 E27-E15 N10-N22
      3 1802 1243 (N18W53 L=261)   400     23        0 (1777/1224-1823/1264)  21.8 W50-W55 N16-N19
      4 1662  917 (S04W42 L=250)   279     37        1 (1643/ 894-1684/ 945) -56.3 W40-W44 S05-S02
      5  628  748 (S25E24 L=184)   239    107        0 ( 584/ 708- 671/ 774)  39.8 E26-E20 S28-S23
      6  317  603 (S37E55 L=153)   771    149        0 ( 269/ 543- 369/ 661) -48.8 E59-E50 S40-S32
      7 1654  437 (S33W62 L=270)   564    122        0 (1598/ 410-1702/ 474)  35.5 W56-W65 S34-S29
      8  857  289 (S57E05 L=202)  2769    458        5 ( 672/ 258-1067/ 322)   9.0 E29-W16 S61-S52
  

    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