Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  686 1735 (N38E37 L=172)   373     31        0 ( 671/1706- 703/1760) -53.1 E39-E34 N36-N40
      2  504 1582 (N25E44 L=164)   294     38        1 ( 488/1552- 521/1609) -45.0 E46-E42 N23-N27
      3  693 1450 (N18E26 L=183)   890    179        0 ( 668/1372- 738/1530) -59.0 E28-E21 N13-N23
      4 1788 1240 (N17W52 L=261)   414     19        0 (1765/1222-1810/1261)  21.8 W49-W54 N15-N18
      5 1648  918 (S04W42 L=250)   249     35        0 (1632/ 893-1666/ 944) -71.6 W40-W43 S05-S02
      6  315  607 (S37E56 L=153)   698    159        0 ( 265/ 549- 367/ 669) -48.8 E59-E51 S40-S32
      7 1647  438 (S33W62 L=271)   697    124        0 (1592/ 409-1696/ 479)  40.6 W56-W65 S35-S29
      8  817  294 (S57E10 L=199)  2342    325        7 ( 671/ 263- 953/ 330)   5.2 E28-W05 S60-S53
      9 1029  295 (S54W12 L=221)   488     53        1 (1001/ 274-1059/ 323) -15.9 W08-W15 S56-S51
  

    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