Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1089 1785 (N44E25 L=184)   525     90        2 (1059/1751-1125/1813)  90.0 E26-E22 N40-N47
      2 1678 1581 (N53W36 L=245)   499     53        1 (1653/1565-1709/1599)  32.5 W31-W42 N49-N56
      3 1836 1338 (N40W54 L=263)   775    106        2 (1785/1311-1879/1363)   7.1 W46-W62 N37-N41
      4  760  609 (S37E04 L=205)   349     69        1 ( 742/ 573- 779/ 643)  90.0 E05-E02 S39-S34
      5 1077  615 (S27W16 L=225)   292     95        0 (1030/ 591-1125/ 639)  59.0 W13-W18 S29-S23
      6 1440  565 (S18W40 L=249)   296    100        1 (1428/ 523-1455/ 621) -45.0 W37-W42 S20-S15
      7  792  514 (S42W02 L=211)   347     56        0 ( 767/ 488- 826/ 535)  68.2 W00-W03 S44-S39
      8 1557  432 (S21W58 L=266)   613    102        1 (1500/ 410-1608/ 453)  32.0 W53-W62 S23-S18
      9 1214  331 (S39W41 L=250)   543     81        1 (1179/ 302-1248/ 361)  90.0 W37-W43 S41-S36
     10 1339  262 (S38W57 L=266)  1037    160        2 (1270/ 215-1402/ 306) -78.7 W54-W60 S43-S33
     11 1145  241 (S47W46 L=255)   340     32        1 (1126/ 222-1162/ 267) -26.6 W42-W48 S48-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