Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  591 1361 (N36E24 L= 90)   419    122        1 ( 536/1325- 640/1391) -29.1 E28-E19 N32-N38
      2 1858 1013 (S13W60 L=174)  2750    407        5 (1808/ 899-1903/1122)  72.9 W54-W67 S20-S07
      3 1775 1030 (S10W51 L=165)   972    158        5 (1735/ 975-1815/1078) -76.0 W46-W54 S14-S06
      4 1073  781 (S08E02 L=112)   509    195        0 (1046/ 702-1116/ 872)  56.3 E04-W02 S12-S03
      5 1616  815 (S18W34 L=148)   271     43        2 (1594/ 787-1632/ 843) -76.0 W32-W35 S20-S16
      6 1248  721 (S16W07 L=122)   631    277        1 (1198/ 619-1340/ 770)  32.5 W03-W14 S21-S12
      7 1612  477 (S43W34 L=149)   435     41        1 (1581/ 460-1638/ 491) -23.2 W30-W37 S45-S41
      8  565  281 (S34E59 L= 55)  1048    116        0 ( 539/ 231- 595/ 331)  19.1 E76-E50 S38-S29
  

    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