Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  381 1414 (N29E44 L= 17)  3760    444        6 ( 269/1246- 480/1600) -31.0 E78-E33 N17-N45
      2 1550 1476 (N22W41 L=102)   755    158        0 (1471/1449-1623/1497) -15.3 W35-W47 N19-N23
      3  168 1223 (N19E65 L=355)  3959    368        7 (  84/1071- 223/1373)  67.2 E81-E58 N10-N29
      4 1896 1300 (N08W73 L=134)  2006     88        2 (1861/1254-1930/1347) -21.8 W68-W80 N05-N11
      5  419  699 (S15E45 L= 16)   844    246        0 ( 355/ 608- 484/ 795)  32.0 E53-E37 S19-S09
      6 1078  661 (S24E00 L= 61)   709    202        1 (1001/ 633-1175/ 710) -26.6 E04-W06 S25-S19
      7  232  588 (S20E71 L=350)  2149    115        2 ( 189/ 534- 267/ 638)  74.1 E79-E65 S23-S16
      8 1480  472 (S41W31 L= 92)   349     50        0 (1463/ 447-1496/ 504) -53.1 W28-W33 S42-S38
      9 1431  453 (S42W27 L= 87)   391     54        2 (1413/ 418-1455/ 483) -78.7 W24-W29 S44-S39
     10 1271  418 (S43W13 L= 74)   279     51        2 (1244/ 393-1292/ 443) -76.0 W10-W14 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