Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  696 1043 (N16E19 L=  5)  2269    509        4 ( 566/ 910- 875/1233) -65.2 E23-E10 N04-N30
      2 1052 1186 (N15W06 L= 29)   730    273        4 ( 924/1127-1150/1228) -45.0 W00-W10 N09-N20
      3  367  988 (N20E43 L=340)   462    163        1 ( 317/ 957- 435/1037) -26.6 E47-E37 N18-N24
      4 1766  825 (S27W44 L= 67)   490     56        1 (1730/ 801-1795/ 849)  -8.1 W40-W47 S28-S25
      5  377  763 (N07E48 L=335)   639    155        2 ( 322/ 724- 418/ 825)  16.7 E53-E43 N05-N09
      6 1527  572 (S35W20 L= 43)  1750    326        3 (1392/ 489-1616/ 677)   8.5 W08-W28 S38-S29
      7 1255  323 (S44E08 L= 16)   684    105        0 (1222/ 286-1292/ 371)  32.0 E11-E03 S46-S40
      8 1184  257 (S47E16 L=  7)   500     55        1 (1164/ 229-1211/ 288)  63.4 E19-E13 S49-S43
      9 1117  224 (S48E24 L=360)   392     51        0 (1093/ 204-1145/ 245) -45.0 E25-E21 S50-S45
  

    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