Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1689 1135 (N20W39 L= 43)   684    155        0 (1637/1109-1752/1162)  11.3 W35-W45 N18-N21
      2  326 1052 (S19E46 L=318)   410    135        0 ( 280/1018- 378/1068)  32.0 E49-E41 S21-S16
      3 1532  978 (N07W31 L= 35)   346     76        0 (1498/ 930-1555/1021)  80.5 W29-W32 N03-N09
      4  888  767 (S22E01 L=  3)   488    150        1 ( 807/ 724- 942/ 811)  58.0 E05-W00 S26-S17
      5 1240  703 (S16W21 L= 25)   310    125        1 (1170/ 682-1299/ 727)  35.5 W17-W24 S18-S13
      6  842  501 (S40W04 L=  8)   972    199        1 ( 768/ 426- 906/ 590) -31.6 E04-W10 S43-S35
      7 1224  557 (S25W26 L= 30)   352    112        0 (1161/ 542-1278/ 573)  23.2 W21-W28 S26-S22
      8 1043  428 (S39W20 L= 24)   337     84        0 (1001/ 409-1085/ 449)  18.4 W17-W23 S39-S37
      9 1533  414 (S24W54 L= 58)   748    211        1 (1463/ 395-1614/ 426)  22.6 W49-W61 S26-S21
  

    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