Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1685 1440 (N20W55 L=105)  1528    390        3 (1564/1317-1749/1573) -80.5 W50-W60 N12-N30
      2  725 1374 (N32E16 L= 34)   242     85        0 ( 686/1363- 769/1384)  -9.5 E18-E12 N31-N32
      3  405 1291 (N30E41 L=  9)  1907    399        4 ( 269/1234- 546/1338)   4.8 E53-E29 N26-N31
      4 1338 1255 (N15W23 L= 73)   224     89        1 (1310/1210-1358/1299) -80.5 W22-W24 N12-N18
      5  869 1142 (N15E08 L= 42)   741    224        3 ( 760/1080- 944/1222) -47.7 E14-E03 N10-N21
      6  326  650 (S10E56 L=354)  3261    581        7 ( 265/ 495- 364/ 817)  49.4 E71-E49 S20-N01
      7  498  688 (S08E39 L= 11)  1310    348        4 ( 453/ 572- 594/ 757)  50.7 E44-E32 S15-S04
      8  456  375 (S30E56 L=354)  1944    411        3 ( 388/ 329- 536/ 451)  16.7 E65-E45 S34-S24
      9 1146  335 (S42E03 L= 47)   306     56        0 (1126/ 311-1162/ 359)  45.0 E04-E01 S44-S40
  

    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