Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1444 1712 (N31W55 L= 65)   461     89        0 (1414/1684-1476/1745)  80.5 W53-W57 N28-N34
      2  604 1618 (N53E14 L=355)  2130    401        6 ( 425/1528- 769/1692)  -1.5 E35-W03 N50-N55
      3 1137 1450 (N26W19 L= 28)   339    111        0 (1109/1404-1161/1487) -90.0 W17-W20 N22-N28
      4  444 1346 (N38E32 L=338)   376     85        2 ( 399/1325- 489/1360) -18.4 E36-E27 N36-N40
      5 1205 1262 (N13W16 L= 26)   479    117        2 (1161/1232-1264/1296) -45.0 W14-W19 N10-N15
      6  222  908 (N16E59 L=311)  5654    734       18 ( 115/ 699- 381/1096) -87.6 E74-E43 N04-N28
      7 1900  997 (S22W63 L= 73)   629     94        1 (1861/ 974-1937/1024) -21.0 W57-W70 S23-S18
      8 1152  752 (S14W00 L= 10)   771    310        1 (1072/ 628-1231/ 871)  35.5 E07-W07 S18-S08
      9 1344  661 (S24W09 L= 19)  1043    291        3 (1216/ 573-1417/ 757)  18.4 W00-W15 S27-S19
     10  899  615 (S15E18 L=352)   397    133        0 ( 847/ 578- 953/ 660) -66.8 E19-E15 S18-S11
     11  580  300 (S28E57 L=313)   672    103        1 ( 555/ 251- 613/ 337)  18.4 E65-E50 S30-S25
  

    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