Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  514 1558 (N31E46 L=180)   995    138        1 ( 459/1507- 570/1604) -12.1 E54-E40 N27-N33
      2 1762 1463 (N35W63 L=289)   728     81        1 (1731/1430-1790/1493) -26.6 W59-W67 N33-N37
      3  191  843 (S17E65 L=162)  1873    155        4 ( 157/ 772- 232/ 895)  63.4 E70-E58 S21-S13
      4  387  820 (S16E43 L=184)   370     59        1 ( 367/ 793- 419/ 844)  76.0 E44-E39 S18-S14
      5  755  795 (S15E15 L=211)   537    173        1 ( 667/ 769- 830/ 817)  -5.2 E21-E10 S15-S13
      6 1289  672 (S18W20 L=247)   307    104        0 (1270/ 623-1315/ 722)  90.0 W19-W21 S21-S14
      7 1415  667 (S17W29 L=256)   282     68        1 (1377/ 645-1446/ 685) -11.3 W26-W31 S18-S16
      8  908  599 (S26E04 L=222)   246     61        0 ( 876/ 585- 948/ 613)   0.0 E06-E01 S27-S25
      9 1514  495 (S28W42 L=269)   368     42        2 (1489/ 479-1541/ 515)  26.6 W39-W44 S29-S26
     10 1407  435 (S34W36 L=262)   452     55        0 (1374/ 418-1441/ 453)  -9.5 W32-W38 S34-S32
     11 1069  427 (S37W09 L=236)   274     53        0 (1035/ 405-1096/ 445)  53.1 W07-W11 S39-S35
     12  488  416 (S45E43 L=183)   531     64        0 ( 461/ 388- 518/ 439)  24.4 E49-E38 S47-S42
  

    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