Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  682 1553 (N48E08 L= 43)   354     58        0 ( 657/1526- 711/1576)  18.4 E10-E04 N46-N49
      2  346  897 (N15E47 L=  4)  2706    428        6 ( 261/ 760- 480/1081) -62.4 E53-E39 N04-N27
      3  709  996 (N13E19 L= 32)   245     78        0 ( 665/ 975- 736/1011) -26.6 E21-E17 N11-N14
      4  646  967 (N13E24 L= 27)   311     90        2 ( 622/ 933- 673/1006) -78.7 E25-E22 N10-N15
      5  677  876 (N07E24 L= 27)   227     68        2 ( 640/ 846- 698/ 901) -56.3 E25-E22 N05-N08
      6 1365  791 (S15W13 L= 64)   204     77        1 (1343/ 761-1382/ 832)  45.0 W11-W14 S16-S13
      7 1833  758 (S34W53 L=104)   487     45        2 (1815/ 733-1853/ 783) -76.0 W49-W55 S35-S31
      8  523  435 (S16E51 L=  0)   975    208        2 ( 465/ 356- 588/ 511)  85.6 E53-E48 S22-S09
      9 1072  307 (S39E21 L= 30)  3286    630        5 ( 896/ 115-1260/ 508)  21.1 E58-E01 S51-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