Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  555 1471 (N17E40 L= 19)   257    105        0 ( 531/1415- 578/1516) -81.9 E40-E39 N12-N19
      2  775 1457 (N20E25 L= 35)   221     39        1 ( 754/1435- 797/1479)  71.6 E25-E23 N18-N21
      3 1149  818 (S11W11 L= 71)   358    102        0 (1132/ 768-1165/ 871) -59.0 W09-W12 S14-S08
      4  405  733 (S31E38 L= 22)   270     66        0 ( 384/ 718- 435/ 757) -14.0 E39-E35 S31-S29
      5 1249  706 (S16W20 L= 79)   277     88        1 (1203/ 685-1291/ 735)  31.0 W17-W22 S17-S13
      6  880  677 (S25E02 L= 57)   327    127        2 ( 808/ 656- 942/ 691)   6.3 E07-W02 S26-S23
      7 1730  521 (S17W65 L=124)  1531    214        4 (1661/ 436-1783/ 591)  67.4 W61-W69 S23-S11
      8 1315  473 (S29W32 L= 92)  2499    622        5 (1067/ 380-1612/ 541)  27.8 W18-W55 S40-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